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ABSTRACT 


A theory has been proposed for conducring model experimentation ' 
i-fhich .permits the use of imperfect models or test conditions that 
do not conform strictly to similitude requirements. It is based on 
a new concept of representing the error states in a multi-dimensional 
Euclidean space when errors in the modeling parameters are small. 

A consideration of the error path in the hyperspace leads to expres- 
sions from which the global effect of these errors may be evaluated 
approximately. Conditions under whicli these expressions would yield 
satisfactory results are given. 

To test the usefulness of the theory, a computer experiment has 
been carried out for the prediction of the transient and steady state 
thermal behavior of a hypothetical spacecraft using perfect, as well 
as imperfect, models. The hypothetical spacecraft has major radia- 
tive and conductive heat flow paths that crudely simulate those of 
the *64 Mariner family of space vehicles. Errors in the conductance 
and capacitance parameters of up to 25 percent and in surface emit- 
tances of up to 20 percent existed in the imperfect models. Exten- 
sive temperature data vfere obtained for the various components of 
the spacecraft v/hen it was subjected to a sudden heating, due to its 
.exposure to the sun's radiation and the povrer dissipation vjithin the 
bus, followed by its attaining the steady state and the subsequent 
cooling as a result of removing the sun's heat. Data were also ac- 
quired under the condition of a cyclically varying thermal environment. 
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Generally speaking, errors of up to 10, 20, and 30°R, originally pres- 
ent in the data of the imperfect models,' were reduced to less than 
3*^R after correction according to the proposed theory. 
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NOMENCLATURE 


ft fr-,2 

A = area, ft 

, .. -a . -=--<;oefficient-“mat-rix -defined in--(3.2-,15) 
i i 

2 

B = diffuse radiosity, Btu/hr-ft 

b = inverse of a , defined in (3.2.17) 

i j • • - • --i i ■ 

3 

C = volumetric heat capacity, Btu/ft -®R 

C = weighting factor for 6, defined in (2.1.5) 
i 

d = thickness, ft 

E = exchange factor 

F = shape factor 

H = irradiation, Btu/hr-ft^ 

k = thermal conductivity, Btu/hr-ft -°R 

Q = heat flow rate, Btu/hr 

q = heat flux, Btu/hr-ft^ 

2 

5 = solar constant, Btu/hr-ft 

s _= .arc length along error path 

T = absolute temperature ; = absolute reference temperature 

t = time, hr 

a = surface absorptance 

a. = error ratio defined in (2.1.9) 

] 

.^3. = error ratio defined in (2,1,14) 

6 = error in modeling parameter 

6 = Kronecker delta function 

i j 

e = surface emittance 
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p = surface reflectance; = diffuse component of surface reflec- 
tance 5 = specular component of surface reflectance 

-8 2 4 

O = Stefan-Boltzmann radiation constant = 0.1713 X 10 Btu/hr-ft -°R 
= diffuse-specular overall radiant exchange factor 

Subscripts 

e = refers to empty space 
» , i = refers to surfaces A , A 

1 ’ j 

m = refers to model 

P = refers to prototype 

Superscripts 

(+) 

( -) 

(Others not found in the list are defined in the text.) 


= refers to solar spectriun: B*, E*, H*, Q*, q*; a*, p|, and p*; 
also refers to modification experiments: , 

and 

= refers to subspace 


= refers to sub space S 
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Chapter 1 


MOTIVATION OF A NEW APPROACH TO THERMAL SCALE MODELING 

Theoretical requirements for the thermal scale modeling of un- 
manned spacecrafts are understood quite well at the present time. 
Difficulties arise in practice when the system to be modeled consists 
of con^onents fabricated from a number of materials and they all par- 
ticipate in influencing the system's thermal performance. This is 
particularly true when information on both the steady state and tran- 
sient behavior is sought. One possible means of satisfying modeling 
requirements is to modify artificially the thermal conductance of 
thin struts, plates, shells, etc., by electroplating vjith a high con- 
ductivity metal, such as copper. The main objective of a previous 
research contract betv:een the University of Illinois and the Jet Pro- 
pulsion Laboratory, JPL No. 951660, was to develop such a possibility. 
Test results obtained with prototypes and models of simple configu- 
rations have demonstrated the technical feasibility of the concept, 
either with or without solar simulation. With the availability of 
adhesive metal foil tapesf in recent years, it is conceivable that 
they can be applied more conveniently than electroplating. In this 
connection , ‘ one also sees the possibility of modifying the materials' 

\ 

heat capacity by using plastic (teflon) tapes. Hovfever, the range 


tone type of copper and several aluminum foil tapes with pressure 
sensitive silicone adhesive, manufactured by Mystik Tape, Inc., of 
Northfield, Illinois, may be suitable for the present purpose. 



of thickness of commercially available tapes is limited. Consequently ,• 
the desired thermal conductance and capacitance of the various struc- 
tural members of the system may not be menufactirred precisely in this 
manner. They can only be approximated. Furthermore, the application 
of a metal foil tape to a surface would also introduce a minor dis- • 
turbance to its heat capacity. Likewise, the use of a plastic tape 
— would produce a minor influence on its conductance. It is also pos- 
sible that the desired surface radiation properties may not be dupli- 
cated accurately. These considerations prompted us to explore and 
examine another concept of model testing in which rnodels that do not 
satisfy the similarity criteria completely are used. The theory , 
which is presented in the following chapter, contains some heuristic, 
but plausible, arguments and, in this sense, it is not mathematically 
rigorous. It was described first, in its rudimentary form, in a pre- 
liminary proposal, "Thermal Scale Modeling with Imperfect Models," 
prepared by the senior author and submitted to JPL for technical evalua 
tion in January 1969. The material contained in Chapter 2 is a gen- 
eralization of the theory presented in that document. It includes 
additional information relating to its practical implementation. 

To test the validity of the theory, it has been applied to a 

hypothetical spacecraft having a geometric configuration that crudely 
\ 

simulates the global radiative and conductive paths of the Mariner 
spacecraft. The results of the present study firmly established the 
usefulness of the proposed new concept and one sees the v/elcomed flexi- 
bility in the use of model testing — a flexibility heretofore impos- 
sible to attain. 
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Chapter 2 


AN APPROXIMATE THEORY OF MODELING WITH IMPERFECT MODELS 

To begin with, it should be noted that, v?hile our concern here 
is the thermal modeling of unmanned spacecrafts , the concept to be 
expounded herein has general applicability to the immense field of 
the technology of model testing. As such, the analysis and discus- 
sion that follow are presented with this viewpoint. 

2.1 BASIC IDEAS 

Consider a physical phenomenon y vihich depends on a number of 
independent variables , • • * > • Symbolically, we may write 

y = y(x^ , x^ , • • • 5 5<p ) (2.1.1) 

If the functional relationship has general validity, it must be in- 
dependent of the units of measurement. Thus, a consideration of the 
requirement of dimensional homogeneity leads to a dimensionless form 

I 

of the relation, 

^ n^, • • •, n^X ' (2.1.2) 

Methods of dimensional analysis, whether based on the governing dif- 
ferential equations of the problem or on the strict algebraic formal- 
ism or on physical intuition, have been vrell documented and need no 
further deliberation here. Suffice it to state that, for perfect 
modeling, the H’s in the model system are made indentical to the cor- 
responding H's in the prototype system. 


3 



In the modeling of the thermal performance of spacecrafts, the 
dependent variables of interest are usually the temperatures of -their 
various components. Thus, the number of (|)'s viith which the experi- 
menter has to deal -could be quite large. Moreover, as has been noted 
previously, the precise satisfaction of all modeling requirements 
may not be 'possible in practice. In fact, even if it -were possible, 
it might not be desirable from the point of view of economy- This 
is, then, the fundamental reason for using imperfect models and some 
deviation in the II 's of the model system from those of the prototype 
will be tolerated. To facilitate further discussion, we rewrite (2.1.2) 
as 



= (n, , n 


n ; fi) 

n 


(2.1.3) 


where 1=1, 2, • • 'jNjN being the number of different components 
of the spacecraft to which a characteristic temperature is assigna- 
ble. The subset of all relevant variables or parameters for which 
the model system is free from error is denoted by fi. Clearly, for 
the problem under consideration all elements of are passive and 
can be Ignored. In the transient thermal modeling of unmanned space- 
crafts, n is usually several times the value of N, as will become 
clear later. 

\ 

If v?e designate the errors in H by H6 and define 




H (1 + 6 ); (2.1.4) 

n n 


then may represent the measured entity using the Imperfect model ; 
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e.g. , the temperature of the component surface A. of the model space- 
craft which does not satisfy the modeling requirements completely. 

The quantity sought is, of course, . V?hen the 5's are sufficiently 
small, the right-hand side of (2.1.4) can be linearized to become 



(2.1.5) 


in which the coefficients C. . are not known. It may be noted that, 
if Taylor's series is used, then 

C. . = H. - (2.1.6) 

I , ] ) 3II. 

J 

However, (2.1.6) does not necessarily give the best estimate of the 
coefficient C. . . In fact, the validity of (2.1.5) does not require 
the usual restriction that needs to be imposed on the function for 
the existence of the Taylor expansion. 

To evaluate from the measured , one would require theo- 

retically the knov7ledge of the n unknovm coefficients Cj . . This 
means that n experiments will have to be conducted with n different 
sets of 6’s. The work may become formidable even when n is only mod- 
erately larger say, 20 or so. Furthermore, the accuracy required 
of the experimental measurements under such circumstance could be 
so high that they become impractical. We, therefore, seek alterna- 
tive solutions with some loss of accuracy. 
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Vfe begin by noting that the 6 ' s associated with the imperfect 

model may be either positive or negative. For reasons to be seen 

later, they are grouped separately. We designate that > ^2 ’ * * * - 

are positive and that 6 ,6 3 * * * j 5 are negative. Purther- 

r + 1 r + 2 n 

more, we define 



(2,1.7a) 


(2.1.7b) 


Clearly, and ^ may be either positive or negative. With the 

i i 

foregoing, (2.1.5) may be rewritten as 

^ = 4 )^ ( 2 . 1 . 8 ) 

i i i i 


If the imperfect model is modified or if the model test condition 
is altered (in either case, the extent of modification must be restricted 
so as not to violate seriously the linear approximation used in the 
theory), then A^ and A^ ^ will change also. A possible practical 
means of effecting model modification is through the use of adhesive 
metal and plastic tapes for altering the thermal conductance and heat 

N. 

capacity as indicated in [l]t. Thus, originally positive errors will 
become more positive, while originally negative errors vjill become 
less negative or positive . In what follows , the errors following 

tNumbers in brackets refer to entries in REFERENCES. 
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modification will be designated with an asterisk. 

Let us consider that the modifications are carried out in two 
steps. In the first, only the positive errors are altered, leaving 
the negative errors intact or essentially unchanged. We denote the 
error ratio by a. Thus, 

6*/6 = a , j = 1, 2, • • • , r . (2.1,9) 

) i. ] 

If <!)/^ represents the experimental data sc obtained, then 


n 






> c .a.' 6 . + >■ C. . 6 . (2.1.10) 

/ f 1 ,] • J 1 X , 1 .) } 


j = 1 


j = r + 1 


Eliminating (j/ from (2.1.5) and (2.1.10) yields 


X' 

^ = y C . 6. (ct - 1) 

i 1 / I 1 , J ] ] 

j = 1 


( 2 . 1 . 11 ) 


If a weighted average of the error ratio is defined according to 


E C 6. a 

_ • TL i.) J 3 

rw A ^ — 


a = 


r 

Ec .6. 

-n 1.3 ] 


(2.1.12) 


j=l 


then (2.1.11) can be rewritten as 



- = (a - 1) 



j = 1 


C. . 6. 

1.3 } 
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and, consequently, there is obtained from (2.1.7a) 


d>‘‘ 

.(+) - 


...( +) 




(2.1.13) 


a - 1 


Obviously, a cannot be evaluated without the prior knowledge of C. . . 

^ I 3 

However, in the event that all a ’s are identical, say, a, then a = a 
.and.no knowledge of C. . is needed, ^ftiile such special condition 
is usually not met in practice , the foregoing observation does sug- 
gest that, in testing with imperfect models, the experimenter should 
seek small deviations in a 's. 

In the second step' of modification, only the negative errors 
are altered. For clarity and for convenience of later discussion, 
the error ratio is to be denoted by 0 . Thus , 


Sr/S. =3. , j=r + l, rt2, 

3 ] ] 


n 


(2.1.14) 


In this instance, the measured data, ^ , can be written as 


n 


<}>:■ 


- cp. + 


Yc. .6, . y 


C. .0 6. 

J I 1 3 3 


j = 1 


j = r + 1 


(2.1.15) 


It follov7s, then. 


n 


^ J > C. . 6. (3 

^ / > » . 3 ] 3 


1 ) 


j = r t 1 


(2.1.16) 


As in the first case, if a weighted average of the error ratio is 
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■ defined j namely. 


n 


3 = 


2 c. .6.3. 


n 




C. . 6. 

> I 3 3 


then 




3-1 


(2'.1,17) 


(2.1.18) 


The observation previously noted for a can likewise be made for 3* 
Thus, without the knowledge of C. . 's, the precise evaluation of both 

^ I 3 

a and 3 is not possible. However, they can be approximated. It is 
this desire of finding the appropriate mean values of the error ra- 
tios that led to the creation of a new concept of model testing with 
imperfect models . 


2.2 GEOMETRIC REPRESENTATION OF ERROR STATES 

I 

The central idea to be explored herein is the representation 

of error states in a multidimensional Euclidean space whose elements 

1 

are 6^ . This will be referred to as the simple error space in con- 
trast to fhe. weighted error space to be discussed later. Such rep- 
resentation is compatible v;ith the lineal- theory considered in the 
present study. 

For j ranging from 1 to n, there exist the following three points 
in the n-dimensional error space: 
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1 ^ ^ ’ i + l ’ ^ n ^ 

Fd^) (5*, 5*, . . ., 6^:, . . ., 6J 

P--^“^ (6^ , 6^, • • ♦, 6 , 6;\^ , ‘ 6*) 

There ere also two subspaces and ^ ^ vjlth their respective ori- 
gin at and (i~^ . The elements of subspace ^ are <5^ = = ' * * 

and those of d ^ are 6 ,6 , • • • , 6 . If the etror ratio a 

r + l ’ r+2 ’ ’ n 

is a constant, the three points 0 , P, and P* are colinear in 

. The straight line which originates at and passes through 


the said points will be labeled ^ 


Similarly, if 3 is a constant. 


0^”^ , P, and P*^“^ are colinear in ^ and the straight line join- 
ing them’ vrill be labeled d ^ • When n = 3 , the foregoing can be rep- 
resented pictorially as is illustrated in Fig. 2.1. No such picto- 
rial representation is possible when n > 3 . 

In , the lengths of the line segments O^^^P and 


0^“^^ P = 


j = 1 


and 



If we regard ss the abscissa of a plane plot with (|). as its ordinate. 
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( +) (— ) 

Figure 2,1 Illustration of subspace S and S in a 3-dimensiona3 
simple error space. ( 6 ^ , 62 += ^3 is is a 

plane while is a line). 





then a point corresponding to P and another point correspond- 
ing to may be located as shown in Fig. 2.2(a). A conjecture 

introduced in [1] is that (j). depends uri'^uely on and varies line- 

arly with it. It follows then 


a(+) - 


- S 

i i 


cjC+) 


a - 1 


p 


( 2 . 2 . 1 ) 


which is identical to (2.1.13) when a = a, as is the case here. 

An analysis based on analogous arguments can be made in subspace 
. Thus, ' 


0^ ) P = 


2 < 

.j = r + 1 


1/2 


0(-)p*(“} 


n 


(6*)^ 

] 


Lj = r + 1 


1 / 2 


3 0^ ^ p 


and, from the vs ^ ^ plot as illustrated in Fig. 2.2(b), it can 
be shovm readily that 


¥ 


) _ I 


q( ) pv’{( ) 


3-1 


( 2 . 2 . 2 ) 


0^ 


which is .precisely (2.1.18) when 3 = &• 
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The foregoing analysis , plausible as it may seem , is subject 

to immediate criticism. Referring to (2.1.5), one recognizes that 

the deviation of A from A depends on the products C. . 6. instead 
1 .1 I , ] J 

of 6^ alone. The coefficients C. . are, in essence, Vfeighting fac- 
tors which prescribe the relative importance of the various errors 
in influencing the dependent variable. Thus, a we-ighted error space 
-whose elements-are ] C. -should be used instead of the simple er- 

ror space and it is more logical to assume that the distance ascer- 
tained in such error space from its origin is an appropriate measure 

of the resultant error in d). . 

1 

. A. quest ion which naturally arises is: V/hy is it, then, that 

for constant a and 3? the analysis made in the sample error space 
leads to correct results? The reason for this may be seen easily. 

In of the weighted error space , the lengths of the line segments 


0^"^^ P and are 


0^'*’^ P = 


r 


-, 1/2 


2 

Lj = 1 


(c. . 6. r 

1 .] 3 


and 


-, 1/2 




V (c. . 6;V 

/Lj 1 1 3 3 

-j = 1 


Hence, v?hen a is a constant, 0 
0^*^ , P, and are co linear. 


P"^ ^ = a 0^ ^ P and the three points 
Since only the vat'io of distances 
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is involved in the expression for , it is seen that (2.2.1) re- 
mains unaltered. For the same reason, (2.2.2) also holds. It is 
significant that the linear correction formulae (2.2.1) and (2.2.2), 
valid for constant a, and g and originally deduced from a conj ecture 
in the simple error space, are formally identical to those based on 
a consideration in the weighted error space. It goes without saying 
that the experimenter has normally no knowledge of C. . and, hence, 
the error state- points like P, , and P*^ ^ cannot actually be 

plotted in the vieighted error space. 

2.3 GLOBAL EFFECT OF POSITIVE AND NEGATIVE 6'S ON (j). —PARABOLIC ER- 
ROR PATH. WORKING FORMULAE. 

We now proceed to examine the more realistic case for which the 

error ratios a, and g are not constants. VJe restrict ourselves to 

j ] 

the consideration that the range of variation among the 's and among 

the g. 's is limited. The analysis which follows closely parallels 

that described in [1] for the evaluation of the global effect of the 

errors based on a consideration made in the simple error space. 

Referring to Fig. 2.3(a.l), the line joining 0^^^ and is 

designated as the -axis. Since a. is not a constant, but varies 

( + ) 

with j 5 P does not lie in this line . Hence , the point P and the ^ 

^axis together determine a plane in S of the weighted error space. 

In this plane , a line 0^ ^ ^ ri^ is erected normal to the ^ -axis 

{ + ) 

as shovm. For convenience, the ri -axis is so oriented that P lies 
on the positive side. It is pertinent that ? and R imply 
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and TX ' since tlie scale of the coordinates in the weighted error 

1 

space' depends on i. For simplicity, the subscript i is omitted-. 

Passing through the three points 0^^^ . P, and p*^^^ , a parabola 
can be constructed. Tracing the points along the parabola and toiirard 
the origin 0 may be interpreted as conducting a series of experi- 
ments with controlled errors of successively diminishing magnitude. 

As the origin is approached, all positive errors become vanishingly 
small. Thus, the are length along the cur>ve measured from is 
an appropriate measure of the global effect of all positive errors 
in . The choice of a parabola as the error path is, to some extent, 
arbitrary, but it is certainly reasonable. Further argument in sup- 
port of the choice will be given later. 

The equation of the parabola in the ^ -Ti plane is 

( 2 . 

where a > 0 and b < 0. The coordinates of are 


'=’p* 


‘p* 


-, 1 / 2 


2 

Lj = 1 


(C. 5*) 

} J 


= 0 


( 2 . 


( 2 . 


'In (2.3.2a), C implies C. . The coordinates of P are 

J 1 I 3 







■*•5 p::/ 


2 

> C 6 6* 
j=l ^ 


f, 

Lj=i ’ ’ J 

1/2 


( 2 . 


3.1) 

3.2a) 

3.2b) 

3.3a) 
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^■*1. = [id'-j p|=’ - c4,*> f]'” 


- j 


r- 

= Z 'S«1 

-j = 1 


(2,3.3b) 


in which the overscore arrow denotes vectors . It can be shown that 


■ 1 ■ , 
a = , D 

■ c 


-N. 

The arc length 0^ ^ ^ ^ ^ is 


r 

!(*>=[ j 

4) 




, ( +) >, 2 1/2 

ds'*’ 


{a(l + + Jin [(1 + a ) ' + a]} (2.3,5a) 


and the arc length 0^^^ P is 


r(+) 

. [‘-(Nj 


- ^{a(l + a^)^^^ + Jin [(1 + a^V'^ + a] 

- aid + - £n [(l.-i- + aA]} (2.3,5b) 


2 . 1/2 


where 


A = 1 + 


2b £.(+) 


(-1 < A < 1) 


(2.3.5c) 
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The ratio of the two arc lengths is 



2 

1 - f(a,X) 


(2.3.6a) 


where 


f (a,X )■ 


, 2 , 2 . 1/2 „ r /-, 2 , 2 . 1/2 

aX (l + a X ) + £n [(1 -r a X ) + aX] 

2 . 1/2 . r /-. 2 . 1/2 

a(l + a ) + Jin [(1 + a ) + a] 


(2.3.6b) 


and a is given by (2.3.4a) and X by (2.3.5c). In the latter two 
equations , we further note that 


{+) 

V (C 6. f 

/ ! J 3 

n 

1/2 



_ 




(2.3.7a) 






C. 6. 6* 

j ] ) 




(6^0' 


(2.3.7b) 


In the foregoing, the summand ^ implies ^ . Consequently, within 

j=l 

the context of the present theory, the correction required for all 
.positive errors is 






- 1 


(2.3.8) 
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. — ■ .(•<)(+) 
and the appropriate mean a is seen to be given by the ratio s“ /s 

• For the practical application of these results , the following 

observations are made . 


(A) Vrhen the range of is restricted, (2.3.7a) and (2.3.7b) 
are weak functions of C . They may be replaced by the following ap- _ 
proximate expressions. 


C’ 

£■( 



(2.3.9a) 


(2.3.9b) 


Alternatively, if the physical nature of the problem is such that 
the dominant C. ’s are of similar magnitude, our experience indicates 
that these approximations are also valid. 

Following a detailed error analysis based on the same arguments, 
but made for the simple error space, it was found that (2.3.6a,b), 
(2.3.4a), and (2.3.5c) remained unaltered provided that, in the latter 
^two expressions 5 the ratios T]^ and ^ were calculated 

according to (2.3.9a) and (2.3.9b). 

(B) If the experimental conditions are so chosen that, in the 
weighted error space, 2(Cp^^ ) = 1, then X = f(a,X) = 0 and 
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s* /s =2 which is independent of the 's. This condition is, 

of course, not generally met in practice. However, if 'a' is small, 
say, less than 0.5, then f(a,X) becomes somewhat insensitive to varia- 
tions in 'a'.' Furthermore, if we require that [Xj be kept small, say, 
vzithin 0.1 or 0.2, then f(a,X) is small compared to 1. Under these 
conditions, the error in the calculated ratio of the arc lengths re- 
sulting from inaccuracies in 'a' and X will likewise be small. Thus, 
the use of the parabolic error path provides guidance to the proper 
selection of the experimental conditions for reducing errors due to 
the lack of knowledge of the weighting factors U . 

The foregoing observations (A) and (B), taken together, suggest 
that, in the application of the linear correction formula (2.3.8), 
the ratio of the arc lengths (or its equivalent a) can be evaluated 
from a consideration of the simple error -space when the range of a 
is restricted and/or when the experimental conditions are chosen as 
explained . 

If, in the modified experiments, the positive errors are made 
less positive, the relative location of P and with respect to 

the origin of S would be as sho^-m in Fig. 2.3(a.2). The ^ -axis 
is now chosen to pass through P. Following a similar analysis, it 
;;as found that 





.•,( +) 


(+) 


(2,3.10a) 
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in v?hich 


s''^ _ 1 - f(a,X) 

(+) 2 
s 


(2.3.10b) 


and f(a,X) is given by (2.3.6b) with 'a' and A redefined as follows 


wherein 


and 


(+) 

^P* 


A = 1 - 2 


-{+) 

^p* 

4+) 

=>p 




rS*'> 


/ I ] ] / yr.. .i } } 


E 


C^6. 6* 

) ] j 


-tl/ 2 


- 1 


+ ) 
* 


4" 


yc^6. 6:^ 

Z—j } } i 


2 .2 
C. 6. 

) ] 


(2.3.10c) 


(2.3. lOd) 


(2.3.10e) 


(2.3.10f) 


Again, if the range of a is restricted and/or if the experimental 
condition is such that |a| is small as compared to unity, the ratios 
4/ 'C' 4/ /4*’ can be calculated from a consideration of 
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the simple erroi? space j ■namely. 


A*) 


nl/2 




( ^ 6 . 6 * 
\ ^C-/ 3 3 


(2.3.11a) 


and 


it 


C’ 


^ 6. 5* 
X 1 1 

YT 


(2.3.11b) 


Clearly, in (2.3.10e,f) and (2.3.11a,b), the simimand 2! implies > 

j=l 

The analysis of the global effect of the negative 6's on 4>. can 
likewise be treated. If the negative errors are Z-ess negative in the 
modified experiments, the error path in ^ ■would be as shown in 
Fig. 2.3(b..l). We shall refrain from presenting the details of the 
calculation and, instead, shall summarize the results as follows. 


A. 

1 


(b. - i'' 
(-} _ J ^ 


*(-) 


1 - 


J- ) 


(2.3.12a) 


where the ratio of the arc lengths along the error path is 


^ 1 - f(a,X) 
„(-) 2 


(2.3.12b) 
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and f(a,A) is again given by (2,3.6b) with 'a' and X, respectively, 
denoting 


-wherein 


a = 


‘"p* 1 


r’ 


A = 1 - 2 


4 '.’ 

4 -’ 


4 ;’ 


4 


-) 

4 ‘ 


Z-J j 3 / Z__/ 3 3 


2 


|l/2 


- 1 


6 . 6 * 

3 3 3 


(2,3.12c) 


(2.3. 12d) 


(2.3.12e) 


and 




6* 

Z.J 3 3 3 


4 '^ ^ 


3 3 


(2.3.12f) 


In the above, ^ implies ^ ratio ^ /s^ ^ clearly denotes 

_ j=r+l 

'the mean p. 

I 

The error path would be as that depicted in Fig. 2,3(b.2) if the 
negative errors become more negative after modification. In this 
case , 
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where 


A - ) 


s 


2 

1 - f(a,A) 


(2.3.13b) 


and f(a,A) is as defined earlier, with 'a' and X, respectively, given 

by (2.3.4a) and, (2.3.5c) when the superscript (+) is replaced by (-). 

The ratios ^ ^ ^ ^ given by (2.3.7a) and (2.3.7b) 

n 

with the summand ^ reinterpreted to mean ^ . The discussion 

j=r+l 

on the validity of using approximations deducible from the simple 
error space also holds. 

Finally, the quantity sought is 



Chapter 3 


APPLICATION OF THE THEORY TO THE THERKiAL MODELING 
OF A HYPOTHETICAL SPACECRAFT-PRELIMINARIES 

3.1 THE HYPOTHETICAL SPACECRAFT 

To ascertain the usefulness of the approximate correction theory 
for imperfect modeling as expounded in the previous chapters we con- 
sider a hypothetical spacecraft shown in Fig. 3.1. The geometric 
configuration chosen is intended to simulate the global radiative 
and conductive paths of the Mariner ' 64 spacecraft but with all sec- 
ondary and minor details omitted. The solar panels and A^-A^ 

are fabricated of an aggregate of aluminum sheets , silicon crystals 
for energy conversion and glass plate covers , all mounted in an alumi- 
num allo^’' celled structure, VJhile the local thermal properties are 
of a highly anisotropic and inhomogeneous character, it is quite' 
feasible to describe its ovevatt behavior by equivalent values of 
thermal conductivity and heat capacity. Based on the information 
transmitted to us from JPLt , it appears that a reasonable value of 
kd for the solar panel is 0.058 Btu/hr-F. Also, the mean specific 
heat of the aggregate is reported to be 0.20 Btu/lb-F and a typical 
Mariner solar panel of 3 ft X 7 ft weighs about 30 lbs. From these, 
it is estimated that the Cd of the solar panel is approximately 0.29 
Btu/ft^-F. 


fThrough Hr. William A. Hageraeyer. 
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Collimated Sunlight 

(Either at Normal / Incidence or 45°) 



Solar Panel: 

Bus. ^3) ^7 5 1 ^9 3 

(A^and A^^ are covered with superinsulation) 
.Antenna: A,^ 

^ 2 

Ail areas are idenlica! and equal to 4 ft. 

Figure 3.1 Hypothetical spacecraft 
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To keep the radiation geometry associated v?ith the diffuse-specu- 
lar surfaces as simple as possible » the bus of the spacecraft is xaken 
to be a cubical box instead of the usual polygonal configuration'. The 
top and bottom faces of the bus are covered vfith multilayer superinsu- 
lation, The collimated sun’s rays are assumed to be either perpendicu- 
lar to the solar panels or slanted at 45° to them, but remain parallel 
to the front and back faces of the bus. TABLE 3.1 lists the materials 
selected for the six component faces of the spacecraft bus , The thick- 
nesses indicated are eqvA-Vatent values; they take into account the 
electronic packaging subassemblies bolted to them. Included is the 
relevant information for the main directional antenna. As V7ill be 
shox-jn in Section 3.4, perfect modeling requires that the corresponding 
surfaces of the model and the prototype have identical radiation proper- 
ties. Furthermore, v;hen the temperature fields are two-dimensional, 
the kd and Cd of all component surfaces of the model must bear a fixed 
relationship to those of the corresponding surfaces of the prototype. 
TABLE 3.2 summarizes al pertinent surface properties and the equiva- 
lent (kd)’s and (Cd)’s for the hypothetical spacecraft. 

3.2 RADIATION HEAT TRANSFER ANALYSIS AND THE DETERMINATION OF THE 
OVERALL EXCHANGE FACTOR ^ 

V/hen the surfaces of a spacecraft are under the irradiation of the 

\ 

sun, the absorbed and reflected radiation are in a wavelength region 
totally different from that of the emitted radiation. Vfhile it is pos- 
sible to perform a general formulation of the problem, there is neither 
the need nor the justification of doing such analysi.s at the present 
time. It is v/ell known that approximately 98 percent of sun’s radiation 
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TABLE 3.1 


MATERIAL LIST FOR THE HYPOTHETICAL SPACECRAFT 


(a) Solar Panel 

Aggregate of silicon, glass, aluminum and its alloy 

(b) Bus 


Panel 

Material 

k 

Btu/hr-ft-F 

C 3 

Btu/ft -F 

Equivalent 
Thickness 
d , i-n . 


Magnesium Alloy 
AN-M-29 

47.0 

25.9 

0.125 


Aluminum Alloy 
5086 

80.2 

37.1 

0..1875 

A 

7 

Titanium Alloy 
A-llO-AT 

4.60 

35.0 

0.10 


Magnesium Alloy 
• A781A(T4) 

29.1 

28.0 

0.1875 


Magnesium .^lloy 
2K51-A 

63.0 

27.6 

0.125 

> 

o 

Titanium Alloy 
•150-A(2) 

8.38 

38,5 

0.025 


(c) Antenna 

Aluminum alloy 2011 (k = 82.3 Btu/hr-ft-F and C ■= 39.3 Btu/ft^-F), 
Thickness - 0.0625 in. 
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TABLE 3.2 


SURFACE PROPERTIES AND OTHER DATA FOR THE HYPOTHETICAL SPACECRAFT 



Solar 

Infra-red 

kd 

(equivalent ) 

Cd 

(equivalent ) 

p* 

p* 

s 


Pd 

Ps 

a(=£) 



Solar 

Panel 

Sunlit Surfaces : A^ , , A^ , Ag 

0 

0.20 

0.80 

0.20 

’ 0 

0.80 

0.058 

0.29 

Surface Away from Sun: Ai' ,A 2 ' ,A 4 ',As' 

(painted with PV-100 white paint) 

0.80 

0 

0.20 • 

0.17 

0 

0.83 

Bus 

As 

exterior (PV-100 paint) 

0.80 

0 

0.20 

0.17 

0 

0.83 

0.489 

0.269 

interior (partially painted) 

— 

— 

— 

0.80 

0 

0.20 

A6 

exterior (PV-100 paint) 

0.80 

0 

0.20 

0.17 

0 

0.83 

l'.25 

0.580 

interior (partially painted) 

— 

— 

— 

0.80 

0 

0.20 

A, 

exterior (superinsulation) 

0 

0.85 

0.15 

0.50 

0 

0.50 

0.0383 

0.292 

interior (black paint) 

— 

— 

— 

0.15 

0 

0.85 

A 3 

exterior (PV-100 paint) 

0.80 

0 

0.20 

0.17 

0 

0.83 

0.455 

0.437 

interior (black paint) 

— 

— 

— 

0.15 

0 

0.85 

1 Ag 

exterior (polished aluminum) 

0 

0. 80 

0.20 

0 

0.95 

0.05 

0.655 

0.287 

interior (black paint) 


— 

— 

0.15 

0 

•0.85 

Alo 

exterior (superinsulation) 

0 

0.85 

0.15 

0.50 

0 

0.50 

0.0175 

0.080 

interior (black paint) 

— 

— 

— 

0.15 

0 

0.85 

Antenna 

All 

sunlit face (green paint) 

0.25 

0.05 

0.70 


0 

0.82 

0.428 

0.204 

face away from sun 

— 

— 

— 

■0.95 

0 

0.05 


All (kd)'s are in Btu/hr-F, (Cd)'s in Btu/ft^-F. 


































































IS contained within the 0-3 ym range.' On the other hand, less than 
0.01 percent of the radiant energy emitted by a black surface at 
530 R (room temperature) is in that wavelength range. Even at 750 Rt‘, 
less than 0.4 percent of the emitted radiation is. of wavelengths be- 
low 3 ym. The semigray or two-band analysis suggested by Bobco [-2] ' 
is based on the foregoing facts. The use of such a two-band model to 
account for the spectral dependence of the surface properties of space- 
crafts has been examined by Plamondon and Landram [3] and found to 
yield very good results . 

■ A simple and often realistic description of the directional prop- 
erties of the surfaces is a subdivision of the hemispherical reflec- 
tance into diffuse and specular components; i .e. , p = . In 

the analysis that follows, we shall adopt such an idealization which, 
together with the semigray approximation and the assumption of diffuse 
emission, makes possible an analytical treatment of the problem in 
a reasonably straightforward manner. The essence of such approach was 
used in a recent paper by Bering 

3.2.1 Radiant Rux 1n the Solar Spectrum 

Vfe first present a general formulation for an enclosure 
of N surfaces, each of which has a uniform temperature and is uniformly 
irradiated. The results v;ill then be applied to the relatively simple 
configuraticn selected for the hypothetical spacecraft . All proper- 
ties and radiant energies associated vath the solar spectrum will 
be designated with an asterisk. 

According to the two-band model, the. emitted radiation in the 

+This far exceeds the ordinary operating temperature range of present 
day spacecrafts. 
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solar spectrum from any surface of the system is completely negligible. 

Hence j the radiosity of a surface A. is given by 

B* = p* . (SE* + H*) (3.2.1) 

i d, 111 

in which E* is a type of exchange factor, denoting the fraction of 
the solar flux incident on A. directly and by all possible specular 
reflections . The irradiation H* is given bjf 



(3.2.2) 


where the exchange factor denotes the fraction of the diffuse 

radiation in the solar spectrum leaving A. which arrives at A. both 
directly and by all possible intervening specular reflections. VJe 
note that, in writing (3.2.2), the reciprocity relation A E* = A E* 

3 ] -1 1 1 -j 

has been used. Eliminating H* from (3.2.1) and (3.2.2) yields 


N 


B* = p* SE* + p* . 

1 d, 1 1 d, I 


E* . B* 

1 -3 3 


j = 1 


(3.2.3) 


which forms a system of N simultaneous, linear, algebraic equations 
\ 

for the N unknown radiositles. Once the B*'s are knovra, the rate 
of solar energy absorbed by A. can be calculated from 
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(3.2.4) 


Q* ■■ A, a/'(SE/' + H*) 

1 111 1 


ar 

A. -T^ B* 
^ Pr • ^ 

d, 1 


Obviously, B* and, consequently, Q/“ are independent of the system 
temperatures . The decoupling of the solar fluxes from those of the 
infrared radiation results in significant simplification of the analy- 


sis . 


If a surface A^^ reflects purely specularly, then = 0 and, 

accordingly , B* vanishes , In this case , 


(3.2.5) 


We now proceed to apply the foregoing results to the hypotheti- 
cal spacecraft when the sun’s rays are parallel to the front and rear 
faces of the bus and are either inclined at 45® or normal to the 
solar panels. Figure 3.2 shovjs the relevant geometries involved. 

To evaluate Q*, Q* and Q*, V7e need only to consider an enclosure 
consisting of A , A , A and an imaginary black surface A at ab- 

solute zero. The latter comprises a rectangle with its long sides 
shown as the dotted line in Fig. 3.2 and two identical right triangles, 
both adjoining the rectangle. Clearly, B* = 0. It is easy to 
see that 


Lf = cos 45®, - 0, - F^ 


(3.2.6a) 


= EL2 


E* 

3-3 


= 1-3 '= = 2-3 
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Designates superinsulation 



Top View 


Figure 3.2 Radiation geometry associated with hypothetical spacecraft 


Consequently 5 (3,2.3) becomes 


n - pf, ° 

P|, 3 <fs-l®f + '■3-2 ®l> - ®? = » 


(3.2.7a) 


When the radiosities are evaluated, Q*, Q* and Q* can be determined 
from (3. 2. 4) or (3.2.5). 

fr/hen the sun’s rays are normal to the solar panel. 


•s. 



Consequently , 


B* _ p* p B* = S 

1 ^d, 1 1 -3 3 1 

B* — p'* F B'> = p* S 

2 - ^d, 2 2-3 3 2 

p* (F B* + F E*) - B* 
^d, 3 ^ 3-11 3-2 2 , 3 






0 


J 


To evaluate Q* , Q* and Q* v?e consider a similar enclosure con- 
4 5 6 

sisting of , A^ , A^ and the associated imaginary black surface 
A^ at absolute zero. Again, referring to Fig. 3.2, we have 


(3.2.6b) 


(3.2.7b) 


= 

cos 

45° , 

^4-4 

= 0, 


= 0, 

= ^4-6 




(1 + 

P* A 

s, 6 

} COS 

45°, 

%-5 

= 0, 

-- ^5 -6 

► 

(3.2.8) 

= 

(1 + 

P* 

, 5 

) COS 

45° , 

^6-6 

= 0 
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and 


B* - p* F B" = p* -S cos 45® 

4 4 4-6 6 ^d, 4 

BiV _ p* p gj; - p-J: (1 + P* ' )S COS 45® 

5 ^d, 5 5 -6 6 . *^d, 5 *^8,6 

~ = -p* (1 + p* )S cos 45® 

'^d, 6 6-4 4 6-5 5 6 "^d, 6 , 5 


As before, when the radiosities are evaluated, Q*, Q* and Q* can 

4 .5 6 

be readily calculated. 

V?hen the sun's rays are normal to the solar panel, the relevant 
radiosities can be evaluated from (3.2.7b) provided that the subscripts 
1, 2 and 3 are replaced by 4, 5 and 6', respectively. 

By following the same procedure, we may evaluate Q* and Q*j^ • 

The results are: 

(a) For oblique solar irradiation 


QJ = A s cos 45« sin 30» 

^d, 9^d, 11 9-13 11-9 


O* - A Ct«> — : 

^11 11 11 1 - p* p* F F 

^d, 9rd, 11 9-n 11-9 


S cos 45® sin 30® 


(3.2,10a) 

(3.2.10b) 


(b) For normal solar irradiation 

It is only necessary to replace cos 45® in (3.2.10a,b) 
by unity. 

The top face A^ of the spacecraft bus is completely shielded from 
the sun's rays by a multilayer superinsulation. The Appendix of' this 
report gives a detailed analysis of the thermal resistance of the 
superinsulation and presents an expression from which the leakage 
flux may be determined. The rear and bottom face Ag and A^^ of the 
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bus do not receive radiation in the solar spectrum » either directly 
or indirectly. Hence, Q* = = 0. 

3.2.2 Radiation Heat Transfer in the Infrared 


Once again for generality, v;e first consider an enclo- 
sure of N surfaces, each having a uniform temperature and being uni- 
formly irradiated. It is assumed that all emitted radiations are 
diffuse and that the surface reflectances can be adequately described 
by the p -p model. Our objective here is to deduce expressions 

d s 

for the overall radiant exchange factor similar to Hottel's diffuse 
exchange factor All surfaces are taken as gray within the spec- 

tral range. Clearly, the results will be directly applicable to the 
determination of similar exchange factors for internal radiation with- 
in the spacecraft bus. 

For any surface A. of the enclosure. 


B. 


e. OT. + p . H. 

1 1 d, 1 • I 


and 


H. = 



j = 1 


(3.2.11) 


(3.2.12) 
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Thus , 


N 

B. = e. 0T^ + / E. .B. ('3. 2. 13) 

I 11 ‘^d, 1 X ! » J 

j = 1 

The system of simultaneous, linear, algebraic equations (3.2.13) can 
— be- rewritten as 


N 



a . B. 


1 3 J 


j = 1 


= e aT. 

1 I 



(3.2.14) 


where , 


a. . = 6. . - . E. . 

1 ] 1 J d, 1 1 - J 


(3.2.15) 


To evaluate the diffuse-specular overall exchange factor-^^ , we 

i *“k 

4 

set OT. = 1 and the temperatures of aVL vemainlng surfaces to zero. 
The heat flow rate is 


i %, = \“k i“k = w 


. . B. 

k -J 1 J 


(3.2.16) 


where the presubscript i is a reminder of the original source of the 

flux. We recognize that . ( j = 1, 2, • • • , N) is the solution 

set of (3.2.14) v;ith the latter's right-hand side set to zero except 

for the- ith equation, for which it is £. . It is most convenient 

to evaluate . B by matrix inversion. 

1 J 
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V7e now proceed to apply the foregoing general results to evalu- 
ate ^ for the spacecraft’s exterior surfaces. V7e first consider 

the enclosure formed by A . and A , , the latter being black 

■^12 3 e,L ° 

and at absolute zero. 


Ai , , A^ are p^-p^ surfaces 

A is black 

e , L 

2(3) 1(3) 


3(1) 

3(2} 


Figure 3.3 Mirror Images of Interacting Surfaces 

in an Enclosure Consisting of A , A , A and A 
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In Fig. 3.3, the surfaces under consideration are identified v/ith 
letters or numerals irmediately adjacent to them. Thus, refers 
to, the upper surface of the solar panel and A^ the exterior surface 
of the left face of the spacecraft bus . The image surfaces are des- 
ignated with two numerals ; the one in parentheses refers to the mirror 
in which a surface in question is seen. For instance, 2(3) is the 
image of A^ as seen in the mirror A^ , With the aid of Fig. 3.3, 
the various exchange factors can be readily written. They are: 



In (3.2.20a,b), the subscript e refers to A . Using the foregoing 
information, the diffuse-specular overall radiant exchange factors 

and can be evaluated immedi- 

1 — .z 1 — 1 — c z — ^ z —e o — e 

ately from (3.2.19). Exchange factors associated vrith A^ , A^ , 
and A can likewise be determined. Surface A sees only the empty 

6 I i\. o 

space, thus the relevant ^ factor becomes its surface emittance. 
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Next, we consider the enclosure formed by (upper sur- 

face of antenna) and the associated black surrounding A . Figure 

e, F 

3.4 shows the two relevant mirror images 9(11)- and 11(9) associated 
with the system. 



Figure 3.4 Mirror Images of Interacting Surfaces in an Enclosure 
Consisting of and A^ ^ 


It is easy to see from Fig. 3.4 that the exchange factors are (with 
subscript e referring to A ) 

e , r 


^9-9 S-11 ^9-11 


^11-11“^’ 1 -e ~ ^1 1 -e 


^9 -re ~ ^9-e , 1 1^9( 1 

^s, 9^11( 9) -e 


1) -e 


(J 


and the corresponding ^ j ’ s are 


^9-9 ■ ®9-ll ^d. 9^9-11’ ®9-e "'^d,9^9-e 


11-9 


Pd, 11^11-9’ ®ll-e Pd, 11^1] 


(^ 


The remaining system of the spacecraft's exterior surfaces for 
which the ^ factors are needed consists of the dark side of the solar 


.2.21a) 


.2.21b) 


41 



panels (A^,, A^,, A^, and A^,), the outermost plastic film A^^^, of the 
superinsulation covering A^^j’the shaded side of the antenna 
and the associated emptj^ space . Since none of the surfaces involved 
has a specular component of reflectances in the infrared region,, the 
evaluation of a_ and, subsequently, ^ requires information on 
shape factors, but not on exchange factors. The formulation presented 
in this section and the expressions (3.2.11) to (3.2.19) remain valid. 
All one needs to do is to replace p, . by p, and E. . by F. . . It 
goes without saying that the factors for interior surfaces of the 
bus can be likewise evaluated. 
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3.3 NODAL HEAT BALANCE EQUATIONS 


To formulate the nodal heat balance equations, we need to evalu- 
ate, in addition to the overall radiant exchange factors the con- 
ductances of the various heat flow paths . VJe shall assume that there 
is no conductive link between the solar panels and the spacecraft bus, 
nor between the bus and the antenna. Further, the eleven sections 
(22 surfaces) which make up the spacecraft are taken to be identical 
squarest as illustrated in Fig. 3.2. Each section may have a differ- 
ent thickness, however. Under the foregoing conditions, the conduc- 
tance between two adjacent areas, A. and , having a common edge is 


(kd). (kd). 

' - 2 ^ ^ - 

\-i (kd). t (kd). 

j j 


(3.3.1) 


Clearly , K = K . If A and A are fabricated from the same ma- 
1 - j J - 1 J ] 

terial and are of identical thickness , like A and A (or A and A ) 

. 12 4 5 

of the solar panel, then (3.3.1) reduces to 

K. J = (kd). = (kd). (3.3.1a) 

1 -] 1 } 


It is pertinent to note that the crude subdivision used in the 
analysis should in no way jeopardize our present objective since 
the model is to be subdivided in precisely the same manner. In this 
regard, vje may further add that, solely from the vie;/point of ascer- 
taining the usefulness of the proposed theory, it makes no difference 

tA significant simplification of the arithmetic results from this 
choice. It is not a restriction necessarily imposed on the analy- 
sis . 
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whether a continuum model or a discrete model is used in the analy- 
sis . 

■The top and bottom faces of the spacecraft bus are covered with 
superinsulation. The exposed layer of the superinsulation for the 
bottom face interacts radiatively with and, thenceforth, with 
Aj,, A^,, A^, and A^,. An exact analysis of the thermal transients as- 
sociated with its many layers is too involved and is not warranted 
for the problem under consideration. If we ignore its heat capacity, 
then a reasonably simple analysis can be made. It is demonstrated 
in APPENDIX .A that the outward leakage flux (from the base plate A^.^ 
toward the outer surface of the insulation) is 


Q 

A 


CJ(T 


10 


- ) 
10'^ 




(3.3.2) 


where = (n - 1)R + R' . All quantities have been defined in 

APPENDIX A. 

We now proceed to present the heat balance equations . The di- 
mensional form will first be given: Set (3.3.3) 


dT, 


(Cd) 


1 dt 


it - '^-2 - ’’ 2 > - 

- )ctT^ - 0 (T^ - t"^ ) 






^5^ 


^ 0(T^ - t"* ) 


- 0 (T^ 

I'-ii' ^ 1 


4 ^ 1-2 

T ) i—f- fT - T ) 

11 ^ A 1 2 ^ 
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^2 = 


(Cd) 


2 dt 


= q|c _ - T^) - ^_^a(T; 


/) 

3 ‘ 




^ ^ + <= 5 ? ) aT "^ - 

- ^ a ( T ^ - t "^) - ^ a ( T ^ 

•^ 2 '- 5 ' '• 2 ^5 ^ ■ ^ 2 '-l O' ° ^ ^2 


- ^ 4 ^ 

T^ ) 

10 '^ 


- ^ a(T^ 
2'-lV ^^2 




K 


2 -1 


( T . 


T ^) 


\-' 


(Cd), 


3 dt 


- ^ Cf(T^ - ^ CffT^ - ) 

•^-8 ^ S 8 *^-9 ^ 3 9 '^ 


- ,„ o ( t ! - t !„) - 


3-10 ' 3 10 


(T - T ) 
A ^ 3 - 7 ^ 


^-8 ^ 9 

-L-L (T - T ) - -— -- (T - T ) 
A 3 8 ^ A 3 9 ^ 


- ^- (T - T ) 
A 3 10 -' 


^4 ‘ 


(Cd) 


5 

4 dt 




- (^_e + i'.-e'K - ^4'-l> 

“ “ t!) 


^ a(T^ 

4 '- 10 ' 


) 


- ^ 

4 '- 1 - 1 ' 


a ( T , 


K 


- 


4-5 


( T . 


T5 ) 
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- T*) - 

- „ oCrt - ,„, a(tf - <„,) 


- ^ a(T^ - ) - 


(T. - T, ) 

5 4 


'J|' + % - ^6-4®<^6 - ’'s' 

- - ’s' - - <> 


- ^ oCt"^ - ) (T “ T ) 

'6-10^''^6 10 •' A ^6 7-' 


!ii (T _ T ) - ^ (T - T ) - (T 

A '■a ^8'' A ^6 9^ A ^6 


A— ( s cos e - aT 

^ \e s 7, 

V p 


) - 


3<'<< - <' 


^ Cf(T - T ) - ■^ oCT(T' - T ) 

7-6 ^ 7 6 -^ ^7-8 ^^7 B 

^ o-(T^ - T^) - ^=5^ - T^ ) 

^7-9 ^ 7 9 7-10 ^7 10^ 

K K , K „ 

A '■7 3^ A ^7 6^ A ^7 


(Ty “ Tg ) 
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^8 = 


(Cd), 


'6 dl- 




<5 




- ^ a (/ - )' 


8-10 ' 8 


K 


8 -3 


('!« - Is ) 


* 8-6 '* 8-7 

(T - T ) (T • “ T ) 

A ^8. ^ 6 -' A '■8 7 ^ 


K 


' 8-10 


^^8 - 




(Cd) 


dT, 
9 dt" 


9 


qi ‘‘ + % 


^ a(x 

9-11^9 


'll 


t;. ) - ^ ot 


,4 

9 - e '^'^9 


‘5^7-s'’«9 - I5’ - 




^ 9 - 7'’«9 - <> - - Ia > 


- ■^ a ( T ^ 
9-10 ^9 


) _ _ 1 :±. (T - T ) 
10 A 9 3 


K. 


^ (T 
A ^9 


. ^ , ^ 9-10 

^6^ ~ — “^7^ --T- 


(T„ 




dT 


10 


2 k 


ct(t‘^ - T ^ ^ a(T - T^O 

10 30'^ -uo-s 10 3 



a(T\ 

- tS - 


a(T^„ 

1 0 —6 

10 

6 

10-7 

10 

^ a(T^ 
■^0-8 10 

- - 

^0-9 10 


<) 

<) 


10-3 


10 


K, 


(T 

A ^10 


T3) 

10-6 

A 

(T, 

"s> 

K 

10-9 

(T, 

A 


10 


10 
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0 


1 




2 k 


4 4 4 4 

a(T’ - T ) - aCT - T’ ) 

10'^ -lO'-i' ^ '•So' *1 ■' 


"‘^0'-2' ■'^2^ 
■^•O'-e' ‘^'^10' 


^ -a(T^ - ) 

^ a(T^ - T^ ) 


1 1 


dT 


(Cd) 


1 1 


11 dt 


= n* - ^ a(T^ - T^) - ffT^ 

Si 11-9 ' 11 9 ^ 11-e 11 




ll'-l'^(^ll 


^ aCT"^ - ’t ') ~ ^ 

-=^1 11 _ 4 ' '•' '^ 1 1 11 ' - 5 ' 


) - ^ 


1 1 ' - 2 ' 


cr(T^, - 

) 

1 1 

2 

1 

b 

N 

) 

1 1 

S 

4 
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Nodal Heat Balance Equations in Dimensionless Form 

To reduce the system of heat balance equations to appropriate 
dimensionless forms , we introduce the following nondimensional vari- 
ables : 


T. 


1 

(3.3.4a) 

o 


it 

C. A 

3 

(3.3.4b) 

3 

ctT a 


o 

(kd). 

(3.3.4c) 


where i = 1, 2, 3, etc. j and is a sui ■cable reference temperature 
for the problem. Its choice is, to some extent, arbitrary. For 
instance, may be taken as the initial uniform temperature of the 
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spacecraft or some representative average of its steady state tempera- 
tures prior to the initiation of thermal transients . Since the physi- 
cal time t must necessarily be the same for all i's, the dimension- 
less time T. for the various component parts’ ^ must be so related 
that: 




C ‘ C 


By using (3.3.4ajb5c) , we obtain the following non-dimensional form 
of the nodal heat balance equations: , set (3.3.6) 


1 ■ dt. 


- “^ 1-2 <>2 - '^.- 4 . <- < - ^.- 40 . < 0 .- <> ll ' 


(kl) T ■ ^2 [^-1 "^-3 ^-e ^ 

2 o 

- "^2-sK - ‘^2'-4' *5 - ^2.-10. «'l ,■ 


1 ^ ^ 

2 - 3^3 2 '- 


2 '-4' ^4 


'^2.-1, 4>k] - <h - h ) 
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di3 

dT_ 


(q| + qj )A 

(kd),T 

3 o 


- h ['“f-l ^ ^3-2 - + ^-6 

^ ^-2 - ‘I-. - t 


'3-9^9 


- ^ 

3-10^10 


” (kd) ^^-7 ^-8 ’*’ ^-9 '*' ^-10 ^'^3 


'*3-7 , , ’*3-8 ^ **S-9 ^ ■ ■*3-10 ^ 

■■ TmT ^7 + TkdT *’« ■" TkTT ♦o ■■ -(kdr *10 

o o 3 


^4 • 


dT. 


q|A 

Ckd^ - 3^ 


'■ -^4-5 ^ ^ M-e M'-l' ^ '^'-2' 


4 -6 ^6 4'- 1' ^1 ^i' -2' ^2 M'-l 0' ^1 O' 


- 


4'-l 1 




K 
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d^3 ^ q*A 
dT " (kd) T ^5 
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K, 


(Sr ‘♦5 “ *0^ 
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dT 


d*, (tf t <1, M ^ t ‘51.3 ■; 

“0?d)Jo 6 L ^ 

_ <5Z - ■??_,'(>* - - ‘^6-9*» ■ '^‘->“ >“J 






\-> 


10 


d<|>. 


\ ' dt. 


a* SA cos 6 
i P 
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51 


iJ 



+ H- )(j)^ - d)^ - ^ 

•^9-7 ^ 9-8 -^-10 ^^9 ^9-11^11 ^9-3^3 

_ ^ ^ J,'^ _ 4? 6 ^ - ^ d,- 1 
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(kd)^Q ‘^6 ■*■ (kd)^Q ^8 ^ (kd)^^ *^9 
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A. 


qf^A 


n 


dx. 


1 1 


(kd)^ T 

II O 


2^ 


~ ^11 [^^ 1-9 ^ll'-l’ ■*■ ‘^11’- 

_ ^ ^ d>^ - ^ d)^ - ^ 6^ • 

"^ii'-io' '^ao'] 


In numerical computations , if we arbitranily select X^ as the "stan- 
dard" dimensionless time, then 


d(J)^ 

/b '=1'! 

dlf)^ 

dx. 

\k, cj 

dr 

1 

' 1 3' 

3 

d<p2 
2 ^ / 



( 

'V 




3.4 MODELING REQUIREMENTS 


An inspection of the dimensionless heat balance equations reveals 
that, for perfect modeling, the following scaling requirements must 
be met. 

1. Model and prototype are geometrically similar in order that 
all configuration factors among corresponding surfaces are 
identical. (Thickness distortion vrould have very minor or 
little Influence on radiation exchange and is thus permitted. ) 

2. Radiation properties of all corresponding surfaces are iden- 
tical. 

The satisfaction of (1) and (2) ensures that the model and the proto- 
type have the same for all corresponding surfaces . 

3. The simulated solar radiation used in model testing has 

a spectral distribution , intensity and direction the same 
as those of the local sun's rays incident on the spacecraft. 
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4 . 


The (kd)'s and (Cd)'s of all corresponding surfaces must 
satisfy the following relationships; 


(kd) 

i) 


m 


a fixed quantity 


(3.3.7a) 


(Cd) 

TcdT 

p 



a fixed quantity 


(3.3.7b) 


The first is associated with the geometric scale ratio and 
the second the time scale ratio. 


5. The internal power dissipation, when expressed in terms of 
the equivalent surface flux, is the same for the prototype 
and the model , i . e . , 




for all corresponding surfaces of the spacecraft bus. 


(3.3.7c) 


Finally, we note that, for complete similitude, the initial 


temperature field of the model and of the prototype must either be 
uniform or have identical distributions. The foregoing modeling 
requirements have previously been stated in the literature, e.g., 

I 

[5], and experimental verifications are available [6a,b,c!]. 
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Chapter 4 

SOLAR FLUXES, DIFFUSE-SPECULAR OVERALL EXCHANGE 
FACTORS AND CONDUCTANCES OF HEAT FLOW PATHS 
-ASSOCIATED WITH THE HYPOTHETICAL SPACECRAFT 

As a prerequisite for the numerical solution of the nodal heat 
balance equations, the solar fluxes Q*/A , the overall radiant exchange 

i I 

— factors ^ and the conductance of the heat flov? paths K associated 

i-j . ■ i-J 

vrith the various surfaces of the hypothetical spacecraft must be 
evaluated. All relevant conflg^u?ation factors can be readily deter- 
mined since all surfaces or their subdivisions are chosen to be iden- 
tical .squares as illustrated in -Fig. 3,1. The results are listed 
in TABLES 4.3 and 4.8 for surfaces facing the sun and in TABLE 4.13 
for surfaces away from the sun. The data were taken either directly 
from an NACA Technical Note by Hamilton and Morgan [ 7 ] or indirectly 
evaluated from the information given therein. The solar exchange 
factors were calculated from (3.2.6a or b) and (3.2.8) and the solar 
radiosities from (3.2.7a or b) and (3.2.9). The irradiation H* is 
then given by (3.2,2) and the required Q* follovjs from (3.2.4) or 
(3.2.5). For surfaces A and A , the solar fluxes can be directly 

9 11 - 

determined from (3.2.10a or b). In these calculations, the solar 

2 

. constant S was arbitrarily assumed to be 442 Btu/hr-ft and the rele- 
vant reflectances and absorptances were taken from TABLE 3.2. The 
solar fluxes so calculated are listed in FABLES 4.1 and 4.2, respec- 
tively, for the normal and oblique incidence of sun’s rays. 
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The procedure used in the determination of the diffuse-specular 

overall exchange factors for the spacecraft's &xteTior> surfaces has 

been explained in detail in Chapter 3. Following the evaluation 

of the exchange factor E ^ , the elements of the coefficient matrix 

a were calculated according to (3.2.15) and those of its inverse 
i-j ■ 

b . were determined by using an available subroutine in our conpu- 
• ter laboratory. For the convenience of the reader, these intermedi- 
ate data are listed in TABLES 4.4 and 4.9 for E. ^ , in TABLES 4.5, 
4.10, and 4.14 for a, and in TABLES 4.6, 4.11, and 4.15 for b. . . 

i-j i-j 

The desired diffuse- specular overall exchange factors ^ were cal- 
culated from (3.2.19) and the results are listed in TABLES 4.7 and 
4.12 for surfaces facing the sun and in TABLE 4.16 for surfaces away 
from the sun. It is pertinent to note that, for an enclosure of 
N surfaces, the overall exchange factors satis:^ the following simple 
relation 


N 

V 4^. . = e. (4.1) 

Z_/ > “i » 

j = 1 

which is the consequence of the requirement of energy conservation. 

Equation (4.1) has been used to check the numerical accuracy of the 

data list in the three tables last mentioned. 

\ / 

For the interior of the spacecraft bus , the radiation v?as as- 
sumed to be devoid of specular component and, hence, only the dif- 
fuse overall exchange factor is involved. TABLES 4.17, 4.18, 4.19, . 
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and 4.20 list, respectively, the data for the configuration factors 
F . , the elements of the coefficient matrix a. . , the elements of 
the inverted matrix b. . and the diffuse overall exchange factors 

i-j 

Finally, the conductances are listed in TABLE 4.21. They 

vfere evaluated from (3.3.1) using information' given in TABLE 3.2. 

It may be noted that there is no conductive link between the solar 
.panels and the spacecraft bus, nor between the bus and the antenna. 
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note: exterior of surface 10 IS DESIGNATED AS SURFACE 12 





TABLE 4.1 SOLAR FLUX AT 

NORMAL INCIDENCE, 

Q4/A, BTU/HR, SQ.FT 



1 

2 

3 

4 5 

6 7 

8 9 10 

11 

12 

353-5999 

353.5999 

0.0 

353.5999 353,5999 

o 

t 

o 

o 

» 

o 

0.0 0.2376 0.0 

154.7000 

0.0 





TABLE 4.2 

SOLAR FLUX 

AT 45 OEG 

INCIDENCE, 

Q^=/A, BTU/HR, SQ.FT 




1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

250.0340 

0.0 

0.0 

257.9092 

298.0498 

75 .0102 

0.0 

0.0 

0.1680 

0. 0 

109.3899 

0.0 


^ ^ SURFAOES ^ Iv 


TABLE 4,3 CONFIGURATION FACTORS FA(j,J) =:=*» 


I / 


2 


3 


E 


1 

0.0 

0, 0 

0.032810 

0.967190 

2 

0-0 

0.0 

0.200040 

0.799960 

3 

0.032810 

0.200040 

0.0 

0.767150 

1(3) 

s[« ^ 4f3{< ^ ^ 


0.032810 

0. 032810 

2( 3) 


<Y^ ^!* 

0, 200040 

0.200040 

3(1) 

0.032810 

;fr 


0-032810 

3(2 ) 


0.200040 

^ ^ ^ j ■# 5jt 

0.200040 


FOLDOUT FRAME | 






FOLDOUT FRAME 








1 0.0 0.0 0.032810 0.967190 

2 0^0 OjJ 0. 200040 0.799960 

3 0.032810 0,200040 0.0 0.767190 





2 0.0 1.000000 -0.040008 -0.159992 

3 -o'. 005578 -0,034007 1. 000000 -0.130422 


E 0_5^0 0^0 1, 000000 



•4- uu ^ fO ^ in vO 
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6 -0.005578 -0.034007 l.OOOOQO -0.130422 

E 0 _ i 0 0^0 0^0 1. 000000 


TABLE 4.6 INVERSE OF COEFFICIENT MATRIX B0U»J ) 


I / J 1 2 3 E 


1 1.000 036 0.000223 0.00657 1 0. 1443 38 

2 0.000223 1.001362 0.040064 ' 0.165478 

3 0.005585 0.034054 I.Q0I399 0.137133 

E 0.0 0.0 0,0 1.000000 


*** TABLE 4.6 CONTINUED 


/ 

j 

4 

5 

6 

E 



1.000036 

0.000223 

0,000223 

1.001362 

0.006571 

0,040064 

0-194338 

0.165478 



0.005585 

0.0 

0.034054 

0.0 

1.001399 

0.0 

0.137133 

1.000000 






TABLE 4.7 DIFFUSE- 

■SPECULAR 

OVER-ALL EXCHANGE FACTORS FFU,J1 

/ 

J 

1 

2 

3 

E 



0.000117 

0.000715 

0.021816 

0,777351 

! 


0.000715 

0,021816 

0.004360 

0.133012 

0.133012 

0-005670 

0.661913 

0.659534 





'll* 

TABLE 4.7 CONTINUED 


/ 

J 

4 

5 

6 

E 

U 


0.000117 

0.000715 

0.021816 

0.777351 

i 


0,000715 

0.021816 

0.004360 

0,133012 

0. 133012 
0,005670 

0,661913 

0.669534 



to 



— 




foldout frame 
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TA3LE 4.8 CONFIGURATION FACTORS FA(ItJI 
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1 / 0 

9 

11 

E 


9 

0.0 

0-021500 

0.978500 


11 

0.021500 

0.0 

0. 978500 


9U1! 

0.489250 

0.489250 

0.0 


11(9) 


0. 021500 

0.021500 






TABLE 

4.9 EXCHANGE FACTORS E(I,J) 







I / 

J 

9 11 

E 






9 


0.0 0.021500 

0.978500 






11 


0.021500 0,0 

0,998925 





*** TABL£ 4.10 COEFFICIENT MATRIX AAil,J) 



I 1 

0 

9 

11 

E 



9 


1.000000 

0.0 

0.0 



11 


-0.003870 

1. 000000 

-0,179806 



E 


0.0 

0.0 

1. 000000 



*** TABLE 4.11 INVERSE OF COEFFICIgNT MATRIX 66(1 , ^}_ 


I / 

J 

9 

11 

E 



9 


1.000000 

0.0 

0.0 



11 


0.003870 

1.000000 

0.179806 



E 


0,0 

0.0 

1. 000000 




TABLE 4,12 DIFFUSE-SPECULAR DVER-ALL EXCHANGE FACTORS PFil.J} 


I / J 9 H E 


9 0.000000 0.000881 0.049118 

11 0.000881 0.0 0.819118 


FO] nniiT / HOLDOUT mWE 

iiS“TABL^4.13 CONFIGURATjaiM FAC TORS ?ATf7jT^ " ' " 
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< - 1 / 

J 

1 ■ 

2' 

10’ 

5' 

4' 

11' 

E 

1' 

2' 


0,0 

0.0 

0. 0 
0.0 

0-0 

0-0 

0-0 

0.0 

0.0 

0-0 

0.002700 
0, 019500 

0.997300 

0.980500 

10' 

5’ 


0.0 

0.0 

0-0 

0.0 

0-0 
0. 0 

0-0 

0.0 

0.0 

0.0 

0-037000 

0.002700 

0.913000 

0.997300 

4* 

11 • 


0.0 

0.002700 

0.0 

0.019500 

0-0 

0.087000 

0.0 

0.002700 

0.0 

0.019500 

0. 019500 
Q.O 

0. 980500 
0.868600 





Jr 

TABLE 4.14 

COEFFICIENT 

MATRIX AA! I 

,J) =>'*=5= 



I / 

J 

1' 

2’ 

10' 

5' 

4' 

11 » 

E 

1’ 


1.000000 

0-0 

0. 0 

0. 0 

0.0 

-0.000459 

-0.169541 

2* 

10 ' 


0.0 

0.0 

1.000000 

0.0 

0-0 

1. 000000 

0-0 

0-0 

0.0 

0.0 

-0.003315 

-0.043500 

-0. 166685 
-0.456500 

5’ 

4' 


0-0 
0. 0 

0.0 

0.0 

0.0 

0.0 

1- 000000 
0.0 

0.0 

1.000000 

-0.000459 

-0,003315 

-0.169541 

-0.166685 

11' 


-0.002565 

-0.018525 

-0. 082650 

-0,002565 

-0,018525 

1.000000 

-0,825170 




TABLE 4.15 

INVERSE OF COEFFICIENT MATRIX BBII.J) 




I / J 

1* 

2* 

10' 

5 ' 

4' 

11' 

E 

1' 

2' 

1.000001 

0.000009 

0.000009 
1 .000061 

0. 00Q038 
0.000275 

0, 000001 
0.000009 

0. 000009 
0.000062 

0,000461 

0.003327 

0.169942 

0.169580 

10* 

5' 

0.000112 

0.000001 

0.000809 

0.000009 

1. 003608 
0.000038 

0.000112 
1, 000001 

0,000809 
0. 000009 

0.043662 

0.000461 

0.494434 

0.169942 

4« 
11 ' 

0.000009 

0.002575 

0.000062 
0.0 18594 

0.000275 
0. 082959 

0.000009 

0-002575 

1.000061 

■0,018594 

0.003327 

1-003734 

0.169580 

0.873194 




-A. 

TABLE 4. 16 

DIFFUSE OVER- 

ALL EXCHANGE 

FACTORS FB(1 



1 / J 

1' 

2' 

10* 

5' 

4» 

11» 

E 

1' 

0.000005 

0.000035 

0.000093 

0-000005 

0. DQ0035 

0.000112 

0.829715 

2' 

10' 

0.000035 

0,000093 

0.000250 

0.000671 

0.000671 
0. 001804 

0.000035 
0- 000093 

0.000250 

0.000671 

0.000812 

0.002183 

0-827947 

0.494483 

5' 

4' 

0,000005 

0.000035 

0.000035 

0.000250 

0.000093 

0.000671 

0.000005 

0.000035 

0-000035 • 
0.000250 

0.000112 

0.000812 

0.829715 

0.327947 

11' 

0.000112 

0.000812 

0,002183 

0.000112 

0.000812 

0-000010 

0,045958 


FOLDOUT FRAMF \ 


FOLQ.OUT EBAiVlE. 



interior surfaces of bus 
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TABLE 

4.17 CONFIGURATION FACTORS FA(I,J) 



I 

/ 

j 

3 

6 

7 

10 

8 

9 

3 



0.0 

0.199820 

0.200040 

0.200040 

0.200040 

0.200040 

6 

7 



0.199820 

0.2Q0040 

0,0 

0. 200040 

0.200040 

0.0 

0.200040 

0.199820 

0.200040 

0.200040 

0. 200040 
0,200040 

10 

8 



0.200040 

0.200040 

0. 200040 
0.200040 

0. 199820 
0.200040 

0.0 

0.200040 

0. 200040 
0,0 

0.200040 

0.199820 

9 



0.200040 

0.200040 

0.200040 

0.2 00040 

0.199820 

0.0 





TABLE 

4.18 COEFFICIENT MATRIX 

AAIItJ) 4* 




I 

/ 

j 

3 

6 

7 

10 

s 

9 

3 

6 



1.000000 

-0.159856 

-0.159856 

1-000000 

-0. 160032 
-0.160032 

-0,160032 

-0.160032 

-0.160032 
-0. 160032 

-0,160032 

-0.160032 

7 

10 



-0.030006 

-0.030006 

-0. 030006 
-0.030006 

■1.000000 

-0.029973 

-0.029973 

1.000000 

-0.030006 

-0.030006 

-0,030006 

-0.030006 

a 

9 



-0.030006 

-0.030006 

-0.030006 
-0. 03D006 

-0,030006 
-0. 030006 

-0.030006 

-0,030006 

1.000000 

-0,029973 

-0.029973 

1-000000 




TABLE 

4.19 INVERSE 

: OF COEFFICIENT MATRIX 

BB ( I J J ) 



I 


J 3 

6 

7 

10 

8 

9 


3 


1.057707 

0. 195533 

0, 220391 

‘ 0.220391 

0.220391 

0.220391 


6 

7 


0.195533 

0.041323 

1.057709 

0.041323 

0.220391 

1.017411 

0.22Q391 

0.046514 

0.220391 

0.046545 

0-220391 

0-046545 


10 

8 


0 .041323 
0-041323 

0.041323 

0.041323 

0-046514 

0.046545 

1.017411 

0-046545 

0.046545 

1-017412 

0.046545 

0.046514 


9 


0.041323 

0-041323 

0. 046545 

0.046545 

0- 046514 

1.017412 





444 TABLE 

4.20 DIFFUSE OVER-ALL 

EXCHANGE FACTORS FF( I.,J 3 

'f' T' 



I 

/ 

J 3 

6 

7 

10 

8 

9 


3 

6 


0.002885 

0.009777 

0.009777 

0.002885 

0.046833 

0.046833 

0.046833 

0.0463,33 

0. 046833 
. 0.046833 

0,046833 

0.046833 


7 

FOLDOUT 10 


0 .046833 
0,046833 

0-046833 

0.046833 

0. 083875 
0.224043 

0.224043 

0.G83875 

0.224193 

0.224193 

0.224193 
0. 224193 


' 8 


0. 046333 

0. 046833 

0.224193 

0.224193 

0-083875 

0-224043 
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9 0.046833 0.046833 0.224193 0,224193 0.224043 0.033875 


44 * TABLE 4.21 DIMENSIONAL CQNDUCTIOM LINKS AKII , J ) , BTU / HR,oeG * 4 * 


I / J 

1 

2 

3 

4 

5 

8 

1 

0.0 

0.058000 

0.0 

0,0 

0.0 

0.0 

2 

0.058000 

0.0 

0‘.0 

0.0 

0 . 0 

0,0 

3 

0.0 

0 . 0 

0 . 0 

0.0 

0.0 

0-0 

4 

0.0 

0,0 

0.0 

0.0 

0 - 058000 

0.0 

5 

0 . 0 

0.0 

0.0 

0.058000 

0-0 

0 . G 

6 

0.0 

0.0 

0.0 

0-0 

0.0 

• 0.0 

7 

0.0 

0.0 

0.071100 

0.0 

0 . 0 

0. 074391 

8 

0.0 

0 . 0 

0,471490 

o ^ o " 

0.0 

0.667263 

9 

0.0 

0.0 

0.560795 

0.0 

0.0 

0.861395 

10 

0 . 0 

0.0 

0.033714 

0.0 

0.0 

0. 034437 

11 

0.0 

0-0 

0.0 

0.0 

0.0 

0.0 

12 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


444 TABLE 4.21 CONTINUED 44 * 


I / J 

7 

8 

9 

10 

11 

12 


1 

0.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 


2 

0.0 

0.0 

0.0 

0 , 0 

0 . 0 

0.0 


3 

0.071100 

0-471490 

0-560795 

0.033714 

0.0 

0.0 


4 

0.0 

0 , 0 

0.0 

0.0 

0.0 

0.0 


5 

0.0 

0.0 

0.0 

0.0 

0.0 

0,0 


6 

0.074391 

0. 667263 

0.861395 

0.034437 

0 ,0 

0.0 


7 

0.0 

0.070706 

0.072435 

0 . 0 

0,0 

0.0 


8 

0.070706 

0.0 

0.0 

0.033626 

0,0 

0.0 


9 

0.072435 

0.0 

0.0 

0.034012 

0.0 

0-0 


10 

0.0 

0.033626 

0.034012 

0.0 

0 , 0 

0 . 0 


11 

0.0 

0-0 

0-0 

0.0 

0.0 

0.0 


12 

0.0 

0.0 

0.0 

0.0 

0,0 

0.0 
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CHAPTER 5 


- NUMERICAL SOLUTION OF NODAL MEAT BALANCE EQUATIONS 
5.1 A RECAPITULATION 


The dimensionless nodal heat balance equations for the hypo- 
thetical spacecraft presented in Chapter 3 and designated by the 
set (3.3.6) are simultaneous, non-linear, first order ordinary- 
differential equations. They can be compactly described by (5.1.1) 
which is written for an arbitrary node i. 


^2 




i_ 

(kd) 


2 'S-i - 2 

j 5^ i j i 


v?here 


(5.1.1) 



M. 



k C. 

0 1 

~ k. C ’ 

1 o 


subscript o refers to a convenient reference 
surface for V7hich the dimensionless time T 
is related to the physical time t according 
to T = (k /C A)t. One of the surfaces in 

O O 

the system may suitably serve for this pur- 
pose. 
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(q* + )A 

Qj = ^ 5 except vfhich is 


1 KJ4 

p 


a*'' S A cos 9 


ZR 


Also, = Qi 2 = 0 


C. . 

1 -j 


conductance of multi-layer insulation, (2 r) 
The three nonvanishing conductances are ^ 


' 10-12 


and C 


1 2- 1 0 


42 '. 

i-j 


K. 

1 -J 


overall radiant exchange factor between A^ and A^. . 

The indicated summand refers to all A. 's vihich have 

) 

direct and/or indirect radiant exchange with A. . 

V/hen approriate, A includes the empty space. A 
similar interpretation should be given "to^^ ^ i* 

j’ 

= conductance betvreen A and A . Here , the indicated 

1 1 

summand refers only to those surfaces v?hich have physi 
cal contact with A. . 


All other quantities have previously been defined in Chapter 3. 

The introduction of the dimensionless parameter is to obtain 
a uniform time scale for all component -surfaces in the system. An 
alternative is to replace the left hand side of (5.1.1) by 


C. A d* 

1 b 
k. dt 

and the integi'ation is to be carried out in the physical time domain. 
The latter is indeed adopted in our computer program for v/hich separate 
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documents v;ere prepared and submitted to JPL. 

5.2 NODAL HEAT BALANCE EQUATIONS IN FINITE DIFFERENCE FORM AND THE 
METHOD OF SOLUTION 

Equation (5.1.1) may be integrated from T to T + At, At being 
"a“small increment , to give 


M. [<!>. (T + At) - <J). (T): = Q.Ax - 

3 j' 

+ fit ^ fit 

3 

3' • * j ?? i ‘ 

/ t+At ■ 

imply J . They will now be replaced by the 


the approximate, two point integration formula, namely. 


T + At 


4> dt = At[( 1 - a)^(T) + dcp(r + Ax)] 


T t At 


J° dr = At[(1 - a)(|)^(T) t (t + At)] 


(5.2.2h) 


in vdiich a is the so-called Integration parameter and is a pure numeral 
between 0 and 1. This scheme of approximation has been used by Strong 
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and Emslie [ 8 ] in their treatment of relatively simple systems in- 
volving combined radiation and conduction heat transport. Upon 
substituting (5.2.2a,b) into (5.2.1), followed by some rearrange- 
ment, one obtains 

f^ . d)^(T+ At) + f, . (t + At) + f _ . =0 

4, 1 U 1,1.1 0, 1 . 


where 


f = aa (c. t V .+ V .,) 

4,1. Z-/ • >-) Z_/ I'-j’/ 

M. 


a 


2 


1, i . (kd). .i— / 1 -3 

j 5^ i 


At 


M. 


^0, . . = -% - * Tkir E *=1-) 

' j H i 


+ a-coe.[c,_, + '2 

j j' 


<|)^(T) 


-fer S 


D ^ i 


- (1- [s-i 2 


+ (j)^, (t)] - a3. [c._.^j^(T + At) 

j’ 


(5.2.3) 

(5.2.3a) 

(5.2.3b) 
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(5.2,3c) 


+ 2 + At) + 2 


Here again, i. , (t) = 4 (t) and 4, (t + A'l) = 4 (t + At). 

The equation set (5.2.3) v.^as solved by a modified Gauss-Seidel 
iterative procedure. Currenr values of <^^(t) vjere used as starting 
approximations for (|). (t + At) and the Newton-Raphson method was used 
for finding the physically meaningful -roots . If we denote the current 
root by $. (t + At) , then the new value of c|). (t + At) is calculated from 

4 (t + At) = 4 (t + At) , . + y ('f + ^t) - 0. (t + At) , , 1 (5.2.4) 

where y is the acceleration factor. The next equation in the set, which 
is for A^ , ^ , is solved in precisely the same manner using improved values 
of 4. (t + At) as soon as they ore generated. It has been found that, in 
general, a and y taken together have a decisive influence on the conver- 
gence rate of the iterative process . 

The steady state solution vfas obtained in a similar maimer except, 
of course, that it does not Involve a. By setting the left hand side 
of (5.1.1) to zero and rearranging, we obtain 




where 


, . = B. (c. . + y + y 

1,1- 1' i-j Zj 1-j 


(5.2.5) 


(5.2. 5a) 
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(5.2.5b) 


"1, i 


(kd). 


Z '5-i 


j i 


g, 


0 , 1 



1 

(kd). 

1 



K . <p. 
1-] 


j i 


g.(c. .4% y y ., ^^) 

1 ^ 1 - J J X -j Z—$ I ' - 3 ' 3 ' 

j j ' 


(5.2.5c) 


The iteration bega3i with a suitable set of guessed values of . For 

each node , the Newton-Raphson procedure of successive approximations 

was again used for determining the root $ of (5.2.5). The new value of 

i 

(|). is calculated according to 


= <j>. 

, new 3 


ol d 


+ y($. 



ol d 


•(5'. 2. 6) 


Thus J for either the transient or steady state computations , there are 
ti'io distinct iterations involved. One is associated with the Newton- 
Raphson method of evaluating the roots of (5.2.3) or (5.2.5) and another 
is associated with the Gauss-Seidel procedure. In the former, when 
the successive 's differ by less than 0.01 percent, the result is 
considered satisfactory and the calculation viill proceed to the next 
node. The Gauss-Seidel iteration will terminate vihen the residual of 
each nodal equation is less than a specified value. This residual is 
a measure of the deviation from the perfect heat balance required for 
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each node, Th\is , we may write 


= §0,i + 


for the steady sxate problems and 


HB. = f „ , + f , . (j). + , <i>, 

1 0, i . li 1 1 4, 1 


(5,2.7a) 


(5.2.7b) 


for the transient case. For the hypothetical spacecraft described in 

2 

Chapter 3, . or f^ , is of the order 10 for all sunlit surfaces and 

. —3 

hence , we have arbitrarily set HB, _< 10 as the condition to stop the 
Iteration. This value was chosen essentially from a consideration of 
the computer time on one hand and the desired accuracy of computed re- 
sults on the other hand. 


9 


5.3 DETERMINATION OF OPTIMUM y AND a 

To ascertain the optimum values of y and a for rapid convergence 
of the iterative process, a series of preliminary computer experi- 
ments was conducted. The results for the steady case v^ith normal in- 
cidence from the siin and with a distribution of electronic package heating 
similar to Mariner '6M- spacecraftt are shown in TABLE 5.1. A redistri- 
bution of the heating soin?ce inside the bus or a change in the relative 
importance of the radiative and conductive mode of heat transfer Vfould 
alter the number of iterations required to achieve the same heat balance j 


tThis calls for the following distribution; Aj-29 percent, A^-24 percent, 
Ag-36 percent, and Ag-11 percent. (See Figs. 3.1 and 3.2 for identifi- 
cation). 
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TABLE 5.1 


INFLUENCE OF ACCELERATION FACTOR ON THE CONVERGENCE 
RATE OF THE ITERATIVE PROCEDURE-STEADY STATE 


Acceleration 
factor Y 

No. of iterations 
required to achieve 
HB. < 10-3 

X 

0.3 

>80 

0.5 • 

65 

0.7 

38 

0.9 

24 

1.1 

14 

1.3 

10 

1.5 

’ 17 

1.7 

■ 33 

1.9 

>80 


Note: Normal incidence from the sun; distribution of heat generation 

inside the bus similar to Mariner '64 spacecraft. 
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“3 

i.e., HB. ^10 . However, these changes, if kept within reasonable 

limits, have only a minor effect on the optimum value of y. which is 
around 1.3. This finding is at variance vfith the suggestion made by 
Strong and Emslie [s]. 

The transient results are presented in TABLE 5.2. They were ob- 
tained for the assumed condition that the spacecraft had initially a 
uniform temperature of 400°R and, at an arbitrarily chosen Instant 
h = 0 , the Internal heating over bus surfaces A„ , A, , A„ , and A„ was 
suddenly raised to their full pov/er and, simultaneously, the space- 
i craft's exterior surfaces were exposed. to sun's irradiation. To a 
large extent, the table is self-explanatory. It is seen that the 
acceleration factor has a dominant influence on the convergence rate 
while the integration parameter plays only a minor role. The opti- 
mum value of the acceleration factor has been found to be 1.1 \mder 

i 

the stated conditions . 

The computed temperature data for the 12 nodal surfaces of the 
hypothetical spacecraft at 10, 100 and 300 minutes after the ther- 
mal disturbance are listed in TABLE 5.3. Included are the data for 
the steady state temperatures. As one might expect, the solar 
panels, antenna and the exterior surface of the multi-layer in- 
sulation (A ) respond very rapidly to the changing environmental 
condition. 
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TABLE 5,2 


INFLUENCE OF INTEGRATION PARAMETER AND ACCELERATION FACTOR ON THE CON- 
VERGENCE RATE OF THE ITERATIVE PROCEDURE -TRANS I ENT CASE 


Integration 

Acceleration 

No. of iterations required 
HB. < 10“^ 

to achieve 

parameter 

factor 

Prototype time interval, rnin. 

a 

T 

0-10 

90-100 

290-300 

0.2 

0.9 

5’ 

5 

4 

0.4 

0.9 

5 

'5 

5 

0.6 

0.9 

6 

6 

5 

0.8 

0.9 

6 

6 

5 

1.0 

0.9 

6 

7 

6 . 

0.2 

1.1 

5 

4 

3 

0.4 

1.1 

5 

4 

3 

0.6 

1.1 

5 

4 

3 

0.8 

1.1 

5 

4 

3 

1.0 

1.1 

5 

4 

3 

0.2 

1,3 

9 

7 

6 

0.4 

1,3 

9 

8 

6 

0.’6 

1.3 

9 

8 

6 

0.8 

1.3 

9 

7 

6 

1,0 

1.-3 

9 

7 

6 

0.2 

1.5 

16 

13 

10 

. 0.4 

1.5 

16 

12 

10 

0.6 

1.5 

16 

11 

9 

0.8 

1.5 

16 

11 

9 

1.0 

1.5 

16 

11 

9' " 

0.2 

1.7 

31 

24 

18 

0.4 

1.7 

. 31 

21 

17 

0.6 

1.7 

31 

26 

19 

0.8 

1.7 

31 

24 

18 

1.0 

1.7 

31 

23 

18 

Note: 1. 

Prototype At = 10 min. 



2. 

Normal incidence 

from the sun; 

distribution 

of heat 


generation inside the bus similar to Mariner '6M spacecraft. 
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TABLE 5.3 


THERMAL RESPONSE BEHAVIOR OF THE 
HYPOTHETICAL SPACECRAFT 

Uniform Initial Temperature = 400°R 





Temperature , °R 


t = 10 min. 

t = 100 min. 

t = 300 min. 

Steady State 


518.1 

597.8 

597.9 

598.0 


520.5 

603.8 

604.7 

604. 8 


429 . 4 • 

519.1 

533.9 

534.4 


518.1 

597.6 

597.9 

597.9 


520.2 

602.8 

604.2 

604.3 


412.2 

498. 8 

524.4 

525.5 


401.5 

479.6 

527.7 

529.5 


420.7 

499.5 

524.0 

525.0 

■ 

'415.9 

511.6 

546.2 

— 547.5 

^o 

403.0 

490.8 

527.7 

529.0 


476.3 

568.0 

568.1 

568.1 

^12 ^ 10 '^ 

152.8 

181.1 

184.8 

185.0 


Note: 1. Normal incidence from the sun, solar constant taken to be 

442 Btu/hr-ft^. 

V 

2. Cvi)ical bus consists of' six 4 ft^ surfaces j total electronic 
package dissipation 1100 Btu/hr distributed over surfaces 
A 3 (29 percent), Ae (24 percent), Ae (36 percent), and A 3 (11 per- 
cent ) . 
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Chapter 6 


RESULTS, DISCUSSIONS AND CONCLUSIONS 

The similitude requirements for the thermal modeling of the hy- 
pothetical spaceez'’aft have been given in Section 3.4. In the present 
study, it vras assumed that the models entail errors in the desired 
values of kd, Cd, e, and e’ (the prime refers to the interior of the 
surfaces). There was no error arising from the requlremenx of geo- 
metric similarity, nor was there error in the heat fluxes q and q*. 

As has been noted previously, this choice of error distribution is 
based primarily on a consideration of practical model fabrication 
with the present-day technology. It has no bearing on the theory 
itself. 

The general theory of imperfect modeling presented in Chapter 
2 is described logically in terms of dimensionless quantities. How- 
ever, in the present instance , it is possible to speak of errors in 
the dimensional quantities like (kd) and (Cd) since not only are they 
the only quantities which may be varied in their respective dimension- 
less groups, but also they are physically apprehended more readily 
by the designer. 

6.1 SPACECRAFT SUBJECTED TO SIMPLE HEATING AND COOLING 

To ascertain the usefulness and liirdtations of the proposed theory , 
several sets of errors in the modeling parameters have been assigned. 
TABLE 6.1 lists the values for the first of such sets. In the table 
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TABLE 6.1 Errors "in modeling parameters 



Percentage Error in Modeling Parameters 

Surface 

Model MI 

Model Mil 

Model Mill 

‘^kd 

X 100 

dcd 
X 100 

X 100 

6gi 

X 100 

djtd 
X 100 

dcd 
X 100 

X 100 

X 100 

X 100 


X 100 

X 100 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

.0 

0 

0 

0 

0 

0 

0 

0 

0 


25.0 


-18.0 

10.0 


52.5 


20.0 

25.0 

25.0 

-28.8 



0 

Hi 

0 

0 

0 

0 

0 

0 

0 

0 

• 

0 

■ 


m 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 ' 

0 


25.0 

29.0 

- 20.0 

10.0 

45.0 

50.0 

- 20.0 

20.0 

25.0 

20.0 

-32.0 

10.0 

K 

i 

20.0 

25.0 

15.0 

- 8.0 

40.0 

55.0 

27.0 

- 8.0 

20.0 

25.0 

15.0 

-14.7 

\ 

25.0 

25.0 

-18.0 

- 8.0 

57.5 

57.5 

-18.0 

- 8.0 

25.0 

25.0 

-25.9 

-11.5 

\ 

25.0 

25.0 

25.0 

- 8.0 

45.0 

45.0 

45.0 

- 8.0 

25.0 

25.0 

25.0 

-16.0 

\o 

20.0 

25.0 

15.0 

- 8.0 

44.0 

60.0 

33.0 

- 8.0 

20.0 

25.0 

15.0 

- 12.8 

• 

0 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

0 

^ 12 ^”^10 ' ^ 

0 

0 

15.0 

-- ----- 

0 

0 

0 

33.0 

0 

0 

0 

15.0 

0 


6 ^. refers to the exterior of the surface in question and 6 ^, 


refers to its interior' 































































and others which follo'w, inode-l MI refers to the 'original' imperfect 
model. Mil refers to that for which the positive errors are modified 
(in this case, they become moTe positive") while the negative errors 
remain unaltered and Mill denotes that- for which the negative errors 
are modified (in this case, they become move negative) while the posi- 
tive errors are unchanged. Thus, the error state of model MI would 
he represented by the point P in the multidimensional hyperspace re- 
ferred to in Chapter 2. Likewise, the- error states of model Mil and 
Mill are represented respectively by and P*^ ^ . 

Initially, the spacecraft is assumed to have a uniform tempera- 
ture of 400°R. At some arbitrarily selected instant, t = 0, inter- 
nal heating within the spacecraft's bus and of a distribution as pre- 
scribed in TABLE 5.3 is assumed to begin instantaneously and, at the 
same instant, its exterior surfaces are exposed to sun's irradiation. 
Subsequent to the establishment of a steady thermal condition , the 
spacecraft is conceived to move into the shade of a planet, thus cut- 
ting off the sun's irradiation. In this manner, a cooling transient 
is created. 

Using the equation set (3.3.3) and the computational procedure 
explained in Chapter 5 , temperature data for the heating transient , 
the steady state, and the cooling transient are obtained for the pro- 

V 

totype (or, equivalently, the perfect model) and for the three models. 

The results are summarized in a table set 6.2 for normal incidence 
of the solar irradiation and in table set 6.3 for 45° oblique inci- 
dence. The prototype temperatures are denoted by Tp and those of 
the models by , and • The errors in the latter temperatures 
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TABLE 6.2(a) Model temperature before and after correction 


Heating Transient: 10 minutes after spacecraft subject to solar irradiation (normal incidence) and internal beating 


Surface , 

Protot3Tpe 
or Perfect 
Model 

Temperature 

■ "p 

■ " ' ' ' ' 

Temperatures of Imperfect Models 

— 


A 

t 

Model Temperature 
after Correction 

_ ^5odel MI 

Model Mil 

Model Mill 

T 

iva 


T 

MI I 


T 

MI I I 

1 1 

T 

MC 



518.1 

518.1 

0 

518.0 

- 0.1 

518.0 

- 0.1 

- O-.Ol 

- 0.05 

WBBI 

' 518.1 

0 


520.5 

520.0 

- 0.5 

520.0 

- 0.5 

519.8 

- 0.7 

- 0 . 05 , 

- 0.32 

- rO .37 

520.4 

- 0.1 

A. 

1 + 29.4 

426,2 

- 3,2 

421.9 

- 7,5 

427.5 


- 3.89 

1.88 

- 2 . or 

428,2 

- 1.2 ■ 


518.1 

• 518.-0 

- 0.1 

518.0 

. - 0.1 

518.0 

- 0.1 

’""0 

- 0.04 

- 0.04 

518.1 

0 

- \ 

520.2 

520.0 

- 0.2 

519,9 

- 0.3 

519.7 

- 0.5 

- 0.2 

- 0.25 

- 0.27 

520.2 

0 


412.2 

411.4 

- 0.8 

409.2 

-3, a 

412.0 

- 0.2 

■ - 1.95 

0.98 ' 

- 0.97 

412 . 3 ' 

0.1 


401.5 

400.9 

- 0.6 

400.6 

- 0.9 

400.9 

- 0,6 

• 

- 0.29 

- 0.11 

- 0.40 

401.3 

- 0.2 


420.7 

419.0 

- 1.7 

415.4 

- 5.3 

420.0 

- 0.7 

- 3.23 

1.37 

- 1.86 

420.9 

0.2 


415.9 

412.9 

- 3.0 

410.9 

- 5.0 

413.0 

- 2.9 

rl .76 

0.23 

- 1.53 

414,4 

- 1.5 


403.0 - 

401.8 

- 1.2 

401.0 

- 2.0 ■ 

401.6 

- 1.4 

- 0.66 

- 0.29 

- 0.95 

402.7 

- 0.3 


476.3 

476.3 

0 

476.3 

0 

476.3 

0 

-0 

~0 

"0 

476.3 

0 


152.8 

is 3.3 

0.5 

153.2 

0.4 

154.1 

1.3 

. - 

- 0 , 07 ' 

1.18 

1.11 • 

152.2 

- 0.6 


All temperatures in °R 

























































TABLE 6.2(b) Model temperatures before and after correction 


Heating Transient: 30 minutes after spacecraft subject to solar irradiation (normal incidence) and internal heating 


Surface 

Prototype 
or Perfect 
Model 

Temperature 

^ 

Temperatures of Imperfect Models 



■ \ 

Model Temperature 
after Correction 

I ; 

Model MI 

Model Mil 

Model Mill 

T E 

M 1 M 

T 


T 

M 1 1 

E 

Ivl 1 1 

T 

MC 

^!V7C 

^ ■ 

590.7 

590.6 

-0.1 

590.6 

-0.1 

590.5 

-0.2 

-0.04 

-0.07 

-0.11 

590.7 

0 


59i^.9 

594.2 

-0.7 

593.9 

-1.0 

593.9 

-1,0 

-0.22 

-0,43. 

-0.65 

594.8 

-0.1 

A„ 

474.5 

469.0 

-5,5 

459.4 

-15.1 

472.0 

■-2.5 

-8.61 

4.47 

-4.14 

473.1 

-1.4 


590.5 

590.5 

0 

590.5 

0 

590.5 

0 

■ -0.02 

-0.04 

• 

• -0.06 

590.5 

0 


593.9 

593.5 

0 

593.4 

-0.1 

593.4 

-0.1 

-0.11 

■-0.20 

-0.31 

■ 

593.9 

0 


440.3 

436.8 

-3.5 

.430.3 

-io.o 

438. 3 

-2.0 

-5.83 

Bi 

-3.59 

440.4 

■ 0.1 

\ 

413.1 

409.3 

-3. 8 

406,5 

-6.6 

409,1 

-4.0 

-2.52 

-0,30 

-2.82 

412.1 

-1„0 


450.6 

448.9 

BB 

441.0 

-9.6 

451.6 

1.0 

-7.10 

4.09 

■ 

-3. 01 

451.9 

^g| 

\ ■ 

445 . 4 

438.7 

B 

432.8 

-12.6 



-5.24 ■ 

1.37 


442.6 

-2.8 I 


424. 3 

418.3 

- 6.0 

412.9 

-11.4 

418.3 

-6.0 

-4.79 

-0.01 

-4.80 

423.1 

-1.2 


547.6 

547.6 

0 

547,6 

0 

■ 547.6 

0 

0.01 


■■ 

547.6 

0 1 


172.2 

172,3 

B 

171.9 

-0.3 

172,9 

0.7 

-0.36 

0.87 

.0,51 

171.7 

-0.5 j 


All temperatures in 














































































TABLE 6.2(c) Model temperatures before and after correction 


Heating Transient; 100 minutes after spacecraft subject to solar irradiation (normal incidence) and internal heating 


Surface 

Prototi^e 
or Perfect 
Model 

Temperature 

p 

Temperatures of Imperfect Models 


' a'*' 


Model Temperature 
after Correction 

Model MI 

Model Mil 

Model Mill 

T 

MI 

T ■ 

Mt I 


T 

M 1 1 

E 

MI r I 

■ T 

lAC 


^ ■ 

597.8 

597.7 

-0.1 

597.6 

-0.2 

537.6' 

-0,2 

-0.06 

-0.08 

-0.14 

597.8 

0 


603.8 

603.1 

-0.7 

602.6 

-1.2 

602.7 

-1.1 

-0.36 

-0.47 

-0.83 

603.9 

0.1 


519.1 

524.3 

5.2 

515 . 0 

-4.1 

532.5, 

13.4 

-8.35 

•12.24 

3.89 

520.4 

1.3 


597.6 

597.5 

-0,1 

597.4 

-0.2 

597.5 

-0.1 

-0.06 

t0.06 

-0.12 


0 


602,8 

602.1 

-0.7 

601,7 

-1.1 

601.9 

-0.9 

-0.37 

. 

-0.33 

-0.70 

■ 

■ 602.8 

0 


49 8. 8 

499.8 

1.0 

488,0 

-10,8 

505.9 

■ 7.1 

-10.52 

9.13 

-1.39 

501.1 

2.3 


479,6 

469.4 

-10.2 

454.1 

-25.5 

471.1 

-8.5 

-13.68 

2.55 • 

-11.13 

480.5 

- 0.9 


. 493.5 

505.0 

• 5.5 

494.8 

BB 

512.3 

12,. 8 

-9.14 

10.81 

1.67 

503.4 

3.3 


511.6 

506.9 

mn 

494.2 

-17.4 

511.9 

0.3 

-11.36 

7.51 

-3.85 

510.8 

-0,8 


490.8 

486.2 

BB 

472.5 

-18.3 

490,0 

-0.8 

-12.25 

5.79 

-6.46 

492.6 

1.8 


568.0 

568.1 

0.1 

568.1 

0.1 

568.1 

0.1 

~0 

~0 

0 

568.0 

0 

^12 0 ' ^ 

181,1 

181.6 

0.5 

180.2 

-0.9 

183.0 

B 

-1.18 

2.22 

1.04 

■ 180,5 

-0,6 


All temperatures in 


























































































































































TABLE 6.2(d^ Model temperatures before and after correction 


Heating Transient: 300 minutes after spacecraft subject to solar irradiation _(normal_ incidence) and internal heating 


Surface ■ 

. 

Prototype 
or Perfect 
Model 

Temperature 

Tp 

Temperatures of Imperfect Models' 




ri 1 * MT-n 

Model Temperature 
after Correction 

Model MI 

Model Mil 

Model Mill 
- 

• 



T 

E 

r/2 

T 

MI I 

I 

T 

"MI I I 

^IVIE I I 

T 

MC 

j 

^ • 

598.0 

597.9 

-0.1' 

597.9 

-0.1 

597.9' 

-0.1 

-0.03 

-0.05 

-0.08 

598.0 

0 1 


604.7 

604.4 

-0.3 

604.2 

-0.5 


-0.5 

-0.15 

-0.30 

• -0.45 . 

604.9 

0.2 1 


533.9 

550.7 

16.8 

547.6 

13.7 

563.6 

29.7 

-2.85 ■ 

19.22 



0.5 1 


597.9 

597.8 ■ 

-0.1 ■ 

597.8 

-0.1 

597.8 

-0.1 

-0.03 

-0.03 

HI 


0 j 


604. 2 

604.1 

-0.1 ' 

603.9 

-0.3 

603.9 

-0.3 

-0.17 

-0.19 

-0.36 

604.4 

0.2 


,524.4 

542.1 

17.7 

538.3 

13.9 . 

555.4 

31.0 

-3.40 

• 19. 89 , 

16.49 

525.6 

1.2 


527.7 

540.6 

12.9 

533.8 

6.1 

550.5 

■ 22.8 

-6.14 - 


14.68 

8.54 

532.1 

4.4 

1 


524.0 

543.0 

19.0 

539.2 

15. '2 

555.3 

31.3 


18.32 

14.86 

528.2 

4.2 


546.2 

558.3 

12.1 

551.7 

5.5 ■ 

569.6 

23.4 

-5.93 

16.84 



1.2 

^o 

527.7 

541.2 

13.5 

535.5 

7.3 

551.3 

23.6 

-5.11 

15.04 

IIIIIIQIIIIIIIIIIII 


3.6 


568.1 

568.2 

0.1 

568.2 

0.1 

568.2 

0.1 

0.02 

0.02 

0.04 

568.1 

0 


184.9 

187.7 

2.8 

186.9 

■ 2.0 

190.6 

5.7 

-0.63 

4. '42 

3.79 

183.9 

-1.0 

f 


All temneratures in °R 





























































































































TABLE 6.2(e) Model temperatures before and after correction 


Heating Transient: 500 minutes after spacecraft subject to solar irradiation (normal incidence) and internal heating 


Surface 

Prototype 
or Perfect 
Model 

Temperature 

T 

P 

Temperatures of Imperfect Models 


A^-> 

^ . 

Model Temperature 
after Correction 

Model MI 

Model Mil 

Model Mill 

T 

E 

-M 

T ■ 
m ! 



E 

m. X I 

T 

MC 



598.0 

h9 

-0.1. 

597.9 

-0.1 

597.9 

-0.1 

-0.01 

-0.05 

-0.06 


0 


604.8 

604.5 

-O'. 3 

604.5 

-0.3 

604.4 

-0.4 

-0.05 

-0.27 . 

-0.32 

604.8 

0 


534.4 

552.8 

18.4 

551.7 



31.9 

-1.00 

20.11 

19.11 

532.7 ■ 

-1.7 . 


597.9 

597.9 

0 

597.9 

0 

597.8 

-0.1 

■-0 . 01 

-0.03 

-0.04 

597.9 

0 

• \ 

604.3 

604.2 

-0.1 • 

604.2’ 

-0.1 

604.1 

‘ 

-0.2 

-0.04 




0.1 


525.3 

545.3 

20.0 

544.4 



34.2 

-0.77 

21.23 

20.46 

523.5 

— 

l.,8 


. 529.4 

545.9 

16.5 

543.8 



27.5 

-1.89 

16.37 



0.7 

. ^.3 

525.0 

546.1 

21.1 

544.9 



IQQIII 

-1.04 

19.40 

18.36 

526.9 

1.3 


547.5 

562.3 

14.8 

559.1 



27.0 

-2.93 

18.19 

15.26 

546.8 

-0.7 ■ 

^0 

529.0 

545.5 

■ 

16.5 

543.5 

14.5 

556.4 

27.4 

-1.78 

16,38 ■ 

14,60 

529.7 

0.7 


568.1 

568.2 

o.i 

568.2 

,0.1 

568.2 

0.1 

0.02 

0.02 

0.04 

568.1 

0 


185.0 

188.2 

3.2 

187.9 

2.9 

191.3 

6.3 

-0.23 

4.70 

4.47 

184.4 

-0.6 


All temperatures in ®R. 

























































































































TABLE 6.2(f) Model temperatures before and after correction 


steady State 


Surface , 

Prototype 
or Perfect 
Model 

Temperature 

■ "P 

Temperatures of Imperfect Models 



' 

Model Temperature 
after Correction 

Model MI 

Model Mil 

Model Mill 

T 

m 

E 

M 

T 

m I 


T 

m 1 1 

^IVII II 

. \c 

V 

vc 


597. S 

597.3 

0 

597,9 

0 

597.9 

0 

-0.01 

-0.05 

-0.06 

598.0 

0.1 


604.8 

604. 5 

-0.3 

604.5 

-0.3 

604.4 

-0.4 

■ -0.03 ■ 

-0.26‘. 

-0.29 

604.8 

0 


534.4 

553.0 

18.6 

552.2 

17.8 ■ 

566.5 

32.1 

*"0 -68 

20.20 

19.52 

533.5 

-0.9 


597.9 

597.9 

0 

597.9 


597.9 

0 

~0 

-0.02 

-0-. d] 

597.9 

0 


604.3 

604.2 

-0.1 


-0.1 

604,2 

-0.1 

-0.01 

-0.14 

-0.15- 

604.4 

= 0.1 


525.3 

54;5.5 

20.2 

545.2 

19.9 

559.8 

34.5 

-0.31 

21.36 

_ 21.05 

524.4 

-0.9 

. ^ 

529.5 

546.3 

16, 8 

545.0 

15.5 

,557.4 

27.9 

-1.15 


15. 39 

530.9 

1.4 


525.0 


21.3 

545;6 

20.6 

559.4 

34.4 

t0x60 

19.51 

18.91 

527.4 

2.4 


5.47.5 

562.6 

15,1 

560.0 

12.5 

574.9 

Mgm 

-2.38 

' 18.33 

15.95 

* 

546,7 

-0.8 

^o■ 

529,0 

545.8 

16. 8 

544.5 

15.5 

556.8 

27,8 

• 

-1.19 

16.52 

15.33 

530,4 

1.4 


568.1 

568.2 

0.1 

568, 2 

0 .1 . 

568.2 

0.1 

0.02 

0.02 

0.04 

568.1 

0 


185.0 

■HI 

3.2 

188.0 

•3.0 

191.4 

6.4 

-0.16 

4.73 

4. 57 

183.6 

-1.4 


All temperatures in '^R. 














































































TABLE 6.2(g) Model temperatures before and after correction 
Cooling Transient: 30 minutes after spacecraft moving dnto shade 


Surface 

' Prototype 
or Perfect 
Model 

Temperature 

Tp 

Temperatures of Imperfect Models 


t 


Model Temperature 
after Correction 

^ ^5odel MI 

Model Mil 

Model Mill 

T 

"M 


T 

■'MI I 


T 

Mt I I 

I I 

T 

MC 



376.5 

376.5 

0 

376.5 

0 

376.4 

-0.1 

■ 0.01 

-0.12 . 

-0.11 


0.1 


390.5 

390.1 

-0.4 

390.2 

-0.3 

389.7 

-0.8 

0.03 

-0.66 

-0.63 

390.8 

0.3 


516.2 

539.6 

23.4 

540.4 

24.2 

554.7 

38.5 

0.78 

22.57 

23.35 

516,2 

0 


376.4 

376.3 

-0.1 

376.4 

0 

376.2 

-0.2 

0.01 

-0.14 

-0,13 

376.5 

0,1 


389.9 

389.5 

-0.4 

389.5 

- 0 . 4 ’ 

388.9 

0 

0.05 

-0.80 

-0.75 

m 

.0.3 


513.0 

536.8 

23.8 

537.9 

24.9 

552.3 

39.3 

0.94 

23.14 

24.08 

512.7 

-0.3 


■ 527.0 

544.5 

17.5 

543.6 

16.6, 

555.7 

28.7 

-0.81 

16.73 

15.92 

528.6 

1,6 

, 


522.6 

544.5 

21.9 

544.2 

21.6 

557. 8 

35.2 . 

-O'. 34 

19.78 

19.44 

525,1 

2.5 


543.4 

559.6 

16.2 

557.4 

MBHBl 

572.2 

28.8 

-2.02 

18.80 

16.78 

542.8 

-0.6 


525.5 

543.2 

17.7 

542.4 

16.9 

554.5 

29.0 

-0.75 

16.86 

16.11 

527.1 

1.6 


395.2 

395.3 

0.1 

395.3 

0.1 

395.3 

imQii 

0.04 

0.04 

0.08 

395.2 

0 

SMB 

150.4 

156.5 

6.1 

156.3 

5.9 

162.1 

11.7 

-0.18 

8.37 

8.19 

148.3 

-2.1 


All temperatures in °R 


















































































































































TABLE 6.2(h) Model temperatures before and after correction 
Cooling Transient: 200 minutes after spacecraft moving into shade . 




Temperatures of Imperfect 

Models 

— 



\ 

• 


Surface , 

Prototype 
or Perfect 
Model 

- ; Model MI 

Model 

Mil 

' 

Model Mill 



\ . 

Model Temperature 
after Correction 


Temperature 

t 



T 

MI I 

I 

T 

MI I I 

'^Ivl II 

■ 


T 

MC 



229.8 

229.7 

- 0.1 

229.8 

0 

, 

229.5 

- 0.3 

0.07 

- 0.36 

- 0.29 

230.0 

0.2 


275.7 

275.4 

- 0.3 

275.7 

0 

274.5 

- 0.2 

0.25 

- 1.33 

- 1.08 

275.5 

0.8 


500.3 

524.6 

24.3 • 

525.1 

24.7 

540.9 

40.6 

0,45 ■ 

24.46 

24.91 

499.7 

- 0.6 , 


229.0 

228.9 

- 0.1 

229.1 

0.1 

228.5 

- 0.5 

0.14 

. - 0.58 

- 0.44 

229.4 

0.4 


272.7 

272.3 

- 0.4 

272.9 

. 0-2 

270.9 

0 

0.50 

- 2,14 

- 1,64 

274.0 

1.3 


491 . 8 - 

518.4 

26,5 

519.6 

27 . 8 

535.6 

43,6 

1.11 

• 25.80 

26.91 

491,4 

- 0.4 


• 508.8 

528.6 

19.8 

528.3 

19.5 

541.3 

32.5 

- 0.24 

19.01 

18.77 

509.8 

1.0 


510 . 4 j 

'f 

533.2 

22.8 

532. '8 

22.4 

547.4 

37.0 

- Q .40 

21.23 

20.‘83 

512.4 

2.0 


525 . e ' 

543.6 

18.0 

541.6 

16.0 

557.7 

32.1 

- 1.83 

21.13 

19.30 

524.3 

- 1.3 

^o 

508.3 

527.8 

19.5 

527.2 

18.9 

540.6 

32.3 

- 0 . 54 - 

19.09 

18.55 

509.3 

1,0 


235.6 

236.0 

0.4 

236.1 

0,5 

236.1 

0.6 

0.16 

0.20 

0.36 

235.6 

0 


134.8 ■ 

142.7 

. 7.9 

142.6 

' 

7.8 

149.8 

15.0 

- 0.15 

10,51 

10.36 

132.4 

- 2.4 


All temperatures in °R. 




TABLE 6.3 (a) Model temperatures before and after correction 
Heating^ Transient: 10 minutes after spacecraft subject to solar irradiation (45® incidence) and internal heating 


Surface 

1 

Prototype 
or Perfect 
Model 

Temperature 

• 

Temperatures of Imperfect Models 

' 

‘ . 

V 

'f 

■ 4 

Model Temperature 
after Correction 

■ 

Model MI 

Model Mil 

Model Mill 

m 

*14: 

X* 

■"mi 

T 

MI I 

1- 

T 

M 1 1 

r I 

T 

MC 

■ ' j 

' i 

477.4 

477.3 

- 0.1 

477.3 

-0.1 

477.3 

-0.1 

- 0.01 

-0.07 

mm 

477.4 

0 


370.1 

369.5 

- 0.6 

369 . .5 

-0.6 

369.3 

-0.8 

-0.06 

-0.44 

-0.50 

370.1 

0 ■ 

1 ^ 

425.7 

423.7 

- 2.0 

419.8 

-5.9 

425.3 

-0.4 

-3.53 

2.36 

-1.17 

424.9 

- -^o.-s :i 

— 

481.0 

481.0 

0 

481.0 

0 

480.9 

-0.1 

0 

0 

1 



481.0 

! 

0 f 

\ 

. 499.3 

49 8. 8 

- 0.5 

498,8 

-0.5 . 

498.6 

-0.7 

- 0 . 06 

-0.33 

, 

-0 . 3^: 

499.2 

i 

( 

- 0.1 


430,4 

426.9 

-3.5 

421.8 

-8.6 

427.7 

-2.7 

-4.57 

1.19 

CO 

CO 

CO 

1 

430.3 

- 0.1 ! 

“7 

401.7 

401.1. 

-0.6 

400.7 

-1.0 

401.0 

-0.7 

-0.34 

-0.11 

-0.45 

.401.6 

- 0.1 ! 

■ \ 

421.3 

419.5 

^’ 1.8 

415.8 

-5.5 

420.5 

-0.8 

-3.33 

1,40 

-1.93 

421.5 

0.2 


416.9 

9 ♦ 

413.8 

msM 

,411.7 

-5.2 

414.0 

-2,9 

-1.93 

0.27 

— 

415.5 

-1.4 

^0 

403.6 

402,2 


401.4 

-2.2 

402.1 

-1.5 

• 

-0.78 

, -0.26 

-1.04 

403.3 

-0.3 


447.6 

447 . 6 

, 0 

447.6 

0 

447.6 

g 

0.01 

0 

■ Q.Ol 

447.6 

0 


142.1 ■ 

142.9 

-0.8 

142.7 

-0.6 

143.9 

- 1.8 

•-0.10 

1.60 

1.50* 

141.4 

-0.7 


All temperatures in ®R 























































































































































TABLE 6.3(b) Model temperatures before and after correction 


Heatins Transienl:: 30 minutes after spacecraft subject to solar irradiation (45® incidence) and internal heating 


— 

Surface 

1 

Prototype 
or Perfect 
Model 

Temperature 

Tn 

p 

^ 

Temperatures of Imperfect Models 


r * 

' A <- > 

. 1 

S Model Temperature 
after Correction ' 

Model MI 

Model Mil 

Model Mill 

rr> 

r * 

T 

Ml I 


T 

M ! 1 1 

E 

IvE 1 1 

T . 
.Vie 

’ 1 


535.6 

535.4 

- 0.2 

535.4 

■EBI 

535.4 

- 0.2 

warn 

- 0.12 

- 0.16 

535.6 

( 

0 1 


335.3 

333.6 

- 1.7 

333 . a 

- 2.2 

, 332.8 

- 2.5 

- 0.37 

- 1.16 

- 1.53 

335.1 

- 0 _^ ...1 


457.7 

456.9 

- 0.8 

448 . 8 

- 8.9 

461.4 

. 3.7 

- 7.18 

-, 6..-79 

- 0.39 

■ 457,2 

i 

,1 

. - 0.._5 _l 


541.9 

541.8 

- 0.1 

541.7 

- 0.2 

541.7 

- 0 * 2 . 

-- 0.06 

- 0.09 

- 0.15 

541.9 

! 


568. 3 ■ 

567.4 

- 0.9 

567.1 

mm 

567.1 

- 1.2 

'- 0.30 

, - 0.48 

- 0.78 

■ 

568.2 

0 1 

- 0.1 1 


480.6 

47 3 .. 2 

- 7.4 

460.8 

- 19,8 

475,7 

1^9 

- 11;18 

, 

3.64 

- 7.54 

480.8 

i 

0.2 ! 

1 

414*9 

410 . 7 . 

- 4.2 

407.5 

MR 

410.6 

- 4.3 

- 2.92 

- 0.23 

- 3.15 

413.9 

! 

- 1.0 j 

■ 

453,5 

451,7 

- 1.8 

443,2 

- 10.2 

454.6 

■ 1.1 

- 7.57 

4.38 

- 3,19 

454.9 

1.4 1 

! 

■ , 450.5 

443.6 

- 6.9 

436,8 

- 13.7 

444.8 

5.7 

.- 6.08 

1.82 

- 4.26 

447.8 

t 

- 2.7 1 


427.6 

421,3 

- 6.3 

415,1 

- 12.5 

421.4 

- 6.2 

- 5.54 

0.29 

- 5.25 

426.5 

- 1.1 1 


498.0 

498.0 

0 

498,0 

0 ' 

498,0 

0 

' 0.01 

~0 

0.01 

• 498,0 

0 { 


155.5 

155.6 

0.1 

155.0 

■ ^ 0,5 

156.6 

1.1 

B 

1.47 

0.90 

154.7 

- 0.8 j 


All temperatures in ®R. 





































































































































TABLE 6.3(c) Model tempev'atures b.efore and after correction 

/ 

Keating Transient: 100 minutes after spacecraft subject to solar irradiation (45° incidence) and internal heating 


Surface 

1 

Protot3TDe 
or Perfect 
Model 

Temperature 

Tp 

Temperatures of Imperfect Models 


j 

■A 

Model Temperature 
after Correction 

■Model MI 

Model Mil 

Model Mill 

rp 

E 

MI 

T 

MI I 


T 

T/t 1 1 

1 1 

T 

MC 



546.3 

546.1 

-0.2 

546.0 

-0.3 

546.0 

-0.3 

' -0.08 

-0.13 

-0.21 

546.3 

0 


296.2 

292.3 

-3.9 

290.3 

-5.9 

290.6 

-5.6 

-1.73 

-2.51 

-4.24 

296.5 

0.3 


498.1 

508.9 

10.8 

500,2 


520.0 

21.-3 

-7.81 

16.67 

8.86 

500,0 ■ 

1.9 


553.4 

553.1 

-0.3 

553.0 

.-0.4 



-0.15 

-0.13 

-o ;28 

553.4 

0 


581.6 

580.2 


579.4 

-2.2 

579.8 

■g 

-0.78 

. -0.68 

-1.46 

581.7 

0,1 


552.8 

556. '0 

3.2 

540.3 

-12.5- 

565.7 

— 

12.9 

■Bl 

14.45 

0.35 

555.7 

2,9 

. \ 

490.6 

481.1, 

-9.5 

463.2 

-27.4 . 

484.2 


-16.04 

4.55 

-11.49 

492.6 

2.0 1 


507.9 

515.1 

7.2 

504.0 

-3.9 

523.7 

15.8 

-9.99 

' 12.83. 

■2.84 

512,3 

4.4 


524.9 

522.7 

-2.2 

508.3 

-16.6 

529.7 

4.8 

-.12.84 

10.49 

-2.35 

525.0 

0.1 


502.6 

499.7 

-2.9 

484.1 

-18.5 

505.4 

2.8 

-14.00 

8.45 

-5.55 

505.3 

2.7 

A I 

520.6 

520.7 

0.1 

520.7 

0.1 

520.7 

0.1 

~0 

0.01 

0.01 

520.7 

0.1 


167.2 

168.1 

0.9 

166.0 


170.6 

B1 

-1.85 

3,69 

1.84 

166.3 

IIIIIIIIQQIIIIIIIIII 


All temperatures in °R. 

































































































































TABLE 6.3(d) Model temperatures before and after correction 

Heating Transient; 300 minutes after spacecraft subject to solar irradiation (45® incidence) and internal beating 


Temperatures of Imperfect Models 


Model Mill 


Surface 



Prototype 
or Perfect 
Model 

Temperature 

T 

p 


546.5 


290.8 


515.2 


553.8 


583.4 


573.4 


■ '540.1 


533.1 


540,0 


520.9 



•Model MI 

Model 

Mil 

T 

1 IV2 

E 

M 

T 

MI I 


546.4 

-0.1 

546.4 

-0.1 

289.8 

-1.0 

288.6 

-2.2 

539.9 

24.7 

538.7 

23.5 

553.7 

-O.l' 

553,6 

.-0.2 

582.9 

-0.5 

582.5 

-0,9 

598.6 

20.2 

588.4 

15.0 

557.6. 

17,5 

551.5 

11.4 

• 555.6 

22.5 

552.5 

19.4 

576.3 

16.7 

5 70 .'3 

10.7 

5.58.0 

18.0 

553.0 

13.0 

1 521.0 
1 

0.1 

521-.0 . 

0.1 


All temneratures in ®R 


176.7 


■175.9 


MU ivnn 




181.2 I 9.0 


- 0.02 


-1.08 


-1.04 2 


^0,07 -0.08 


-0.33 I -0.38- -0.71 


•20.23 


-4.62 j 24.85 


-5.42 18.60 


- 2..11 21.38 


-5.34 '• 20,76 


-4.46 18.87 


0.02 


7 ^ 0.73 


13.18 


18.61 


15.42 


Model Temperature 
after Correction 



553.8 




573.4 


544.4 


537.0 


560.9 


543.6 


520.9 


170.7 I -1.5 
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TABLE 6.3(e) Model temperatures before and after correction 


Heating Transient: 500 minutes after spacecraft subject to solar irradiation (45° incidence) and internal heating 


Surface 
1 ' ' 

'Prototype 
or Perfect 
Model 

Temperature 

T 

P 

Temperatures of Imperfect Models • 

-A(-) 

. 

A<-> 

A .. 

I Model Temperature 
j after Correction 

'Model MI 

Model Mil ‘ ■ 

Model Mill 


' 

T 

Mil 


T 

Mi 1 1 

■^Ivl 1 1 

T 

MC 



545.5 

546.5 

0 

546.5 

0 

546.4 

-0.1 

0.01 

-0.06 

-0.05 

54.65 

0 


291.4 

291.7 

0.3 

291.9 

0.5 

291.0 

-0.4 

0.19 

. 

-1.05 

-0.86 

29.26 

1.2 


515.7 

541.7 

26.0 

542.7 

27.0 

558.7 

43.0 

0.85 

25.41 

26.26 

515.5 

-0.2 i 


553,8 

553.7 

-0.1 

553.7 

-0.1 

553.7 

-0.1 

-0.03 

-0.07 

■■ 

553.8 

1 


583.4 

. 

583.1 

-0.3 

582.9 

-0.5 

582.8 

-0:6 

-0.18 

-0.34, 

-0.52 

583.6 

0.. 


573.9 

595.7 

21,8 

592.9 

19.0 

612.8 

38.9 

-2.50 

25.66 

23.16 

572.5 

-1.4 

’ ■ 1 

541.3 

561.8 • 

20.5 

560.2 

18.9 

575.0 

33.3 

-1.40 

19.79 

18.39 

543.4 

! 

2.1 { 

^ ■ 

533.7 

558.1 

24.4 . 



572.9 

39.2 

mmm 

22,16 

21.69 

536.4 

2,7 


560.5 

579.4 

18.9 

576.6 

imm 

593.9 

33.4 

-2.53 

21.69 

19.16 

560.3 

-0.2 

^0 

541.0 

561.3 

20.3 

559.8 

18.8 

574.6 

33.6 

-1.33 

. 19,81 

18.48 

542,8 

1,8 

A: 

520.9 

521.0 

.0.1 

521.0 

0.1 

521.0 

0.1 

0.03 

0.04 

.0.07 ’ 

520.9 

0 


172.3 

177.2 

4.9 ' 

177.0 

' 

mm 

181.9 

9.6 

-:0.23 

6.96 

6.73 

170.5 

-1,8 


All temperatures in °R. 









































































TABLE 6.3(f) ^iodel temperatures before and after, correction 


CO 


Steady State 


Surface 

f 

t 

Prototype 
or Perfect 
Model 

Temperature 
• Tp • 

Temperatures of Imperfect Models 

■'a'"’ 


■ • 

Model Temperature 
after Correction 

•Model MI 

Model Mil 

Model Mill 

rn 

*r/3 


T 

m I 


T 

M 1 1 

E 

lvn 1 1 

T 

,ViC 

^r/jc 1 


546.5 

546.5 

0 

546.5 

0 

546.4 



-0.06 

-0,04 


■ 0 1 


231.5 

291.9 

0.4 

292.4 

0,9 

291.3 

-0.2 

0.45 

-0.98 

-0.53 

292,5 

1.0 1 

A, 

515.8 

541.9 

26.1 

543.1 

27.3 

558*9 

43 .■! 

1.10 

25.46 

26.56 

515.3 

i 

-0.5 { 


553.8 

553.7 

-0.1 

553.7 

-0.1 

553 '. 7 

-0.1 . 

-0.03 

t0.07 

■ -0.10 

553.8 

1 

0 


583.4 ■ 

583.1 

-0.3 

582.9 

- 0:5 

582.8 

-0.6 

-0.15 

-0.34 

-0.49 

583.6 

0.2 

r 

573.9 

595.. 8 

21.9 

593.3 

19.4 

613.0 

39.1 

-2.23 

25.72 

23.49 

572.3 

-1.6 



562.0 

20.6 

561,0 

19.6- 

.575.3 

33.9 

-0.88 

19.87 

18,99 

543,0 

1.6 1 


533.8 

558.2 

24.4 

558.0 

24.2 

573.1 

39.3 

-0.’17- 

22.21. 

22.04 

536.2 



j||||||SQ^^ 

579.6 

19.0 

577.2 

16.6 

594.2 

33.6 

-2.15 

21.75 

19.60 

560.0 

70,6 

^o■ 

541.0 

561.5 

20.5 

560.5 

19.5 

574.8 

33.8 

-0.92 

19.87 

18.95 

542.5 

1.5 


520.9 

5-21.0 

. 0 . i • 

521.0’ 

0.1 

521.0 

o;i 

0.03 

0.04 

0.07 

520.9 • 

. 0 


172.3' ■■ 

177.3 

5,0 

177.1 

4.8 

181.9 

’ 9,6 

7 O.I 6 

6.98 

6.82 

170.4 

- 1.9 

( 


All temperatures in 






























































































































TABLE 6.3(g) Model temperatures before and after correction 



Cooling Transient: 30 minutes after spacecraft moving into shade 


Surface 

Prototype 
or Perfect 
Model 

Temperature 

T 

"p 

Temperatures of Imperfect Models 



A 

t 

- . 

Model MI 

Model, Mil 

Model Mill 

Model Temperature 
after Correction 

— 

T 

r/3 

r»* 

hZ 

T 

MI I 


T P 

mil j nvnn 

T 

MC 


^ ■ 

366.5 

368.5 

0 

366.5 

0 

386.4 

-0.1 

0.03 

-0.10 

-0.07 

366.5 

0 

2 

288.6 

289.0 

0.4 

289.6 

1.0 

288.3 

-0.3 

0.48- 

-1-.04 

-0.56 

289.6 

l.'O 


512.5 

539.1 

26.6 

540,8 

28.3 

556.2 

HD 

1.46 

25.5 

27.01 

512.1 

1 

-0.4 ■! 


370.5 

370.4 

-0.1 

3*70.5 

0 

370.3 

-0.2 

0.02 

-0.15 

-0.13 

370.6 

0.1 


390 . 4' 

390.0 

-0.4 

390.2 

-0.2 

389.5 

-0,9 

0.16 

-0.81 

-0.65 

390.7 

• 0,3 


529.5 

557., 7 

28.2 

560.7 

31.2- 

575.4 

45.9 

2.72 

'26 . 51 

,29.23 

528.5 

■ -1.0 


, 536,7 

558,2. 

' 

21.5 

558.1 

21.4 

■571.4 

34.7 

CO 

o 

o . 

19,75 

• — 1 
CO 

CT> 

H 

538'. 5 

1.8 

^ ■ 

529.0 

554.0 

25.0 

554.'5 

25.5 

568.8 

39.8 

0 .-if.O 

22.11 

22.60 

531.4 

2,^ 




552.6 

572.6 

mm 

571,3 

18|. 7 

587.2 

Bi 

-1.19 

21.75 

20.56 

552.1 

-0.5 

^0 

h 

534.3 

555.7 

21.4 

555.8 

21.5 

568.8 

34.5 

0113 

19.66 

19.79 

535.9 

1.6 


383.1 

383.3 

0.2 

383.'3 

0.2 

383.3 

0,2 

0.05 

0.06 

0.11 

383.1 

0 


149.7 

157.1 

7.4 

.157.1 

■ 7.4 

.163,5 

CO 

« 

CO 

o 

o 

9.64 

9.65 


-2,3 


All temperatures in °R. 
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TABLE 6.3(h) Model temperatures before' and after correction 


Cooling Transient: 200 minutes after spacecraft moving into shad4 


■Surface • 

Prototj-pe 
or Perfect 
Model 

Temperature • 

Temperatures of Imperfect Models • 

/ ■ 
' 

\ 

A 

Model Temperature 
after Correction 

■Model MI 

Model Mil 

Model Mill 


E 

Mt 

T 

MI I 


T 

MI I I 

■ E 

: Mil 

1 

T 

MC 

■ 


227.8 

227.9 

0.1 

228.1 

0.3 

227. 7 

' -0.1 

0.16 ' 

-0.33- . 

-0.17 

223.0 

0.2 


269.9 

270.0 

0.1 

270-.7 

0.8 

269.2 

-0.7 

mam 

-1,131 

-0.71 

270.7 

0.8 


500.5 

525.3 

24; 8 

526.4 

25.9 

541.9 

'41,4 

0.97 

■M 

25,77 

499 . 6 

1 

-0.9 ( 


228.8 

228. 8 

0 

229.1 

0.3 

228.5 

■ -0.3 

0.26 

-0.52 

-0.26 




273.4 ■ 

273.5 

0.1 

274.5 

0.1 

272.2 

-0.2 

0.90 

-1.90 

-1.0 

mm 

'1.1 


492,4 

519; 7 

27.3 

521.7 

29.3 

537.2 

44.8 

1.85 

26.26 . 

28.11 

, 

491.6 

'-0.8 


509.6 

530.6 

21.0 

531.4 

21.8 

543.7 

34.1 

0.78 

19.58 

20.36 

. 510.2 

i 

1 

0.6 

• \ 

510.9 

534.4 

23.5 

534.6 

23.7 

548. 8 

37.9- 

' O'. 23 

21,60 

21.83 

512.5 

1.6 


526.2 

545.1 

18.9 

544 ;0 

17.8 

559.6 

33. 


21.60 

20.54 

524.6 

■ -1.6 

^o 

509.0 

529.4 

20.4 

529. 8 

20.8 

542.5 

33.5 

0.30 

19.55 

19.85 

509.6 

0.6 


234.0 

234.4 

0.4 

234.6 

0.6 

234.6 

0.6 

0.19 

...0.22 . 

0.41 

234.0 

0 


BSSIIIIIIII^ 

143.1 

8.3 

_ 143:2 

8.4 

150.2 

15.4 

0.06 

10.67 

10.73 

132.3 

-2.5 

—■■■'? 


All temperatures in 



















































































































































are designated by ^ » and • Equation (2,3.8) is used 

to evaluate A , which is the correction for all positive errors , 
and A^ ^ , the correction for all negative errors , is evaluated from 
(2.3.13). The combined correction A^ is simply the algebraic sum 
of A^^^ and A^^^ . Finally, the correcxed model temperatures calcu- 
lated according to (2.3.14) are denoted by and their errors are 
. In either set of the foregoing tables , data for the heating 
transients are listed in tables (a) to (e), the steady state in ta- 
ble (f), and the cooling transients in tables (g) and (h) . 

An inspection of these tables shows that the proposed correction 
theory is indeed capable of providing good results when errors in 
kd and in Cd of the 'original' model MI are up to 25 percent and er- 
rors in surface emittance up to 20 percent (corresponding errors in 
models Mil and Mill are much larger). As an example, take the case 
of heating transient at t = 100 minutes after the sudden change of 
the thermal condition. Uncorrected model temperatures of errors up 
to approximately 10°R in MI, 26°R in Mil, and 13°R in Mill are brought 
to within 4®R of the prototype value after correction [see TABLE 6.2(c)3 
For the steady state, corresponding errors of up to 21.3, 20.6, and 
34.5®R are reduced to less than 2.5*^R [see TABLE 6.2(f)]. At a time 
30 minutes after cooling begins, uncorrected model temperatures ex- 
'hibit errors of up to 23.8, 24.9, and 39.3°R, respectively, for the 
three models. After correction, the ma^mum error is 2.5®R [see TA- 
BLE 6.2(g)]. Data for the 45° solar incidence show similar results. 
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To test further the performance of the theory, a second set of 
errors in kd, Cd, e, and e’ vras studied. They are listed in TABLE 
6.4. The magnitudes of the various errors are approximately the same 
as before. However, both positive and negative errors now occur in 
kd and Cd. Moreover, the modifications for model Mill are carried 
out in such a manner that the negative errors become Zess negative. 

The results of this study are summarized in TABLES 6.5(a)-(h) when 
the spacecraft receives the sun*s radiation at normal incidence and 
in TABLES 6.6(a)-(h) for the 45° oblique incidence. The spacecraft's 
initial thermal condition, the manner of the initiation of the heat- 
ing and cooling transients, the relative distribution of internal 
povrer dissipation within the bus, etc., are identical to those con- 
sidered in the earlier study. All data again lend support to the 
effectiveness of the proposed correction technique. The errors 
whiph remain in the corrected model temperatures are, in-general, 
smaller than those shorn in TABLES 6.2 and 6.3. Our limited experi- 
ence seems to indicate that the theory would yield better results 
when positive and negative errors coexist in each set of the model- 
ing parameters. 

6.2 SPACECRAFT SUBJECTED TO CYCLIC HEATING AND COOLING 

s 

For an orbiting spacecraft , the temperatures of its various com- 
ponents undergo continuous cyclic changes. Under such condition, 
there was the concern that the error in the model temperatures might 
accumulate and it was not certain "that the proposed correction technique 


96 



TABLE 6,4 Errors in modeling parameters 


Surface 

Percentage error in modeling parameters 

Model MI 

Model MU 

Model Mill 

X M> 

X 100 

6 

X 1^0 

X 1^0 

'^kd 
X 100 

'cd 
X 100 

6 

X 1^0 

X ioo 

/.d ' 

XlOO 

X 100 

6 

X ioo 

< ioo 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 . 

0 

^2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


-20.0 

20.0 

-28.8 

10:0 

-20.0 

42.0 

-28.8 

20.0 

-11.0 

20.0 

-18.0 

10,0 


0 

Q 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

. 0 

0 

0 

0 

0 

0 


18.0 

-15,0 

-32.0 

10.0 

39.6 

-15.0 

-32,0. 

20.0 

18.0 

-8,2 

-20.0 

10.0 


20.0 

15.0 

15.0 

-14.7 

42.0 

33.0 

27.0 

-14.7 

20.0 

15.0 

15,0 

-8.0 


8.0 

-15.0 

-25.9 

-11.5 

19.2 

-15.0 

-25.9 

-11.5 

8.0 

-9.0' 

-18.0 

-8.0 


-10.. 0 

-15.0 

25.0 

-16.0 

-10.0 

-15.0 

45.0 

-16.0 

-6.5 

-8.2 

25.0 

-8.0 


-18.0 

15.0 

15.0 

-12.8 

-18.0 

36.0 

33.0 

-12.8 

-10.8 

15.0 

15.0 

-8.0 

■^11 

0 

0 

0 

imiiQ^i 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

15.0 

0 

0 

0 

33,0 

0 

0 

0 

' 15.0 

0 


(S refers to the exterior of the surface in question and 6 , refers to its interior 
0 S 













































































TABLE 6.5(a) Model temperatures before and after correction 


Heating_ Transient : 10 minutes after spacecraft stibject to solar irradiation (normal incidence) and internal heating 


Surface , 

Protot3rpe 
or Perfect ' 
Model 

Temoerature 

Temperatures of Imperfect Models 




Model Temperature 
after Correction 

Model MI 

Model Mil 

Model Mill 

T 

r/2 

T 

M[ I 


T 

MI X I 

E 

XvE XI 

T 

NxC 

^IvIC 

^ ■ 

518.1 

518.0 

-0.1 

518.0 

-0.1 

518.1 

0 

-0.01 

-0.09 

-0.10 

518.1 

0 


520.5 

519.9 

-0.6 

519.8 

-0.7 

520.1 

-0.4 

-0.04 

-0.54 

-0.58 

520.4 

-0.1 


429.4 

429.2 

-0.2 

425.3 

BB 

427,7 



3.78 • 

. 0.19 

429.0 

-0.4 


518.1 

■518.0 

-0.1 

518.0 

-O.l' 

518.0 


~0 

-0.07 

. -0.07 

518.1 

0 


520.2 

519.8 

-0.4 

519.8 

-0.4 

519.9 

-0.3 

~0 

-0.41 

-0.41 

520.2 

0 


412.2 

416.7 

4.5 

416.7 

■jU 

414.7 

2.5 

-0^01 

BB 

Bsa 

411.6 

. -0.6 


. 401.5 

401.5 

0 

401.3 

-0.2 

401.4 

-0.1 

-0.16 

0.15 

-0 .01 

401.5 

0 

• ^ 

420.7 

427.8 

mam 

427.6 

6.9 

425.0 

4.3 

-0.19 

7.09^ 

6.90 

■ 420.9 

0.2 


415.9 

418.3 

2.4 

417.8 



■a 

-0.43 

3.55 

3.12 

415.2 

-0.7 


403.0 

403.0 

• 0 

402.5 

1 -0.5 

402.9 

-0.1 

-0.44 

0.40‘ 

-0.04 

403.1 

0.1 


476.3 

476.3 

0 

478.3 

• 0 

476.3 

0 

~0 

~0 

0 

476.3 

0 


152. 8 

153,4 

0.6 

153.4; 

0.6 

153.1 

0.3 

-0.05 

0.90 

0.85 

152.6 

■ -0.2 


All temperatures in °R 



































































































































TABLE 6.5(b) Model temperatures* before and after correction 


Heating Transient; 30 minutes after spacecraft subject to solar. Irradiat.ion Cnormal ?)-nd internal heatins 


Surface 

Protot 3 ?pe 
or Perfect 
Model 

Temperature 

T, _ 



Temperatures of Imperfect Models 

' — — 

. 

A- 

Model Temperature 
after Correction.. 

/Model MI- 

Model Mil 

Model Mill 



T 


T 

MI 1 1 

E 

IvOI 

T 

.ViC 


A, 

590.7 



590.5 

-0.2 

590.6 

-0.1 

• -0.03 

-0.12 

-0.15 

590.7 

0 

A 

''■2 

594.9 

594.0 

-0.9 

593., 8 

-1.1 

594.3 

-0.6 

-0,17 

. 

-0*70- 

-0.87 

594.9 

0 . 


474.5 


3.5 

470.2 


473.9 

-O'. 6 

-7.19 

10.52 

3:33 

474.7 

0.2 


590.5 

590.5 

0 

590.5 

0 

590.5 

0 

'~0 

-0.04 

-0.04 

590.5 

0 


593.9 ■ 

593.6 

-0.3 

593.6 

-0.3, 

593.7 

- 0.2 

~0 

-0.24 

-0.24 

' 

593.9 

0 

h 

440.3 

451 . 6 

11.3, 

'451.4 

11.1 

446,4 

6.1 

-0.24 

13.27 

13.03 

438.6 

■-1.7 

. ^ 

. 413.1 

413.4 

0.3 

411.8 

mm 

412.8 

-0.3 

-1.49 

1.60 

0.11 

413.3 

0.2 

! 


450.6 

' 

466.2 

15.6 

46 5. a 

14.5 

459.9 

9.3 

-0..99 


14.89 

451.3 

0.7 1 

K 

445.4 

' 



■450 ..6 

5.2 

448.4 

3.0 

-1.73 

10.25 

8.52 

443.9 

-1.5 1 


424.3 



422.8 

-1.5 

424.3 

0 

-3,00 

4.43 

1.43 

424.6 

0.3 1 

A I 

547.6 

547.6 

0 

547.6 

0 ' 

547.6 

0 

0.01 

-0 

0.01 

547.6 

0 } 

7, *'"^1 0 ' ^ 

172.2 

172.8 

0.6 

172.6 

0.4 

172.5 

. 0.3 


0.98 

. -0.75 

172.1 

'- 0.1 1 
. [ 


All temperatures in °R 
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TABLE 6.5(c) Model temperatures before and after correction 


Heating Transient: 100 minutes after spacecraft subject to solar irradiation (normal incidence) and internal heating 


j : 

Surface 

Prototype 
or Perfect 
Model 

Temperature 

Temperatures of Imperfect Models 



4 .. 

Model Temperature 
after Correction 

Model MI 

Model Mil 

Model Mill 

T 

17 

T 

-Mil 


T 

MI 1 1 

^IvH 1 1 

T 

MC 

^r/ic 

A, • 

. 597.8 

597.7 

-0.1 

597.6 

-0.2 

597.7 

-0.1 

■ -0.03 

-0.10 

-0.13 

597.8 

0 


603.8 

603.1 

-0.7 

602.9 

-0.9 

603.4 

-0.4 

-0.19 

-0.60 

-0.79 

603.9 . 

0.1 


519.1 

542.0 

22,9 

537.4 

18.3 

531.3 

12.2 

-4.26 

27.24 

22„98 

519.0 

-0.1 ■ 


597.6 

597.6 

0 

597.6 

0 

597.6. 

0 

-0.01 

-0.04 

-0.05 

597.6 

0 


602.8' 

602.5 

-0.3 

602.5 

-0.3 

602.6 

-0.2 

-0.07 

• 


-0.30 

602.8 



498. 8 

. 523.7 

24.9 

521.6 

22.8 

512.4 

13.6 

-1.96 

'28.72 

26.76 

496.9 

-1.9 


■ 479.6 

490.1 

10.5 

482. 5 _ 

28.9 

484.0 

4.4 

-7.00 

15.38 

• 8.38 

481.7 

2.1 


499.5 

526.1 

26.6 

523.5 

24.0 

515.4 

15.9’ 

-2--.45^ 

27.20, 

24.75 

501.4 

1.9 


511.6 

530,6 

20.0 

525.8 

14.2 

520,9 

9.3 

-4.41 

24.56 

20.15 

510.5 

-1.1 


490.8 

507.6 

16.8 

. 

501.9 

11.1 

499.3 

8.5 

-5.21 

20.97 

15.76 

491.8 

1.0 


568.0 

568.1 

0.1 

568.1 

0.1 

568.1 

0.1 

0.01 

0.02 

0.03 

568.0 

0 


. 181.1 

183.7 

2.6 

183.1 

• 2.0 

182.4 

1.3 

-0.56 

3.46 , 

2.90 

180.8 

'-0.3 


All temperatures in 
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TABLE 6.5(d) Model temperatures before and after .correction 
Heating Transient: 300 minutes after spacecraft subject to solar irradiation (noiroal incidence) and internal heating 


Surface 

' 1 

Prototype 
or ; Perfect 
Model 

Temperature 
"i Tp 

Temperatures of Imperfect Models , 

' 


5 

1 . 

At ! ■ 

1 

• 

Model MI 

Model Mil 

Model Mill 

Model Temperature 
after Correction 


E ' 

T 

.MI I 

•p 

• T 

MI 1 1 

; M XI 

T 

MC 

V ’ 

"I/C 

A 

598.0 

597.9 

-0.1 

597.9 

-0.1 

597.9 

‘ "0,1 

-0.01 

-0.07 ■ 

-0.08 

5.9..8..-Q— 

1 


604.7 

604.3 

-0.4 

604.3 

-0.4 

604.5 

' -0.2 

•-0.06 

-0;43 

-0.49 

604.8 

0.1 1 


533.9 

565.9 

32.0 

564.6 

30.7 

552.2 

18.3 

-1.13 

34.62 

33. .49 

532.4 

-1.5 ,| 

A 

597.9 

5'97.8 

-0.1 

597.8 

.-6.1 

597.9 

0 

-0.01 

-0.04 

.-0.04 

597.9 

0 


604. 1 

604.0 

-0.2 

604.0 

-0.2 


' 0 

-0.05 

-0.21 

n 

604.3 

0.1 


524.4 

558.4 

34.0 

557.3 

32.9 

543.9 

■ 19.5 

-1.03 

36.66 

■35.63 

522.8 

-1.6 j 

• \ 

527.7 

556.5 

28.8 

554.2 

26.5 

544.7 

i 17.0 

-2.17 

29.87 

27.70 

528.8 

w 


IIIIIIIIIQQIIII^ 



557.2 

33.2 

545.2 

- 21.2 


33.60- 

32.41 

526.1 




574.9 

28.7 

571.8 

25.6 

562.2 

16.0 

— 

32.09 

29.23 

545.7 

-0.5 

o 

< 

527.7 



554.4 

26.7 

544.8 

^ 17.1 . 

-1.98 

■■ 30.00 

28.02 

528.6 

0.9 


568.1 

568.2 

0.1 

568.2 

0.1 

568.2 


0.02 

0.04 

G.06 

568.1 

0 


184.9 

189.5 

■B 

189.3 

4. 4 

187.4 

2.5 

-0 . 26 

5.53 

5.27 

184.3 

-0.6 


All temperatures in °R. 
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TABLE 6.5(e) Model 


Heatins Tr’an‘?T<=>Tvt: SOO minutes after spacecraft subject 


Surface 

Prototype 
or Perfect - 
Model 

Temperature 

T 

Temperatures of Imperfect 

■•Model MI 

Model Mil 

#r* 

E^T 

tilt 

T 

MI. I 


■ 

598.0 

597.9 

-0.1 

597.9 

-0.1 


604.8 

604,4 

-0.4 

604.3 

-0.5 


534.4 

566.9 

32.5 

566.2 

31.8 


597.9 

597.9 

0 

597.8 

-0.1 


604.3' 

604.2 

-0.1 

'604.1 

-0.2 


525.3 

559.9 

34.6 

559.2 

33.9 , 

'•V 

529,4 

559.2 

29.8 

557 ; 9 

28.5 ; 


525.0 

560.0 

35.0 

559.1 

34.1 


5.47.5 

■ 

576.8 

29.4 

574.4 

26.9 ■ 

^0 

529.0 

558.7 

29.7 

1^1 

28.3 


568.1 


0.1 

568.9 ■ 

0.1 


185.0 
■ ■ 

189.8 

4.8 

189.6 

mn 


All temperatures in °R 


before and after correction 


to solar irradiation (normal incidence) and internal heating 


Models 

. 

a'->- 

A . 

t , ' 

Model Temperature 
after Correction 

Model Mill 

• T 

hn 1 1 

1 1 

T 

MC 


597.9 

-0.1 

-0.01 

-0.07 

-0.08 

598.0 

0 ■ 1 

604.6 

-0.2 

-0.03 

^0.42 

-0.45 

604-9 

0.1 j 

553,2 

18.8 

-0.72 

34.91 

34.19 

533.7 

-0.7 .j 

597.9 

0 

~o 

-0.04 

-0.-04 

597.9 1 0 . j 

604.2 

-0.1 

,-0.03 

-0.21 

-0.24 

1 

604,4 

• 0.1 

545.3 

20.0 

-0.60 ■ 

36.95 



• -0.5 

547.3 

17.9 

-1.20 

30.31 

29.11 

531.5 

2.1 

546.6 

,21.6. 

-0.75 

33.83 

33.08 


2.7 

564.1 

16.6 

-2.28 

32,34 

30.06 

547.1 

-0.4 

546.7 

17.7 

-1.28 

• 30,29 

29.01 

530.9 

1.9 

568.2 

0.1 



Op 06 

568.1 

0 


187.6 


2.6 


-0.18 


5.44 


183.7 
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Steady State 


TABLE 6.5(f) Model temperatures before and after correction 


Surface 

Prototy^je 
or Perfect 
Model 

Temperature 

• 

Temperatures of Imperfect Models 

■ . 

' 



Model Temperature 
after Correction 

Model MI 

Model Mil 

Model Mill 

T 

T? 

T 

MI I 

I. 

T 

f/t 1 1 

E ■ 

M r I 

T 

.ViG 

1 

A, 

597.9 

597.9 

0 

597.9 

0 

597.9 

0 

-0.01 

-0.07 

-0.08 

598.0 

0.1 


604.8 

604.4 

-0.4 

604.4 

-0.4 

604.6 

-0.2 

-0.03 

-0.42 

-0.45 

‘ 

604.8 

0 


534.4 

567.0 

32.6 

566.2 

31.8 

. 

553.2 

18.8 

-0.69 

34.91 

34.22 

532.8 

-1.6 , 


597.9 

597.9 

0 

597.8 

-0.1 

597.9 

0 

. 

-0 

,-0.04 

-0.04 

597.9 

i 

0 

-s 

604.3 ■ 

604.2 

-0.1 

604.1 

-0.2 

604.2 

-0.1 

•-0.03 

-0,21 

-0.25 

604.4 

0.1 


525.3 

559.9 

34.6 

559.3 

o 

CO 

545.3 

20.0 

-0.57 


36.39 

523.5 

-1.8 


529.5 

559.4 

29.9 

558.1 

28.6 

547.4 

17.9 

-1.13 

30.32 

29.19 

530.2 

0.7 


525.0 

560.0 

35.0 

559.2 

34.2 

546.7 

21.7 

1 

o 

33.83 

33.12 

526.9 

1*9 

miiQii 

. 547.5 

576.9 

29.4 

574.5 

27.0 

564.2 

16.7 

-2.23 

32.34 

30.11 

546. 8 

-0.7 


529.0 

558.8 

•29.8 

557,5 

28.5 

546.8 

17.8 

-1.22 

30.30 

29.08 

529,7 

0.7 


568.1 

568.2 

0.1 

568.2 

0.1 

568.2 

0.1 

0,02 

0.04 

0.06 

568,1 

0 


185.0 

, 

189.8 

4.8 

189.6 

mn 

187.6 

2.6 

-0.17 

5.62 

•5.45 

184.4 

■-0.6 
. ■ 1 


All temperatures in 
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TABLE 6.5(g) Model temperatures before and after correction 


Cooling Transient: 30 minutes after spacecraft moving into shade 


Surface 

Prototype 
or Perfect 
Model 

Temperature 

Tp 

Temperatures of Imperfect Models 

— 

. 

4 ■ 

Model Temperature i 
after Correction | 

Model MI 

Model Mil 

Model Mill 

T 

W 

MI 

T 

MI I 


T 

MI 1 1 

1 1 

m 

■^MC 



376.5 

376.4 

-0.1 

376.4 

-0.1 

376.5 

0 

~o 

-0.19 ■ 

-0.19 

376.6 

0.1 

A 

390.5 

389.7 

-0.8 

389:7 

-0.8 

390.1 

-0,4 

~o 

-1.04 

' 

-1.04 

390.7 

0.2 


516.2 

554.4 

38.2 

554.7 

38.5 

539.0 

22/8 

0.33 

38. 89 

39.21 

515.1 

- 0.9 'i 

^ . 



H 

376,2 

0 

376.3 

- 0.1 

-0.01 

-0.23 

-0.24 

376.5 

. t 

0.1 

1 ^ 

389.9 

388.8 

- 0.1 

388.8 

-0.1 

389.3 

- 0.6 

-0.05 

-1.29 

-1.34 

390.2 

0.3 


• 

513,0 

550.4 

37.4 

549.9 

36.9' 

534.3 

21.8 

-0*43 

39.44 

39,01' 

. 511.4 

■■ 

. ^ 

527.0 

557. if 

30.4 

556.3 

29,3 

545.3 

18.3 

- 0.95 

30.61 

29.66 

527.7 

0.7 


522.6 

557.9 

35.3 

557 .‘2 

34.6 

. 

544.5 

21.9 

-o'. 6 9 

34.11' 

33.42 

524.5 

1.9 


543.4 

573.6 

30.2 

571 .'3 

27.9 

560.6 

17.2 

-2.14 

32.92 

30'. 78 

542.8 

-0,6 

'So 

525.5 

556.0 

30.. 5 

554.9 

1 

29.4 

543.8 



30.74 ■ 

29.72 

526.3 

0.8 


395.2 

395.3 

0.1 

MHi 

0.2 

395.3 

0,1 

-0.04 

-0.07 

-0.11 . 

395.2 

0 


150,4 

159.2 

8.8 

158.9 

8.5 

155.3 



9.96 ■ 

9.72 

149.5 

-0.9 


All temperatures in °R. 
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TABLE 6.5(h) Model temperatures before and after correction 


Cooling Transient: 200 minutes after spacecraft moving’ into shade 




Surface! • 

Prototype 
or Perfect 
Model 

Temperature^ 

Tp 

Temperatures of Imperfect Models 


• 

4 .. 

1 

Model Temperature! 
• after Correction 

- ■ Model MI 

Model MU 

^?odel Mill 

m 

T« 

T 

MI I 


T 

1 1 

E 

M 1 1 

T 

MC 


^ ■ 

229.8 

229.4 

-0.4 

229.4 

-0.4 

229.6 

-0.2 

. I I 

' -0.01 

-0,57 ■ 

-0.58 

230.0 

0.2 

'V- 

275.7 

274.2 

-1.5 

274.1 

-1.6 

275.1 

-0.6 

-0,08 

CO 

o 

CN 

1 

-2.16 

276.4 

0.7. 


500.3 

540.3 

‘40.0 

539.8 

39.5 

523.7 

23.4 

-0.39 

41.84 

’41.45 

498.8 ‘ 

1.5 . 


229.0 

228.4 

-0.6 

228.4 

-0.6 

228.7 

-0.3 

’ -0.01 

'-0.92 

■ -0.93 

229.3 

0.3 


272.7 

270.3 

B 

270.3 

B 

271.7 

-1.0 

-0.05 

-3.38 

-3.43 ' 

273,8 

•1.1 


491.8 

534.5 

42.7 

534.3 

42.5 

517.0 

25.2 

-0.16 

44.34 

44.18 

490.3 

-1,5 


508. 8 

542. 6’ 

33.8 


32.6 

523.0 

20.2 

-1.01 

34.30 

'33,29 

509.3 

0.5 


510.4 

547.8 

37.4 

547 ;o 

36.5 


■23.0 

-Cf.81 

36.53 

35,72 

512.1 

1.7 ■ 

‘^0 

525.6 

559.6 

34.0 

557.2 

31.6 

545.0 

19.4 

-2.15 

36.92 

34.77 

524.8 

-0.8 


508.3 

542.0 

0- 

CO 

CO 

540.7 

32.4 

528.5 

20.2 

-1.20 

34.40 

33.20 

508. 8 

0.5 


235.6 

236.1 

0.5 

236.3 

0.7 

236.0 

0.4 

0.19 

0,34 

0.53 

235.6 

0 


134.8 


11.5 

145.9 

11,1 

141.2 

6.4 

-0.33 

12.90 

. 12.57 

133.7 

-1.1 


All temperatures in °R. 
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TABLE 6.6(a) Model temperatures, -.before and after correction 


Heating Transient: 10 minutes after spacecraft subject to solar irradiation (45® incidence) and internal heating 


Surface 

Prototype 
or Perfect 
Model 

Temperature 

Temperatures of Imperfect Models 



, 

■ Model Temperature 
after Correction 

-• Model MI 

Model Mil 

Model Mill 

T 

r.'l 

T 

T 

MI 1 

r 

T 

Mt 1 1 

E ■ 

M ri 

T 

MC 



477,4 

477.3 

-0.1 

477.3 

-0.1 

477.3 

-0.1 

-0.01 

-0.12 

-0.13 

477.4 

0 


370.1 

369.3 

-0.8 

369.2 

-0.9 

369.6 

-0.5 

-0.05 

-0.74 . 

-0.79 

370.0 

-0.1 

’a 

425.7 

427.0 

1.3 

423.4 

-2.3 . 

425.1 

-0.6 

-3.30 

4.63 • 

1.33 

425.6 

-0.1 


' 481.0 

480.9 

-0.1 

480.9 

-0.1 

481.0 

0 

■~o 

-0.09 

-0.09 

481.0 

0 


499.3 

498.7 

-0.6 

498.7 

-0.6 

498.8 

-0.5 

~0 

-0.50 

-0.50 

499.2 

-•<•0.1 


430.4 

438.1 

7.7 

438.0 

7.6 

434.5 

4.1 

-0.15 

• 9.22 

9.07 

429.1 

-1.3 


401.7 

401,8 

0.1 

401.6 

-0.1 

401.7 

0 

-0.16 

0.23 

0.07 

401.7 

0 


421.3 

428.7 

7,4 

428.6 

7.3 

425 . 8 

4.5 

-0.09 

7.45 

7.36 

' 421.3 

0 


416.9 • 

419.8 

2.9 

419.4 

■ 2.5 

418.2 

1.3 

-0.36 

4.07 

3.71 

415.1 

-0.8 

■ 

403.6 

403.8 

0.2 

403.2 

-'0.4 

403.5 

-0.1 

-0..48 ■ 

0.61 

0,13 

403.6 

0 


447.6 

447.6 

0 ^ 

447.6 

0 

447,6 

0 

0.01 

~o 

0.01 

447.6 

0 

> 

n 

o 

142.1 

143.0 

0.9 

143.0 

0.9 

142.5 

0.4 

-0.06 - 

1.25 

1.19 

141.8 

-0.3 


All temperatures in °R. 




TABLE 6.6(b) Model temperatures before and after correction 

Heating Transient: 10 minutes after spacecraft subject* to solar irradiation (45° incidence) and internal heating 


Surface 

' 

Prototype 
or Perfect 
Model 

Temperature 

T. 

Temperatures of Imperfect Models 


, 

* 1 

. ‘ 

} 

Model Temperature 
after Correction 

•Model MI 

Model Mil 

Model Hill 

T 

T* 


^ivn I 

T 

M I ! 

X I 

T 

MC ^ 




535.4 

-0.2 

535.3 

-0.3 


-0,2 

-0,03 

- 0.21 

-0,24 

. 535.6 

0 


335.3 

333.0 

-2.3 

332.7 

-2.6 

333.8 

-1.5 

-0.30 

-1.90 

■ -2,20 

335.2 

- 0.1 


457.7 

467.0 

9.3 

460.4 

B 

461.2 

3.5 

' -6,06 ■ 

14.68 

8.62 

458.4 

0.7 >! 


541.9 

541.8 

-0.1 

541.8 

-0.1 

541.9 

0 

~0 

-0.10 

-O'. 10 

541.9 

0 


568.3 


-0.6 

567.7 

-0.6 

567.9 

-0.4 

-0.02 

-0.54 

-0.56 

568.3 

0 


480.6 

499,2 

18.6 

498.1 

17.5' 

490.4 

9.8 

-0.98 

22.26 

,21,28 

477.9 

- 2.7 

■''7 

414.9 ■ 

416.0 

H 

414.4 

-0.5 , 

415.1 

0.2 

-1.55 

2,46 

,0.91 

415.1 

0.2 


453.5 

471,1 

17.6 

470.5 

17.0 

464.1 

10.6 

■ -0.54 

17,94 

17.40 

453.7 

0.2 


450.5 

460.4 

9.9 

458.9 

8.4 

455.2 

IQI 

t1.42 

13,25 

11.83 

448.6 

- 1.9 


427.6 

431.2 

3,6 

427.8 

' 0.2 

428.7 

1.1 

-3.11 

. 6.35 

3.24 

427.9 

0.3 

^ix 

498.0 


0 

498.0 

0 

498.0 

0 

0.01 

0.01 

p.02 

498.0 

0 


155,5 

156.7 

■B 

156.3 

0.8 

156.0 

0.5 

-0.33 

1.76 

1.43 

155.2 

- 0.3 


All temperature in °R. ' 
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TABLE 6.6(c) Model temperatures before and after correction 


Heating Transient: 100 minutes after spacecraft siibject to solar irradiation (45° incidence) and internal heating 


Surface 

Prototype 
or Perfect 
Model 

Temperature 
. . Tp .... 

Temperatures of Imperfect Models 


. t 


’ ' * 1 

Model Temperature , j 
after Correction j 

' Model MI 

_ Model Mil 

Model Mill 

T 

‘1/Z 


T 

Ml I 


T 

m 1 1 

E 

!vt 1 1 

T 

.ViC 

' i 

< 

^ ■ 

546.3 

546.1 


546.1 

-0.2 

546.2 

-0.1 

'-0.04 

-0.17 . 

HU 

546.0 

0 


296.2 

292.4 

-3.8 

291. ’3 

-4.9 

293.7 

-2.5 

-1,05 

-3:i9 

-4.24 

296.7 

0.5 


498.1 

529.3 

31.2 

525,0 

26.9 

515.3 

17.2 

-3.88 

35.30 

31.4:^ 

497.9 

-0.2 i 

1 

1 

I 

553.4 

553.3 

-0.1 

553.2 

-0.2 

553.3 

-0.1 

-»0 .03 

-0.11 

-0.14 . 

553.4 

■1 

.0 1 


581.6 

580.9 

-0,7 

580.7 

-0.9 

581.1 

-0.5 

-0,.18 

-0 . 55. 

-0.73 

581.7 

0.1 


552.8 

. 

585.1 

33.3 

582.5 

29.7 

570.8 

18.0 

-3.31 

CO 

CO 

CO 

CO 

.35.57 

550.6 

, 

-2.2 

• ^ 

490.6 

507.5 

16.9 

499.6 

9.0 



mm 

22.17 

14,88 

492.6 

1 

2,0 


507.9 

539.8 

31.9 

537.6 

29.7 

527.0 

19.1 

-2’.:04.. 

32.56 

30.52 

509.3 

1.4 


524.9 

550.8 

25.9 

546.1 

21.2 

538. 3 

13.4 

•74.27 

31.66 

27.39 

523.4 

-1.5 


502.6 

525.1 

23.5 

520.6 

18.0 

515.0 

12.4 

-5.05 

. 28.01 

22.96 

503.1 

0.5 1 


520.6 

520.7 

,0.1 ■ 

520.7 

B 

520.7 

0.1 • 

0.02 

0.04 

0.06 

520.6 

1 

0 

\ 2 ^ 0 ' ^ 

167.2 

171,8 

B 

171.0 

■ 

3.8 

169.4 

2.2 

-0.74 

5.95 

5.21 

. 

166.6 

-0.6 j 


All temperatures in °R. 
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TABLE 6.6(d) Model temperatures before and after correction 
Heating Transient: 300 minutes after spacecraft siibject to solar irradiation (45® incidence) and internal heating 


Surface • 

1 

1 

, 

• Prototype 
or Perfect 
Model 

Temperature 

T. 

Temperatures of Imperfect Models 

! 

! 

■ 1 

1 

1 

( 


Model Temperature 
after Correction 

•Model MI 

Model Mil 

— 1 

Model •Mill 

T 

r 

T 

wi I 

1 , 

T 

M 1 1 

E 

ivn r I 

T 

wc 


A ' 

546.5 

546,4 

-0.1 

546 . 4 

-0.1 

546.4 

-0.1 

-o.Oi 

-0.11 

-0,12 

546.5 

0 



289.4 

-1.4 

289.0 

-1.8 

290.2 

-0.6 

-0.37 

-1.79 

-^2.16' 

291.6 

0.8 

. 

^ ■ 

515.2 

556.4 

41.2 

556.1 


mm 

.24.-2 

-0.31 

43.10 


513.7 

-1.5 


553. 8 

553.7 

-0.1 

553.6 

. -0-2 

553.7 

-0.1. 

~p.03 

-0.10- 

1 

-0.-13 

553.8 

0 ! 


583.4 ■ 

582.9 

-0.5 

582. 8 


583.1 

mm 

-0 . 13 

-0.50 

' -0.63 

583.5 

0.1 1 


573.4 \ 
' 

614.0 

40.6 

611.8 

38.4 

596.4 

23.0 

-2.p4 

44.73 

42.69 

571.3 

-2.1 1 


540.1 

575.3 . 


573.5 

33*4 

561.0 

20.9 

-1.66 

36.08 

34.42 

540.9 

i 

0.8 1 

■S 

533.1 

572.6 

39.5 

571.9 

38.8 

557.4 

.24; 3. 

-0.69 

38.55 

37.86 

534.8 

MOK 

K 

559.6 

594.5 

34.9 

591.8 

32.2 

579.4 

19.8 

-2.50 

38.20 

35.70 

558.8 

-0.8 


540,0 

575.2 

35.2 


33.5 

560.9 

20.9 

-1.54 

36.12 

34.58 

540.6 

0.6 


520.9 



521.1 

Wm 

521.0 

0.1 

0.03 

0.07 

q.io 

■520,9 

0, 


•• 172.2 
. 

179.4 

7.2 

179.1 

■ 6.9 

176.1 

3.9 

-0.28 

8.44 

8.16 

171.3 

-0.9 

! 


All temperatures in ®R. 
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TABLE 6.6(e) Model temperatures before and after correction 


Keating Transient: 500 minutes after spacecraft subject to solar irradiation (45® incidence) and internal heating 


Surface 

1 

1 

Prototype . 
or Perfect 
Model 

Temperature 

-fp 

Temperatures of Imperfect Models ’ 

! 

t 

’ a'-' 

( 

■ A - •• . - 

Model Temperature.! 
after Correction 

Model MI 

Model Mil 

Model Mill 

m 

r 

T 

. Vi ! 

T 

.VJ I 


T 

MI I r 

1 1 

T 

.MC 

■ P 

f 

! ^ 1 546.5 

546.4 

“ 0.1 

546.4 

- 0.1 

546.4 

- 0.1 

~o 

- 0.10 

- 0.10 

546.5 

0 

1 291.4 

290.6 

- 0.8 

290 ..6 

- 0.8 

291 , 2 

- 0.2 

~0 

- 1.57 

- 1.57 

292.1 

0.7 


515.7 

557 . 3 

41.6 

557.4 

41.7 

540.2 

24,5 

• 0.09 

43.23 

43.32 

514.0 

- 1.7 , 


553.8 

553.7 

- 0.1 

553.7 

. - 0.1 

553.7 

- 0.1 

- 0.02 

- 0.10 

- 0;12 

■ 553.8 

0 




- 0.4 

582.8 

- 0.6 

583.2 

- 0.2 

. - 0 .11 

■ - 0.51 

- O . C '2 . 

583.6 

0,2 

■Bi 

573.9 

614.8 

40.9 

613.0 

39.1 

597.1 

23.2 



43.09 

571.7 

■ - 2.2 

S ' 

‘ 541.3 

577.2 

35.9 

576.2 

34.9 

562.9 

21.6 

- 0.87 

36.18 

35,31 

541.9 

0.6 1 
t 

• 

533.7 

573.7 

• 40.0 

573.3 

39.6 

558.4 

24.7 

- C .32 

38 . 60 . 

• 38.28 

iQimii 

bHI 

^ . 

, 560.5 

595.9 

35.4 

593..7 

33.2 

580.8 

20.3 

- 2. ‘02 

38.23 

36.21 

559.7 

- 0.8 

1 0 

541.0 

576.6 

; 35.6 

575.6 

34.6 

562.4 

21.4 

• 

- 0.97 

36,16 

‘ 35.19 

541.4 

0.4 

A., 

a . 



0,1 ■ 

521.1 

0.2 

521.0 

0.1 

0.04 

0.07 

0.11 • 

520.9 

0 


172.3 

179.7 

7.4 

179 .'5 

■H 

176.3 

4.0 

.- 0.18 

8.48 

8.30 

171.4 

jjiimmu 


All temperatures in ®R. 
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TABLE 6.6(f) Model temperatures before and after correction 


steady State 


Surface 

Prototype 
or Perfect 
Model 

Temperatures of Imperfect Models 

« 

. t 

A ■..- 

Model Temperature 7 
after Correction 

Model MI 

Model Mil 

Model Mill 

Temperature 

i • 

m 

~uz 


T 

•VJ I 


T 

hn 1 1 

n 

IQf 

i 

1 ^ 

■ 546^5 

546.4 

-0.1 

546.4 

-0.1 

546.4 

-0.1 

0.01 

- 0.10 

-0.09 

546.5 

0 

( 

291.5 

290.7 

-0.8 

290.7 

-0.8 

291.3 

-0.2 

0,04 

-1.56 

-1.52 

292.2 

0.-7 


515.8 

557.3 

41.5 

557.5 

41.7 

540.3 

24.5 

0.11 - . 

43.23 

4-3.34 



A, 

553.8 

553.7 

rO.l 

553.7 

-0.1 

553.7 

-O.-l 

-0.02 

.-0.10 

-0,12 

■ 1 

553.8 

0 1 

i 

t 

583'. 4 

583.0 

-0.4 

582.8 

-0.6 

583.2 

-0.2- 

-0.11 

‘-0.51 

-0,62 

583.6 

0.2 j 


573.9 

' 

614 ..9 

41.0 ■ 

613.0 

39.1 

597.2 

23.3 

BB 

44.82 

43.11 

571.8 

- 2.1 

!■ V 

541.4 

577.2-, 

35.8 

576.3 

34.9 

. 

562.9 

21.5 

mmm 

36.18 

35.35 

541.9 

0.5 j 


533.8 

573.7 

3-9.9 

573.4 

39.6 

553.5 ■ 

24.7. 

-0*30 

38.60. 

38.30 

535.4 

1,6 



595.9 . 

35.3 

593.7 

33.1 

580.8 

20.2 

-1,9'9 

38.22 

36.23 

559.7 

- 0.9 

■^o 

541.0 

576.7 

35.7 

575.7 

34 . 7 ' 

562.4 

21.4 

-0.93 

36.16 

35.23 

541.5 

0.5 


520.9 

521.0 

0.1 

521.1 

0.2 

521.0 

0.1 

0.04 

0.07 

0.11 

520.9 

0 

'■■> = '=■''10 ■> 

172.3 

■ 

179.7 

mmi 

179.5 

mm 

176.4 

4.1 

-0.17 

— 

rnmmmi 

8.31 

171.4 

- 0.9 

r 

1 


All temperatures in °R. 
























































































































































X12 


TABLE 6.6(g) Model temperatures before and after correction 


Cooling Transient: 30 minutes after spacecraft moving -into shade 




Temperatures of Imperfect 

Models 

• 

, 

, 



. 

Surface 

• Prototype 
or Perfect 

Model MI 

Model 

Mil 

Model Mill 




Model Tei 
after Co: 

npe nature 
erection 

Model 

Temperature 

T„ 

T 


T 

,MI I 

!■ 

T 

Ml 1 1 

E 

ZvI 1 1 



rn 

1 


' ■ J7 1 ... 1— 

366.5 

366.4 

-0.1 

366.4 

-0.1 

366.4 

-0.1 

'~0 

-0.16- ■. 

-0.16 

366.5 

0 


288.6 . 

287.7 

-0.9 

287.7 

-0.9 

^88.3 

-0.3 

0.05 

-1.68 

-1.63 

289.3 

0.7 


512 ..5 

554.2 

41.7 

554.4 

41.9 

537.2 

24. -7 

0.25 

43.03 

43.28 

510.9 

-1.6 '1 

1 ^ 

370.5 

370.2 

- 0 , 3 ' 

370.2 

, -0,3 

370.3 

-0.2 

-0.03 

-0.30 

-0,33 

370.6 

0.1 


390.4' 

389.0 


388.8 


389.6 

-0.8 

-0.16 

-1.56 

1.72 

390.7 

■ 0.3 


529.5 

568,7 

39,2 

567.6 

38. 1 

552.3 

22,8 

-1.02 . 

41.62 

40.60 

528.1 

'-1.4 

. ^ 

'536.7 

571,8 

35.1 

571.2 

34.5 

.558.0 

21.3 

-0.57 

35.16 

34.59 

537.2 

0.5 


529.0 

567.4 

38.4 

566.8 

37.8 

552.8 

23,8 

-0.56 

37.09- 

36.53 

530.9 


% 

5.52.6 

1 ^ 

586,4 

33.8 

584.1 

31.5 

572.0 

19,4 

-2 . 12 

36.51 

34.39 

552.1 

- 0.5 1 


534.4 

568.7 

.34'.4 

567.9 

33.5 

555.1 

20.7 

-0.75 

. 34.54 

mamm 

534.9 

0,5 


383.1 

383.3 ' 

. 0.2 ■ 

383.4 

0.3 

383.3 

0.2 

0.06 

0,11 

0.17 

383.1 

0 

•\ 2 ^ 0 ' ^ 

: 149.7 

1 . 

160.0 

10.3 

159.8 

■ 10.1 

155.4 


-0.19 ■ 

11.48 

J 

11.29 

148.7 

- 1.0 


Ail temperatures in °R 
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TABLE 6.6(h) Model temperatures before and after correction 




Cooling Transient: 200 minutes after spacecraft moving. into shade 


j Surface 

! 

i 

> 

Prototype 
or Perfect 
Model 

Temperature 

T 

p 

Temperatures of Imperfect Models 

''a'*’. 

' A<->- 

, 1 


Model Temperature ' 
after Correction 

r 

Model MI 

Model Mil 

Model Mill 

m 

‘uz 


T 

I 

I- 

T 

r/t 1 1 

r t 

m 

■^MC 

T? 


227.8 

227.4 

-0.4 

227.4 

-0.4 

227.7 

-0.1 

■ ~,o 

-0.57 

-0.57 

22.80 

0.2 

A 

^'2 

269.9 ■ 

268.4 

-1.5 

268.3 

-1.6 

269.3 

-0.6 

-0.03 

-2.27 

-2.30 

270.7 

• 0.8 

j A, 1 500.5 

t ^ 1 

540.7 

40.2 

540.3 

39.8 

524.1 

23.-6 

-0.29. 

42.00 

41.71 



1 

228. 8 

228.1 

-0.7 

228.1 

-0.7 


-0.3 

‘ -0.01 

vO.94 

■ -0;95 

229.1 

0.3 1 


273.4 


B 

271.0 

-2.4 

272.4 

-1.0 

•-0.03 

• -3.36 

-3.39 

274.4 

1.0 


492.4 

535.2 

42.8 

mma 

42 . 7- 

517.6 

25.2 


44.46 

,44.38 

490.8 

-1.6 


509.6 

543,7 

34.1 

542.8 

33.2 

530.1 

20.5 

-0.82 

34.50 

33.68 

510.0 

0.4 

. 1 


510.9 

548.5 

37.6 

547.7 

36.8 

534.0 

23;1 

■-0,73 

IIIIIQ 

35.92 

512.5 

—1 

", 

526.2 

560.11 

34.2 

558.2 

32.0 

•545.8 

19.6 

.-2.04 

37.04 

35.00 

525.4 

-0.8 


509.0 

542.9 

.33.9 

541.8 

32; 8 

529.3 

20.3 

-1.06 

34.54 

33.48 

509.5 

6.5 


234.0 

234.6 

0.6 

234., 8 

B 

234.4 

0.4 

. 0.20 

0.36 

0 . 56 

234.0 

0 


• 134.8 

146.5 

11.7 

146.1 

■ 11.3 

141.3 

6.5 

-0.30 ■ 

12.97 

12.67 

133.8 

-1.0 

! 


1 , 


All temperatures in °R 





































































































































vfould still be capable of providing satisfactory results. To obtain 
some answers to these questions , additional computer experimentation 
was carried out in v/hich the hypothetical spacecraft was subjected 
to direct sun's irradiation at normal incidence for 200 minutes, fol- 
lowed by its 'entry into the shade for 100 minutes and the cycle re- - 
peated itself thereafter. The internal power dissipation was left 
unaltered at all times. At the very beginning, the spacecraft had 
a uniform temperature of 400°R throughout and the errors in the mod- 
eling parame'ters were taken as those listed in TABLE 6.4. Other con- 
ditions were the same as those used in the earlier’ studies . The. com- 
putation was continued for five complete cycles. Because of the mas- 
sive data generated, only temperatures at 30 and 200 minutes from 
the start of each cycle and those at the end of each cycle are pre- 
sented herein. In addition, data for the third and fourth cycle are 
omitted in the interest of conserving space. This is illustrated 
in Fig. 6.1, Details of the temperature data are given in TABLES 
6.7(a)-{l). It goes without saying that the results are very encour- 
aging indeed. A second experiment of the same kind was also conducted 
in which the heating and cooling time viere doubled; so vras the period 
of- the cycle. The results were equally good. Details may be supplied 
upon request, 

N. 

The foregoing studies have established, beyond doubt, the use- 
fulness of the proposed concept of imperfect modeling v/hen the errors 
in the modeling parameters are restricted such that the linear theory 
is applicable and i/hen they are properly controlled, as pointed out in 


.114 . 



remperature 



figure 6.1 Cyclic heating and cooling of the hypothetical spacecraft (open cii-cle o 
denotes instant for which data are presented in TABLE 6.7) 
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TABLE 6.7(a) Model temperatures before and after correction 

(Spacecraft subject to cyclic heating and cooling) 


30 minutes from start of the first cycle 


'Surface 

Prototype 
or Perfect 
Model 

Temperature 

T_ 

Temperatures of Imperfect Models ■ 


1 


Model Temperature 
after Correction 

Model MI 

Model Mil ' 

Model Mill 

T 

T 

T/il 

T 

MI I 


T 

MI 1 1 

E 

Ivit r I 

T 

MC 


^ ■ 

590.7 


-0,1 

590.5 

-0.2 

590.6 

-0.1 

■-0.03 

-0.12 . 

BQi 

590.7 

0 


594.9 

594.0 

-0.9 

593 .,8 

f 

594.3 

-0.6 

-0.17 

-0.70 

CO 

o 

1 

594.9 

0 


474,5 

478.0 

3.5 

470.2 

|||M| 

473.9 

-0.6 

-7.19 ■ 

10.52 

3. '33 

,474.7 

i 

0.2 '1 


590.5 

590.5 

0 

590.5 

, 0 

590.5 

0 - 

~0 

-0.04 

-o;o4 

590.5 

.0 


593.9' 

593 . 6 

-0.3 

593,6 

-0.3 

593.-7 

-0.2 

~0 

n 

-0,24 

593,9 

0 

A.' 

momiii 

451. .6 

11.3 

451,4 

11.1 

446.4 

6.1 

■-0.24 

■ 13.27 

.13.03 

438.6 

-1.7 


• 413.1 

413.4 

0.3 

411.8 

|HHH 

•412.8 

-0.3 

-1.49 

1.60 

0.11 

413,3 

0.2 


450.6 

466.2 

15.6 

465.1 

14.5 

459.9 

9.3 

. 

-0.99 

15. 88' 

14.89 

451.3 

0.7 

^ . 

' 4.45.4 

452.5 

7,1 

450.6 

5.2 

448.4 

3.0 

-1.73 

10.25 

8.53 

443.9 

'-1.5 

^o 

424.3 

426.1 

■1.8' 

miQQi 

-1.5 

424.3 

0 , 

-3.00 

. 4.43 

MQI 

424.6 

' 0.3 


547.6 

547.6 

0 

547.6 

0 

547.6 

0 

0.01- 

"0 

0.01 

547.5 

0 

> 

11 

o 

172.2 

172.8 

0.6 

172.6 

0.4 

172.5 

0.3 

-0.24 

0.98 

0.74 

172.1 

-0.1 


All temperatures in °R. 
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TABLE 6.7(b) Model temperatures before and after correction 


200 minutes from start of the first cycle 


■ 

Surface 

1 

Prototype 
or Perfect 
Model 

Temoerature 
* T, 

Temperatures of Imperfect Models 


1 

K 

Model Temperature 
after Correction 

• 'Model MI 

~ Model Mil 

Model Mill 

— 

rn 

E 

T 

MI I 


T 

1 1 

E 

!vE 1 1 

m 


A, ■ 

597.9 

597.8 

-0.1 

597.8 

-0.1 

597.9 

0 

-0..O1 

-0.08 

'-0.09 

597.9 

■ ° .1 


604.6 

604.1 

-0.5 

604.0 

-0.6 

604.. 3 

-0.3 

-o.io 

-0.46 

-0.56 

604.7 

0.1 

A. 

w> 

531.3 

561.5 

30.2 

559.2 

27.9 

548.3 

17.0 

-2.07 

33.31 

31.24 

530.2 

i 

-1.1 

1 

A, 

597.8 

597.8 

0 

597.8 

0 

597.8 

0 

-0..01 

-0.04 

-O.O'S 

597,8 

i 

0 1 


604-0 

603.8 

-0.2 

603.7 

-0.3' 

603.9 

-0.1 

-0.08 

-0.22 . 

■ 

-0.30 

604.1 

0,1 j 

A 

520.3 

552.3 

31.9 

550.3 

30.0 

538.3 

18.0 

-1.78 

35,22 . 

33.44 

518.8 

i 

-1.5 j 

A„ 

/ 


544.8 

24.9 

540. Q 

20.1 

534,0 

14.1 

-4.42 

27.53 

23.11 

521..7 

t 

1.8 

: 

519.8 

552.4 

32.6 

550.2 

30.4 

539.6 

19.8 . 

-2 . ca 

32-43 

30.42 

522.0 

2', 2 

N 

540.4 

566.8 

26.4 

562.5. 

22.1 

554.6 

14.2 

n 

30.77 


539.9 

-0..5 

A, 

521.5 

547.7 

26.2 

544.0 

* 

22.5 

536.5 

15'. 0 

-3.46 ' 



522.7 

1.2 

A: 

568.1 

568.2 

■ 

0.1 

568.2 

0.1 

568.2 

0.1 

0.02 

0.03 

0:05 

568.1 

0 



188.4 

Di 

187.9 

3.7 

186.4 

2.2 

-0.44 

5.14 

4.70 

183.7 

ran 


All temperatures in °R. 
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TABLE 6.7(c) Model temperatures before and after correction 


End of first cycle 


• 

' 

Surface 

1 

Prototype 
or Perfect 
Model 

Temperature 



• • J 

Tempo rat lire s of 

Imperfect Models | ‘ 


] *. 


Model Temperature | 
after Correction | 

Model MI 

B 




• 

i‘-> ■ 

.-■I' 

*T* 


T 

.■vs I 


I. 

T 

MI I I 

^i'vt 1 1 

m 

^^'iC 

Bill 

^ ■ 

274.7 

274.4 

-0.3 

274.4 

1 


274.5 

-0.2 

-0.03 

1 


-0.41 

274.8 

0.1 

i 2 

304.2 

302.8 

-1.4 


gH 

303.5 

-O'. 7 

-0.14 


-1.89 

304.7 

0.5 

i K 

503 . 3 

542.3 

39.0 

Era 

38.4 

526.2 

22.9 

-0.56 

i40.85- 

i 

40.29 

502.0 

-1.3 1 



274.4 

273.9 

-0.5 

■BMHI 

-0.6 

274.1 

omn 



-0.62 

274.5 

0.1 1 


302.6 ■ 

300.4 

QH 

300.2 ■ 
1 

I 

BBI 

301.5 




-2.88 '' 

303.3 

■.'0.7 j 

i 

496.4 

537.-1 

40.7 

mm 

40.0 

520.2 

23.8 


IBEfl 

42.06 • 

495.0 

~ j 

-1.4 1 

' . 

514.1- 

545.5- 

31.4 

Bra 

29.5 

■532.8 

18.7 

-1.75 

j32.13 

J 

30.38 

515.1 

1 

1.0 1 


513.6 

549.7 

36.1 


34.8 

535.8 

22,2 



34.11 

515.6 

•2.0 


530.2 

562.1 

31.9 

rai 

29.0 

548.3 

18.1 

-2.66 

35.09 
! 1 - 

32.43 

529.7 

-0.-5 


513.0 

544.6 

31.6 

HQraj 

MHIBI 

29.7 

531.8 

18.8 

-1.73 


30.81. 

513.8 

■ 

0.8 


289.7 

290.0 

0.3 


0.4 

289.9 

0.2 

j 

0,09 

0.18 

f 

0.27 

289.7 

0 


138.3 

J 

148.8 

10.5 

148.4 

■10.1 

144,1 

5.8 

rO.45 

11.95 

11.50 

137.3 

-1.0 


All temperatures in °R. 
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TABLE 6.7(d) Model temperatures before and after correction 


30 minutes from start of the second cycle 


{ Surface 

1 

i 

Prototype 
or Perfect 
Model 

Temperature 

■ ■ ■* 

Temperatures of Imperfect Models’ 

* 

A ^’^' 

1 

• •• . ' 

Model Temperature 
after Correction 

•Model MI 

Model Mil 

Model Mill 

( 

M 

T 

.MI I 


T 

MI 1 1 

I Mr 1 1 

fr» 

\ v.c 


! ^ 

y 

580*6 

. 

580.5 

- 0.1 

580.5 

- 0.1 

580.5 

- 0.1 ■ 

- 0.01 

- 0.09 

- 0.10 


0 

■■ A , . 

590.5 

590.1 

- 0.4 

590.0 

- 0.5 

590.3 

- 0.2 

- 0.08 

- 0.40 

- 0.48 

590.6 

0.1 


518.3 

552.3 

33.9 

550,5 

32.2 

537.4 

19.1 

- 1.54 

37.56 

36 . 02 

516 '. 2 

- 2.1 1 


580.4 

580.4 

0 

580.3 

- 0.1 

580.4 

0 ■ 

•- o '. 01 

- 0.09 

• - 0.10 

580.4 

“ 1 

• K 

589.6 

589.5 

- 0.1 

589.5 

- 0.1 

583.6 

0 


- 0.30 

mm 

589.9 


\ ■ 

• 504.5 

543.7 

39.2 

543.1 

38.6 

527.4 


- 0,56 

• 41.28 

40.72 

503.0 

- 1.5 I 

! 


' 513.5 

545.4 

32.0 

543.6 

30.2 

532.5 

19.0 

- 1.66 

32.84 

• 31 . 18 

514 ; 3 

0,8 ( 


513.9 

550.2 

36.3 

549:0 

35.1 

536.2 

22.3 

- 1 , 08 . 




. 1.8 


531.1 

563.3 

32.2 

_ 560.5 

29.4 

549,3 

18.2 

72.58 

35.44 

32.86 

530.4 

- 0.7 

^0 

513.4 

545.3 

31.9 

543.4 

30.0 

532.3 

18 '. 9 

- 1.75 

33.01 


514.0 

0.6 


521.2 

521.4 


521 . 4 . 

0.2 

521.3 

0.1 

0.05 

0.11 



0 


176,4 

, 

181.8 

wm 

181.5 

wm 

179 . 3 ' 

2.9 

- 0 . 25 - 

6.27 

6.02 

175.8 

- 0.6 

1 


All temperature in °R, 
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TABLE 6.7(e) Mcidel temperatures before and after correction 


200 TEiinutes from start of the second cycle 


Surface 

i 

Prototype 
or Perfect 
Model 

Temperature 

j Temperatures oi 

1 1 1 ,1 

Imoerfect Models ■ 

* 

f 

' 1 

i . i' 

A 

j 

j 

Model Temperature 
after Correction 

Model MI 

Mode 

1 Mil ' 

1 ’ 
Model Mill 

I ■' '! 

1 

m 

Nil 

*M[ I 


1 - 

T 

M I z 

\ 

^ivi r I 

m 

\v.c 

■ j 

1 ' 

598.0 

597.9 

-0.1 

597. 9| 
1 

-0.1 

597.9 

-0.1 

•-0 . 01 


-0.08 

598.0 

0 

i ^ 

604.7 

604.4 

-0.3 

mi 

-0.4 1 604.5 

HB 

-0.04 

f 

= -0 . 41f 
1 

-0.45 

604.8 

. 0.1 1 

1 ^ 

534.0 

566.5 

32.5 

n 

31.6 

552.7 

18 .-7 

-0.83 

‘ 35 . 00 ; 

1 • 

34,17 

532.3 

i 

-1.7 ,j 


597.9 

597.8 

-0.1 

597. Sj 
1 

-0.1 

597.9 

0 

-0.01 

imHi 

-0.04 

597.9 

! 

0 f 


604.2 ■ 

604.1 

-0.1 

604. li 
i 

-0.1 

604.2 

_0 

-0.03 

IDB 

-0.22 

604.3 

.0.1 


524.6 

559.2 

34.6 

HI 

33.9 

544.6 

20.0 

-0.70 

■37.14; 
1 1 

36.44 

522.8 

r-1.8 j 

< 

\ 

528.1 

558.1 

30.0 

556.5! 

28.4 

545.9 

17.8 

■m 

-30.71 

1 

29.27 

528.'8 

! 

0.7 

% 

524.2 


3-5.1 

558:4 1 
1 

34.2 

545.9 

.21.7 

-0.85 

34.0.5' 

• i 

33.20 

1 

526.1 

1.9 


546,5 

576,0 

29.5 

573-, 4.' 
1 

26 1 9 

1 

563.1 

n 

-2.41 

32.66 

30.25 

545.7 

-0.8 1 
1 

^0 

527.9 

557.8 

29.9 


28.3 

545.7 

17.8 

-1.44 

IBI 

29.19 

528.6 

1 

0.7 I 


568.1 

568.2 

0.1 

568.2. 

0.1 

568.2 

0.1 

0.02 

0.04 

■ 0.06 

568.1 

1 

i 


184.9 

189.7 

4.8 

189.5 

m 

187.5 

2.6 

-0.20 

5.63 

5.43 

184.3 

j 

-0.6 ! 
( 


All temperatures in °R. 
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TABLE 6.7(f) Model temperatures before and after correction 


End of second cycle 


Surface 

! 

1 

Prototype 
or Perfect 
Model 

Temoarature ; 


Temperatures of Imperfect Models 


' A<->- 

1 


.Model Temperature- 
after Correction 

•Model MI 

Model Mil ' 

Model Hill 

m 


T 

•w I 


T 

Nil 1 1 

1 1 


*^r/c 

^ ■ 


274.5 

-0.3 

274.5 

-0.3 

274.6 

-0.2 

•~0 

-0.36 

-0.36 

274.9 

0.1 


304.4 

303.2 

mgm 

303.3 

-1.1 

303.9 

-0.5 

0.01 

-1.66 

-1.65 

304.9 

0.5 


503.9 

543.5 

39.6 

543.3 

39.4 

527.2 

23,3 

- 0.11 . 

41.16 

41 ..05 

502.4 

- 1.5 1 

*[ 


274.5 

274.0 

- 0.5 

274.0 

- 0.5 

274.2 

-0.3 

- 0.01 

-0,57 

-0.-58 

274 . 6 

mgm 


303.0 



300.9 

- 2.1 

302.0 

-1.0 

-0-06 

-2.62 


303:7 

1 

0.7 . i 

. 

497.3 

538;6 

41.3 

538.4 

4i.i 

521.6 

24; 3 

.-0.16 

43.00 

42.84 

495.8 

jggm 

A 

• 

'515.7 

548.2 

32.5 

547.3 

31.6 

535.3 

' 19.6 

-0.78 

32.66 

31.88 

516:3 

0.5 ’1 


' 514.5 

551.1 

36.6 

550.3 

35.8 

537 .1 

22.6 

-0.73 

35.57 

34,84 

516:3 

, 1.8 


5.31.4 

• 

564.1 

32.7 

561.8 

mm 

550.1 

18.7 

-2.0’7 

35.41 

33 . 34 

530,7 

jBm 

■ 

514.2 

546.7 

.32.5 

545.6 

31.4 

533.7 

19.5 

-1.02 

32.90 

31.88 

514.8 

0.6 

•^T 

289.8 

290.0 

0.2 

290.1 ' 

0.3 

290.0 

0.2 

0.10 

0.18 

0.28 

289.8 

0 


, 138.6 



149.4 

10.8 

149.1 

• 10.5 

144.6 

6.0 

-0.28 

12.07 

11.79 

137.6 

-1.0 


All temperatures in °R 
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TABLE 6.7(g) Model temperatures before and after correction 


30 minutes from start 

of the fifth cycle 

i :• 




Temperatures o, 

' Imperfect 

Models 

i 

i ’ 

• , 

t 


' 


Surface 

Prototype 

or Perfect 

Temperature 

T 

p 

•Model MI 

Mod 

5l 

Mil 

Model mItII 

■■ i 

i 


i , * * • 


1 

Model Temperature 
after Correction 

m 

r» 


1 

I. 

T 

M I r 


, 

'* 

^ 1 
i •, : ’ 

rrr 

^r/c 

A 

580.6 

580.5 

-0.1 


-0,1 

580.5 

tO.I 

1 

-0.01 


-0.09, 

■580.6 

0 

/V 

590.5 

590.2 

-0.-3 

590. i 

-0.4 

590.3 

-0.2 

! 

-0.04 

: -0.38 

i 

-0.42 

590.6 

0.1 . 


518.6 


34.4 


33.0 

538.1 

19,5 

-1.22 . 

. 37,74 

j 

36 . 52 

516:4 

-2.2 i 

j 

A ■ 

580.4 

580.4 

0 

n 

0 

580.4 

■ 0 

~0 

, -0.08 
J 

-0..08 

580.5 

0,1 j 


589.7 

589.6 

-0.1 

■1 

-0.1 

589.8 

n 

-0.01 

Mlfl 


589.9 

. 0.2 ! 


505,1 

' 

544.7 

39.6 ' 

ES 9 

39.3 

528.3 

23.2 


41.48 

41.25 

503.4 

■-1.7 

< 

1 

514.5 

547.2 

32.7 

Qim 

31.6 

5 . 34.1 

19.6 

-0.96 


32.19 

515.0 

0.5 i 


514.4 

551.2 

36.8 

550:3 

35.9 

537.0 

. 22.6 

-0.76 

35.77 . 

35.01 

516.1 

• 1.7 

'V 

531.8 

564.5 

32.7 

562.2 

30.4 

550.5 

18.7 

-2.14 

35.61 

33.47 

531.1 

-0.7' 


514.2 

: 

546.7 

32 . '5 

545.3 

> 

31.1 

533.5 

19.3 

• 

-1.23 


31.99 

514.7 

0.5 


521.2 

521.4 

0.2 

521,4 

1 

1 

0.2 

521.3 

B 

0.06 


0.17 

521.2 

0 


176.5 

• 

182.0 

5.5 


5.3 

179.5 

3.0 

-0.18. 

6.33 . 

6.15 

175.8 

-0.7 . i 
1 


All temperatures in °R. 
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TABLE 6.7(h) Model 


200 minutes from start of the fifth cycle 


Surface 

1 

. 

■ 

Prototype 
or Perfect 
Model 

Temperature 

T 

Temperatures of Imperfect 

Model MI 

Model Mil 

m 

.'/2 

T 

.MI I 



HI^^qIIIIIIII 

597.9 

-0.1 

597.9 

-0.1 


604,7 

604.4 

-0.3 

604.3 

-0.4 



566.5 

32.5 

565.7 

CO 


597 . 9 

597.8 

-0.1 

597.8 

-0.1 


604.2 

’ 

604.1 

-0.1 

604.1 

-0.1 

i 


559.3 

34.7 

558,6 

34.0 

'V 

523.1 

558,2 

30.1 

556.7 

28.6 


524.3 

559,4 

35,1 

558.5 

34.2 

N 

546.5 

576.1 

29.6 

573,5 

27.0 • 

^o 

528.0 

557.9 

29.9 

556.4 

26.4 


568.1 

568.2 

.0.1 

568,2, ’ 

0.1 

1 

■ 184.9 

189.7 

4.8 

189.5 

4.6 


All temperatures in °R 


atures before and after correction 


Models 


Model Mill 



604.5 


557.9 


E 

Ivt n 

1 -0 

.1 

■ 

.2 

18 

.1 

0 

0 

20 

.0 

00 

.0 

21 

. 6 

16 

,7 

17 

.8 

i ° 

.1 

2 

,6 


« 

a'-> 

* 

** ' 

Model Temperature 
■ after Correction 

rrl 

■‘ViC 


-0.01 

-0,07 , 

-0.08 

598.0 

0 

-0.04 

-0.41 

-0.45 

604.8 

0.1 

-0.79 ■ 

35.01 

34.22- 

532.3 

-1.7 

~0 

-0.03 

-0.03 

597.9 

0 

*1 

O 

o 

CO 

■-0.19 • 

-0.22 

604.3 

0.1 

-0.66 

37.15 

36.49 

522.8 

-1.8 

-1.36 

30.71. 

29.35 

528.8 

0.7 

-0.81, 

34.05 

33.24 

526.1 

1.8 

1 

1 

i 

to 

00 

<M 

1 

32.66 

30 . 3,1 

545.8 

-0.7' 

-1.38 

30.63 

29.25 

528.6 

0.6 

0.02 

IHRMl 

0.06 

568.1 

0 

-0.66 

37.15 ■ 

36.49 

184.3 

-0.6 
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I • \ , 

j TABLE 6.7(i) Model temperatures before and after correction 

End. of the fifth cycle • i . . 


I ' I 


Surface 

1 

1 

Prototype 
or Perfect 
Model 

Temperature 

T 

p 

. 

Temperatures of Imperfect Models j 

, 

• 

t 

. 

! j ■ ’ 

' a'->- 
1 1 

! • "i • 


Model Temperature 
after Correbtion 

■ Model MI 

.Model Mil 

i 

Model Mill 

« t 

' t 

m 

‘J/t 

iiZ 

T ■' ■ 
.VI I, 


T 

1 1 

i 

E 

j iva 1 1 

m 

• \v,c 

Bi 

i 

274.8 

274.5 

-0.3 

274.5 

-0.3 

274.6 

mm 

~o- 

5 \ 

i -0.36 

-0.36 

274.9 

.. 0..1 

! ^ 

304.4 

303.. 2 


303.3 

-1.1 

303.9 

n 

, 0.01 

! ] 

! - 0.17 

-0,16 

304.9 

0.5 

! ■ ^ 

503.9 

543.5 

39.6 

543.4 

39.5 

527.2 

2^.3 

} 

-0.10 ■ 

; 

41.07 

502.4 

i 

-1.5 '! 


274.5 

274.0 

-0.5 

274.0 

-0.5 

274.2 

DB 

-0.01 

■ - 0.57 

-0.58 


1 

0.1 1 

■ 

303.0 • 
' 

301.0 

-2.0 

300.9 

mm 

302.0 

-1.0. 

mm 

! -2,62 
1 

-2.68 

303.7 

0.7 


497.3 

538'6 

41.3 

538.5 

41.2' 

521.6 

24.3 

■-0.15 

wm 

42.85 

-_4.9.5..,8 


' . \ 

515.7 

548.2 

32.5 

547.4 

31.7 

535.3 

1§.6 

-0,77 

1 32.66 

31.89 

HHi 

0.6 

\ ■ 

514.5 

551.1 

36.6 

550.4 

35.9 

537.1 

WM 

-0.72 

; 35.57' 

34.85 

516.3 

'l .8 

^ . 

. 5.31.4 

564.1 

32.7 

561. 8 

HDI 

550.1 

18.7 

Hymniii 

1 35, 4d 

33.34 

530.7 

-0.7 

^0 

514.2 

546.7 

■32.5 

545.6 

31.4 

533.7 

19.5 . 

-1.01 

32.90 

31.89 

514.8 

0.6 

. 

289.8 


.0.2 

290.1 ' 

0.3 

290.0 

6.2 

O'.'lO 

0.18 

0.28 

289.8 

6 

1 

138.6' ‘ 

149.4 

10.8 

149.1 

' 

10.5 

144.6 

6.0 

-0.27 

12.07 

11.80 

137.6 

-1,0 


All temperatures in °R 








































































































































Chapter 2. During the course of the present investigation, additional 
observations, some of which are quite essential for the successful 
application of the theoiy, have been made. They are discussed briefly 
in the section V7hlch follows . 

.6.3 SOME RELEVANT OBSERVATIONS 

6.3.1 Influence of Decreasing and Increasing Model Errors on 
Prediction Reliability 

The errors in the various modeling parameters as listed 
in TABLES 6.1 and 6.4 are of the order of 20 to 25 percent in kd and 
Cd and of 10 to 20 percent in e and e ' for the ' original ' model MI . 

It V70uld be of interest to know how the reliability of the final cor- 
rected temperature data be altered by increasing and by decreasing 
such errors. To this end, the computations described in Section 6.1 
were repeated, first with all errors reduced to one-half of their 
.respective values as shown in the tables. It was found that the im- 
provement in the accuracy of the corrected model temperatures of 
was totally insignificant. In no case did it exceed 1°R. Thus, if 
the present theory is used for the correction of data gathered v?ith 
imperfect models, there is no real need of striving for the undue 
reduction of modeling errors , provided that they are kept vrithin 
limits . This welcoming fact allows additional flexibility in model 
fabrication. On the other hand, •sihen errors in the modeling parame- 
ters are doubled, the reliability of the correction procedure dete- 
riorates rapidly. Errors in as large as 20°R have been noted. 
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This points to the need of a development of a more general theory 
-which could accommodate greater errors in the modeling parameters. 


6,3.2 Proper Control of Model Errors and Selection of Experi- 
mental Conditions 


—liT'has been "demonstrated "in 'Chapter 2 that the best re- 
sults obtainable from the present theory, arise when the ratios 

and ^ ^ close to 1/2 Csee Fig. 2.3(a.l) and (b.2)] 

or 2 [see Fig. 2.3(a.2) and (b.l)], depending on x-rhether, in the modi- 
fication experiments, the magn.'itude of the 6's increases or decreases. 
These ratios, evaluated for the simple error space, are, respectively, 
0.474 and 0.532 for models prescribed in TABLE 6,1 and are 0.474 and 
1.64 for models prescribed in TABLE 6.4. The ranges of errors shown 
in these tables are, in fact, selected with the guidance of the said 
theoretical requirements . ’ The question which naturally arises is ; 

In the modification experiments with models Mil and Mill, would it 
not be desirable to 'strive for smaller errors? That is, to make the 
positive errors become less positive and negative errors less nega- 
tive, t To see if this were the case, the computations described un- 
der Section 6.1 and pertaining to models prescribed in TABLE 6.4 v;ere 
repeated, but with both positive errors in model Mil and negative 


'errors in Mill kept less than 5 percent. The change gave rise to 
rather unfavorable ratios in namely, ~ 0.31 and 




0.22. It was found that the computed data forl^^ shoxved 


tSuch modification scheme may be unrealistic in practice. However, • 
this is HOT the point in question here. 
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littl'e or no improvement over those given in TABLES 6.5 and 6.6. 

In the discussion of error paths in the multi-dimensional error 
space, four possibilities are .considered as illustrated in Fig. 2.3. 

In the first modification experiment, all positive errors remain posi- 
tive after modification, X'lhether they become smaller or larger. Like- 
wise, in the second modification experiment, all negative errors re- 

main. .negative . -Clearly, such are -not t-he -only -possibilities ; -it -is 

physically feasible to have mixed errors (positive and negative) af- 
ter modification, although they should always be avoided for reasons 
to be seen shortl 5 '-. Consider, for instance, the experimentation with 
Mill. Some of th.e..,ori.ginally negative errors may .become positive 
after modification while others remain negative. Under such circum- 
stance, the parabola joining the error states P, ^ and 0^ ^ in 
subspace ^ would have a relatively sharp curvature, thus creating 
a condition that is in direct conflict with the fundamental assump- 
tion upon which the theory is built. 

To provide a pictorial illustration, we consider a three-dimensional 
simple error space shown in Fig. 6,2. The coordinates of the point 
P which is associated with the original model MI are all negative. 

Suppose that, after modification, 6* and 6* remain negative, but 6* 

lb ^ 

become positive. Thus, the point P*^ ^ vrould be as shovm in Fig. 6.2(a). 
The error path joining P, P* and 0 makes a sharp bend as' illus- 
trated and it would totally invalidate the linear theory developed in 
[Ij and in this report. On the other hand, if 6- remains negative 
as sho;-m in Fig. 6.2(b), the error path vTould have a gentle curvature 
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which is small everywhere along the path. If all 6*’s are made posi- 
tive, it is again possible to have an error path v^hich has a uniformly 
small curvature. This is illustrated in Fig. 6.2(c). The three points 
through which the parabola passes must novr be of the sequence P, 0^ ^ , 
and . The extension of the foregoing consideration to the design 

of model Mil experimentation is obvious. 

, _ .Extensive numerical calculations have been carried out, not only 
for the hypothetical spacecraft considered in the present report, 
but also -for a totally different problem, namely, the correlation 
of boiling heat transfer data, and the results conclusively demonstrate 
the validity of the arguments just given. The nature of the error 
paths has a controlling Influence on the reliability of the theory. 
Error paths exhibiting small and sharply changing radius of curvature, 
such as that of Fig. 6.2(a), should not be considered in any circum- 
stance. . 

6.4 CONCLUSIONS AND RECOMMENDATION 

Based on the evidence provided by a computer experimentation 

of the application of the proposed new theory of imperfect modeling 

to the study of the thermal performance of a hypothetical spacecraft, 

the follovfing conclusions may be drawn. The hypothetical spacecraft 
\ 

has major radiative and conductive heat flow paths that crudely simu- 
late those of the '64 Mariner family of space vehicles. 

(1) Vhien errors in kd and Cd of up to 25 percent and in surface 
emittance of up to 20 percent exist in the model and when 
the spacecraft is subject to either a simple heating and 
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cooling transient with an intervening steady condition or 
when, it experiences a continuous cyclic transient, the in- 
dicated model temperatures would entail errors of 10 , 20 , 
and 30®R.and higher. These errors were reduced to less 
then 3°R after a correction according to the theory. There 
is no evidence of error accumulation under the cyclic con- 
dition investigated. 

(2) VJhile the theory is deduced under the basic assumption that 
the errors in the modeling parameters are small and that 

. they have independent effects on the dependent variable 
or variables , there is no need to strive for the undue , 
costly reduction of these errors. Rather, in the modifi- 
cation experiments, the model design and/or the test con- 
ditions should be properly controlled in order that the 
ratios /5p*^ do not deviate excessively 

from 1/2 when the magnitudes of the 6's generally increase 
or 2 when they generally decrease. 

(3) The reliability of the theory deteriorates rapidly when 
errors in the modeling parameters become large. At the 
present time , it is not feasible to describe quantitatively 
what is meant by 'large' or 'small' error. 

In view of the success, as well as the limitations found in connec- 
tion with the performance of the theory, it is recommended that fur- 
ther investigation be launched with emphasis on developing a means 
of treating large errors. 
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APPENDIX A 


THERMAL RESISTANCE OF MULTI-LAYER INSULATION 


The superinsulation used in today’s spacecraft is a multi-layer 

aggregate of aluminized mylar or teflon. In the ideal condition, - 

the layers form parallel planes and have no physical contact with 

each other. For the purpose of our present analysis, such condition 

is assumed. Here again, we adopt a semi-gray, two-region spectral 

subdivision of properties; namely, the solar and the infrared region. 

The directional properties of the surface are taken into account 

by the simple diffuse-specular model; i.e.,p=p +p. All emit- 

d s 

ted radiant energies are assumed diffuse. 




T 


(Temp, of the Base 


Plate) 


Fig. A.l A Base Plate Covered with Superinsulation 

Figure A.l depicts a base plate at temperature T^ covered with 
a superinsulation of n-layers . The exposed surface is plastic and 


.132 



is absorbing sun^s radiation at a rate; 


0 *-^ t 

~ = a^^'S cos 0 

A p s 


(A.l) 



solar constant (at the outer fringes of earth's atmosphere, S = 442 
2 

Btu/hr-ft. ). All properties shown in Fig. A.l without an asterisk 
are for the infrared range and the siibscript 'a' refers to aluminum 
and 'p' refers to plastic. Appropriate values of reflectances and 
absorptances of aluminized mylar or teflon currently used in sp'ace- 
craft applications are listed in Table A.l. 

TABLE A.l 

REFLECTANCES AND ABSORPTANCES OF ALUMINIZED MYLARt 



Solar Range (0.25p-2.5y) 

Infrared Range 

(> 5u) 


p* p* a* 

Ps 

Pd 

a(= e) 

Aluminized side 

0.85 0 0.15 

0.95 

0.00 

0.05 

Plastic side 

0.85 0 0.15 

0.00 

0.50 

0.50 


Consider a typical section of two adjacent layers of the super- 
insulation; e.g., A^ and A^ of Fig. A.l. V?e assume that the tempera- 
ture of each layer is uniform and that the edge effects are negli- 
gible. Figure A. 2 shows the equivalent resistance netvfork when steady 
condition prevails . In the figure , is the black body emissive 
power, B denotes the diffuse radiosity and E is the exchange factor. 

A O 

1'Data transmitted to the author by Mr. W. A. Hagemeyer of JPL. 


133 




The Slabs cpipts 2 and 3 have the usual meaning. 




B- 


-b,2 

o — -vwwvv- 




1-0 

Os,p 


b,3 


■A^WWVV 

* 


— AAAWsAr- 

'j.p 


s,p 


A2E,3U-£.)(I-p_) 


F1g. A. 2 Radiation Network for Infinite Parallel Planes 


since A = A = A, the thermal 'resistance' per unit area is 
^ <6 


R = 


^b, 3. ^b, 2 

Q/A 


d, a 


d, p 


e (1 - P ) E (1 - p ) 

a s , a p s , p 


E_ (1 - P )(1 - P ) 

■ 23 s ,.a s , p 


(A. 2) 


The exchange factor E can be easily evaluated using the image tech- 
nique. Recognizing that the separation distance between layers is 
very small as compared with the lateral extent of the superinsula- 
tion, one finds 


E„„ = 1 + p P + 


23 


s,a s,p s,a s,p 


1 - P P 

s , a 5 , p 


(A. 3) 


Hence , 


R 


d, a 


d, p 


1 - P P 

s , a s , p 


£(1-P.) E(1-P ) Tir^^-p )(1 - P ) 


P 


s , a 


s, p 


(A. 4) 
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If the surfaces are diffuse, p = p = and Eq. (A, 4) simpli- 

s , a s I P 

fies to 



..which is a vrell-known result . 

The resistance per unit area R' between the alurainuni surface 
of the last layer of the superinsulation and the protected surface 
of the base plate can also be evaluated from Eq. (A. 4), provided 
'that the three properties associated with the plastic, namely, , 
p^ p and p^ ^ , and replaced by the corresponding properties of the 
base plate, e. , p, , and p , . 

The inward 'leakage* heat flux through the superinsulation ^ 
can be determined from a consideration of the heat balance at its 
outermost exposed surface. It is 


— - OP'S cos 0 - £ pT^ 

A p P exp 


4 - 4 

CTT - C7T 

exp b 



(A. 5a) 


with 



(n - 1)R + R’ 


(A. 5b) 


n being the total number of layers in the superinsulation. In prac- 
tice, for the sunlit surfaces of the spacecraft. 
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(A. 6) 



a*S cos 6 
P s 


or 


e oT 

p exp 


and, hence, a valid approximation for Eq. (A. 5a) is 


Q 

A 



ct* 

p 

£ 

P 


S COS 0 

s 



(A. 7) 


which may be either positive or negative. When the superinsulation 
sees only the empty space, any possible heat leak is away from the 
surface. A straightforward calculation leads to 



p 


(A. 8) 


where ^ R is again given by Eq. (A. 5b). 

While the foregoing results were obtained under the assumption 
of steady heat flow, they may be used for transient analysis with- 
out entailing significant errors since the heat capacity of multi- 
layer insulation is usually quite small and, in spacecraft applica- 
tion, the transients are seldom very rapid. 

^ To ascertain the effectiveness of of superinsulation in reduc- 
ing surface heat flux, we consider an aluminum base plate, with and 
without the superinsulation , vrhen it is exposed to 

2 

a. sun's irradiation of a local intensity of 300 Btu/hr-ft 
and 
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1). the black empty space. 

The plate has a temperature of 528°R and is either highly polished 
or is sprayed with PV-100 white paint, Bv using Eq. (A. 4) and the 
relevant data in TABLE A.l, one finds 

■ R =21 for aluminized mylar or teflon , plastic side out 

38.8 when the aluminum base plate is highly polished 

R’ = J 

I 20.2 when the base plate is covered with PV-100 
^ white paint 

Then, it follows from Eq. (A. 5b) that, for a five-layer superinsu- 
lation, ~ 122.8 and 104.2, respectively, for the highly polished 

and painted base plate. The corresponding values. for a’ ten-layer 
Insulation are 227.8 and 209.2; and for a twenty-layer Insulation, 
437.8 and 419.2. 

With the foregoing information, the leakage surface flux can 
be readily calculated from Eq, (A. 7) and Eq. (A. 8). The results 
are shown in Table A. 2. In practice, the effectiveness of insula- 
tion is somewhat less due to the unavoidable conduction leaks. 
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TABLE A. 2 


Effect of Multi-Layer Insulation on Surface Heat Flux, Btu/hr-ft^ , 
of an Aluminum Base Plate at 528°R 



Based Plate Highly Polished 

Base Plate Covered with PV-100 White Paint 


Without 

Superinsulation 

With Superinsulation 

Without 

Superinsulation 

With Superinsulation 


(number of lay 

ers ) 

(number of layers) 



5 

10 

20 



Surface Exposed 
to Sun's Irradiation 

53.3 

-0.35 

-0.19 

-0.10 

- 50.6 

-0.41 -0.21 -0.10 

Surface Exposed 
to Empty Space at 0°R 

- 6.7 

-1.07 

-0.58 

I 

o 

CO 

o 

-111 

-1.25 -0.63 -0.32 


The positive value is for heat flow into the base plate; all negative values are for heat flow away from the base 
plate . 





















APPENDIX B 


A NUMERICAL EXAMPLE 

As an example, we consider the case in which errors in the mod- 
eling parameters are as shown in TABLE 6.1 and the sun's rays strike 
the solar panels at normal incidence. Specifically, we refer to data 
listed in TABLE 6.2(d) which is for a heating transient, 300 minutes 
after the spacecraft is subjected to solar radiation and internal 
heating. The choice of TABLE 6.2(d) for illustration is totally ar- 
bitrary, any other tabulated data given in the report may serve the 
purpose equally well. 

First, we note that the perfect model temperatures listed under 

T are calculated from the equation set (3.3.6) using the theoreti- 
p 

cally correct values of the modeling parameters . This is achieved 
by the combined Nevjton-Raphson and Gauss-Seidel procedure as explained 
in Chapter 5. The various temperatures of the three imperfect models 

computed from the same equation set using 
the same procedure, but with modeling parameters entailing errors 
as listed in TABLE 6.1, 

To Illustrate hov7 the various corrections are computed, it is 
sufficient to work out the details for one of the surfaces, say A . 

At t = 300 minutes after the start of the heating transient, 

\ii ” 542.1^'R, = 538, 3°R, and ^ ^ - 555. 4°R. These may be 

compared with the theoretically correct temperature of 524. 4°R and, 
hence, their respective errors are 17.7, 13,9, and 31.0°R. All 
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temperatures cited here may he read from TABLE 6.2(d) in entries fol- 
lowing A . In passing, we note that the three temperatures associ- 
6 * 

at.ed with the three imperfect models cor--respond to the function 
of Chapter 2 as follows : 


T 





T 



Mil 


T 

Mill 




DETERMINATION OF CORRECTION FOR POSITIVE ERRORS, 

From TABLE 6.1, it is seen that the various posi-t-ive errors are: 


(a) for Model MI 



A 

1 

A 

2 

A 

3 

A 

4 

A 

5 

A 

6 

A 

7 



A 

10 

^I 

^2 


— 

— 

0.25 

Bj 

B 


0.20 

0.25 

0.25 

0.20 

B 

— 


— 

— 

0,25 

B 

B 


0,25 

0.25 

0.25 

0.25 

1 

— 

6 

e 

— 

— 

— 

B 

B 


0.15 ■ 

— 

0.25 

0.15 

1 

0.15 

'e' 

— 

— 

0.10 

B 

B 


— 

— 

— 


B 

— 


(b) for Model Mil 



A 

1 

A 

2 

A 

3 

A 

4 



A 

7 


A 

9 


1- 

lO 

6* 

kd 

— 

— 

0.55 

— 

— 

0,45 

0.40 

0.575 

0.45 

0.44 

B 

■■ 

6* 

Cd 

— 

— 

0.52E 

— 

— 

0.50 

0.55 

0.575 

0.45 

0.60 

1 


6* 

e 

— 


— 

— 

— 

— 

0.27 

— 

0.45 

0.33 

1 

0.33 

6*, 

e 

— 

— 

0.20 

— 

— 

0.20 

— 

— 

— 

— 

— 

— 


140 








































The asterisk designates errors after modification. In the present 
instance, all errors become more positive after modification. Hence, 
the error path is as illustrated in Fig. 2.3(a.l) and (2.3.8) is the 
appropriate formula for determining . 

From the foregoing data, we may readily calculate 



0.25 + 

0.25 + 

0.10 + 

0.25 + 

0.20 + 

0.10 

2 

2 

2 

2 

2 


+ 0.20 

+ 0.25 

+ 0.15 

+ 0.25 

+ 0.25 

+ 0.25' 

+ 0.25^ 
0.833 

+ 0.25^ 

+ 0.20^ 

2 

+ 0.25 

+ 0.15^ 

+ 0.15‘ 


2 <6*)^ = 0.55^ + 0.525^ + 0.20^ + 0.45^ + 0.50^ + 0.20^ 

+ 0.40^ + 0.55^ + 0.27^ + 0.575^ + 0.575^ + 0.45^ 
+ 0.45^ + 0.45^ + 0.44^ + 0.60^ + 0.33^ + 0.33^ 

= 3.687 

= 0.25 X 0.55 + 0.25 X 0.525 + 0.10 X 0.20 + 0.25 
X'0.45 + 0.20 X 0.50 t 0,10 X 0.20 + 0.20 
X 0.40 + 0.25 X 0.55 + 0.15 X 0.27 + 0.25 X 0.575 
+ 0.25 X 0.575 + 0.25 X 0.45 + 0.25 X 0.45 + 0.25 
X 0.45 + 0.20 X 0.44 + 0.25. X 0.60 + 0.15 X 0.33 
+ 0.15 X 0.33 = 1.740 

Using (2.3.9a), we find 


E 6. 6?^ 

] J 
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T1 


.(+) 


p 


0.833 X 3.687 


- 1 


1.740' 

p _ 1.740 


1 / 2 


= 0.112 


j.{+) 3.687 

^p* 


0.472 


Also, - V 3,687' = 1.32. Hence, from (2.3.4a5b), we determine 

P* 


a = 0.112 X 


1 - 0.472 


= 0.212 


T- — ^ - n 1 1 n 


and, from (2.3.5c), 


X = 1 - 2 X 0.472 = 0.056 


Substituting the foregoing values of a and X in (2.3.6b) gives. 


f(a,X) = 0,0549 


and, hence, the ratio of arc lengths is 




(+) 1 - 0.0549 

• s 


= 2.116 


as according to (2.3.6a). Finally, the required correction for all 
positive errors is given by (2.3.8) 




538.3 - 542.1 
2.116 - 1 


-3.40®R 


which is listed in TABLE 6.2(d). 
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DETERMINATION OF CORRECTION FOR NEGATIVE ERRORS, 


From TABLE 6.1 , it is seen that the various negative errors are: 
(a) for Model MI 



and 6 = -0.102 for superinsulation covering A and A . 

7 10 

(b) for Model Mill 



and 6* = -0.205 for superinsulation covering A and A . Again, the 

2jK * 7 10 

asterisk denotes errors after modification. Since all errors become 
moTe negative after modification, the error path is as illustrated 
in Fig. 2.3(b.2). From the data, we find 



(-0.18)^ + (-0.20)^ + (-0.08)^ + C-0.18)^ + (-0.08)^ 
+ C-0.08)^t (-0.08)^ + 2x(-0.102)^ = 0.1512 
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2(6*f = (-0.288^ + (-0.32)^ + (-0.147)^ + (-0.259^ 

3 

+ (-0.115)^ -3- (-0,36)^ -1- (-0.128)^ 

•+ 2x(-0.205)^ = 0.412 

a 6 6* = (-0,18)(-0.288) + (-0.20)(-0.32) 

3 3 

+ (-0.08)(-0.147) + (-0.18)(-0.259) 

+ (-0.08)(-0.115) + (-0.08)(-0.16) 

+ (-0.08)(-0.128) + 2(-0.102)(-0.205) = 0.247 





0.395 


0.395 

0.642 


-0.615 


X ,= 1 - 2 X 0.601 = -0.202 
f(a,X) = -0.198 (from 2.3.6b) 

■ Vs^ = 1-669 (from 2.3.13b) 

1 T o.xyo 
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and' 


^ = 19.890R (from 2.3.13a) 

1.669 - 1 

which is given in TABLE 6.2(d). The total required correction is 

- -3.40 +19.89 = 16.49°R 
and, hence, the corrected model temperature is 
= 542.1 - 16.49 = 525. 6°R 

Upon comparing with the theoretically, correct value of 524. 4°R, one 
sees that the error in is 1.2°R. 
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